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ABSTRACT 

(Distribution  Limitation  Statement  No.  1) 


The  use  of  laser  detection  and  ranging  (ladar)  systems  in  the  atmosphere  is 
limited  by  the  severe  attenuation  of  the  signal  due  to  absorption  by  atmos¬ 
pheric  gases  and/or  scattering  by  cloud  and  fog  particles.  Furthermore,  the 
cloud  and  fog  particles  may  scatter  the  transmitted  radiation  toward  the 
receiver,  thus  producing  a  false  return,  which  must  be  distinguished  from  the 
true  target  return.  Described  herein  is  an  optical  polarization  technique  for 
discriminating  between  a  cloud  return  and  a  target  return.  Also  presented  are 
the  results  of  a  research  program  which  demonstrated  the  validity  of  the  dis¬ 
crimination  technique  both  in  the  laboratory  and  in  the  field.  The  discrimina¬ 
tion  technique  is  based  on  the  theory  that  if  a  linearly  polarized  laser  beam 
propagating  through  an  atmospheric  cloud  or  fog  is  incident  on  a  diffuse  sur¬ 
face,  then  according  to  Hie  theory  the  radiation  scattered  from  the  cloud  is 
either  linearly  or  elliptically  polarized  while  the  radiation  reflected  from 
the  surface  is  depolarized.  The  results  of  the  investigation  described  have 
shown  that  if  the  transmitter  output  of  a  ladar  system  is  linearly  polarized 
and  an  analyzer  with  its  plane  of  polarization  perpendicular  to  that  of  the 
transmitted  wave  is  placed  at  the  input  to  the  receiver  (which  is  adjacent  to 
the  transmitter),  then  the  cloud  return  will  be  up  to  94  percent  less  than  the 
value  obtained  when  the  analyzer's  plane  of  polarization  is  parallel  to  that  of 
the  transmitted  signal.  The  target  return  remains  unchanged.  The  conclusion 
is  that  the  cloud/target  discrimination  technique  will  greatly  improve  the 
capability  of  a  ladar  system  to  distinguish  between  an  atmospheric  cloud  and  a 
target. 
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SECTION  I 


INTRODUCTION 

The  development  of  the  first  laser  in  1960  has  stimulated  multi¬ 
disciplinary  research  in  almost  all  areas  of  science  and  engineering. 

The  list  of  applications  is  seemingly  unlimited,  running  tt>  gamut  from 
retinal  "welding"  and  "bloodless"  surgery  to  power  transmission  systems. 
However,  the  list  will  not  be  expounded  upon  in  this  report.  The  reader  is 
referred  to  Reference  39  for  a  comprehensive  review  of  many  laser  applica¬ 
tions.  While  considerable  attention  has  been  devoted  to  the  application  of 
lasers  in  communication  and  radar  systems  and  to  a  lesser  extent  in 
surveillance  and  reconnaissance  work,  the  effective  use  of  such  systems 
in  the  atmosphere  is  seriously  limited  by  the  severe  attenuation  of  the 
laser  signal  due  to  absorption  by  atmospheric  gases  and/or  scattering 
by  cloud  and  fog  particles.  Furthermore,  in  detection  and  ranging 
applications  (which  broadly  includes  surveillance  and  reconnaissance  appli¬ 
cations,  etc.)  the  cloud  and  fog  particles  may  scatter  the  transmitted 
radiation  in  the  direction  of  the  receiver,  thus  producing  a  false 
return  which  must  be  distinguished  from  the  true  target  return.  It  is  the 
primary  concern  of  this  report  to  present  a  solution  to  the  problem  of 
distinguishing  between  the  target  return  and  the  cloud  return. 

Most  present  day  laser  detection  and  ranging  systems  use  some  form 
of  range  gating  (coincident  pulse  discrimination,  binocular  ranging, 
etc.)  to  eliminate  (not  solve)  the  cloud  return  problem.  While  for 
some  situations  this  may  be  a  satisfactory  compromise  between  the 
advantages  and  disadvantages  of  using  a  laser  system,  the  use  of  range¬ 
gating  techniques  to  eliminate  cloud  returns  is  unacceptable  in  the 
majority  of  potential  applications.  Any  electronic  discrimination 
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technique  based  on  determining  the  difference  in  characteristics  between 
the  cloud  and  target  returns  would  probably  be  very  complicated  and 
extremely  expensive  if  not  impossible,  to  Implement.  What  is  needed, 
then,  is  a  simple,  reliable,  inexpensive  method  for  discriminating 
between  cloud  and  target  returns. 

In  1908  Gustav  Mie  (Ref.  42]  presented  the  results  of  his  research 
on  the  interaction  of  a  linearly  polarized  monochromatic  plane  wave  of 
light  with  a  single  homogeneous  sphere  of  any  description  situated  in 
a  homogeneous  medium.  He  found  that  if  the  diameter  of  the  sphere  is 
of  the  same  order  of  magnitude  as  the  wavelength  of  the  incident  radi¬ 
ation,  then  the  scattered  radiation  is,  in  general,  either  linearly 
or  elliptically  polarized  depending  on  the  plane  of  observation.  At 
laser  wavelengths,  then,  clouds  actually  represent  a  medium  composed 
of  Mie  scattering  particles.  On  the  other  hand,  most  targets — earth  ter¬ 
rain  including  man-made  irregularities,  airborne  objects  such  as  aircraft, 
etc., — can  be  considered  to  be  diffuse  scatterers  at  these  wavelengths. 
Thus,  one  possible  approach  to  developing  a  cloud/target  discrimination 
technique  that  is  both  simple  and  reliable  while  being  relatively 
inexpensive  to  implement  could  be  based  on  the  difference  in  the  pol¬ 
arization  characteristics  of  each  return. 

The  purpose  of  this  report  is  to  describe  an  optical  technique  for 
discriminating  between  a  cloud  return  and  a  target  return,  and  to  present 
the  results  of  a  research  program  to  determine,  experimentally,  the 
validity  of  the  discrimination  technique.  The  feasibility  evaluation 
of  the  discrimination  technique  encompasses  such  an  enormous  amount  of 
multidisciplinary  reference  material  that  we  shall  restrict  our  discussion 
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to  those  subjects  that  are  directly  applicable  to  the  investigation. 
Inasmuch  as  clouds  present  the  most  serious  obstacle  to  the  use  of  a 
ladar  (laser  detection  and  ranging)  system  in  the  atmosphere,  we  begin 
our  report  in  Section  II  with  a  comprehensive  discussion  of  the  processes 
involved  in  the  formation  of  an  atmospheric  cloud.  Section  III  presents 
our  derivation  and  a  discussion  of  the  Mle  scattering  equations.  Section 
IV  presents  a  qualitative  description  and  the  mathematical  derivation 
of  the  discrimination  technique.  Section  V  describes  the  cloud-chamber 
facility  used  in  the  laboratory  verification  phase  of  the  program. 

Section  VI  presents  the  results  of  the  laboratory  verification  and 
field  evaluation  phases  of  the  program.  In  Section  VII  a  discussion  of 
the  feasibility  of  the  discrimination  technique  is  presented. 
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SECTION  II 


THE  PHYSICS  OF  ATMOSPHERIC  AEROSOLS 


2.1  INTRODUCTION 

A  cloud  or  a  fog  can  be  described  as  a  colloidal  suspension  in  air  of 
liquid  water  droplets  or  ice  crystals,  or  a  mixture  of  droplets  and 
crystals,  which  is  frequently  referred  to  as  an  atmospheric  aerosol. 

The  formation  of  an  atmospheric  aerosol  is  dependent  upon  the  conversion 
of  some  of  the  water  vapor  in  the  atmosphere  to  either  the  liquid  or 
solid  state  or  a  combination  thereof.  The  conversion  process  involving 
a  direct  change  from  the  vapor  state  to  the  liquid  state  is  defined  as 
condensation,  while  a  direct  change  from  the  vapor  state  to  the  solid 
state  is  called  sublimation.  As  demonstrated  by  several  researchers 
toward  the  end  of  the  last  century,  the  changes  of  state  of  water  vapor 
in  the  atmosphere  that  result  in  the  formation  of  clouds  are  dependent 
upon,  (a)  the  numbers  and  types  of  condensation  and  sublimation  nuclei 
present,  and  (b)  the  state  of  saturation  of  the  surrounding  air.  The 
degree  of  saturation  necessary  for  the  formation  of  an  aerosol  can  vary 
considerably  depending  on  the  hygroscopic  nature  and  abundance  of  the 
condensation  and  sublimation  nuclei.  The  process  of  condensation  or 
sublimation  of  water  vapor  on  suitable  nuclei  is  called  nucleation. 

For  the  nucleation  process  to  proceed,  it  is  necessary  for  the  air  to  be 
maintained  in  a  saturated  or  nearly  saturated  state  which  in  essence 
requires  some  means  of  continuously  supplying  the  necessary  water  vapor. 

A  state  of  saturation  may  be  achieved  either  by  decreasing  the  tempera¬ 
ture  of  the  air  to  the  point  where  its  capacity  for  water  vapor  is 
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greatly  reduced,  or  by  adding  water  vapor  to  the  air  to  the  extent 
necessary  to  Increase  the  content  to  the  maximum  capacity  possible  at 
the  prevailing  temperature.  A  combination  of  these  processes  may  also 
result  in  a  state  of  saturation.  In  general,  a  cloud  is  formed  when 
air  containing  water  vapor  rises  to  the  higher  altitudes  where  a  con¬ 
tinuous  decrease  in  atmospheric  pressure  results  in  the  expansion  and 
subsequent  cooling  of  the  moist  air  until  some  of  the  water  vapor  con¬ 
denses  or  sublimes  on  the  available  nuclei. 

2.2  THE  CLASSIFICATION  OF  FOGS  AND  CLOUDS 

An  internationally  accepted  classification  system  for  describing 
clouds  is  based  primarily  on  visual  appearance,  and  is  essentially  the 
same  system  as  that  proposed  in  1803  by  a  British  pharmacist,  Luke 
Howard,  [Ref.  25]  and  later  modified  by  Renou  [Ref.  49]  and  Hlldebrandsson 
[Ref.  19]  to  include  considerations  of  altitude.  Since  the  close  of  the 
19th  century,  the  World  Meteorological  Organization  has  been  responsible 
for  any  alterations  to  the  basic  system.  It  publishes  in  the  form  of  an 
atlas  [Ref.  62]  the  definitions  and  photographs  of  the  different  types 
of  clouds.  Ten  genera  of  clouds  are  arranged  into  four  groups  according 
to  their  height  above  the  ground,  as  follows: 

High-Altitude  Clouds 

1.  Cirrus  -  Temperate  Regions:  20,000  to  40,000  feet 

-  Polar  Regions:  Close  to  ground  up  to 
40,000  feet 

2.  Cirrostratus  20,000  to  40,000  feet 

3.  Cirrocuraulue  20,000  to  40,000  feet 
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Middle-Altitude  Clouds  -  6000  feet  up  to  20,000  feet 

4.  Altocumulus 

5.  Altostratus 

Low-Altitude  Clouds  -  close  to  ground  up  to  6,000  feet 

6.  Stratocumulus 

7.  Stratus 

8.  Nimbostratus 

Clouds  with  Vertical  Development 

9.  Cumulus 

a.  Cumulus  humilis  (fair-weather  cumulus) 

Bases:  500  to  10,000  feet 

Tops:  1,000  to  15,000  feet 

l 

b.  Cumulus  congestus 

Bases:  1,000  to  6,000  feet 
Tops:  6,000  to  15,000  feet 

10.  Cumulonimbus 

Bases:  1,000  to  6,000  feet 
Tops:  10,000  to  40,000  feet 

While  an  airflow  crossing  a  hill,  mountain  or  ridge  may  cause, 
orographic  clouds  to  form  that  appear  to  be  different  than  any  of  the 
ten  genera,  the  physical  constitution  of  an  orographic  cloud  is  actually 
similar  to  Altocumulus,  Stratocumulus  or  Cumulus  Clouds  and  they  are 
therefore  classified  under  one  or  other  of  the  ten  genera. 

These  ten  types  of  clouds  may  be  further  subdivided  into  species 
and  varieties  according  to  their  structure  and  transparency.  An  example 
of  two  common  species  of  cumulus  clouds  is  shown  in  the  above  listing. 
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The  following  definitions  of  the  ten  genera  of  clouds,  including  a 
description  of  their  physical  constitution  are  taken  verbatim  from 
Volume  I  of  the  World  Meteorological  Organization's  1956  International 
Cloud  Atlas. 


1.  Cirrus 


(1)  Definition 

Detached  clouds  in  the  form  of  white,  delicate  fila¬ 
ments,  or  white  or  mostly  white  patches  or  narrow  bands. 

These  clouds  have  a  fibrous  (hair-like)  appearance,  or  a  silky 
sheen,  or  both. 

(2)  Physical  Constitution 

Cirrus  is  composed  almost  exclusively  of  ice  crystals. 
These  crystals  are  in  general  very  small,  a  fact  which,  together 
with  their  sparseness,  accounts  for  the  transparency  of  most 
Cirrus  clouds. 

Dense  Cirrus  patches  or  Cirrus  in  tufts  may  never¬ 
theless  contain  ice  crystals  large  enough  to  acquire  an  appreci¬ 
able  terminal  velocity,  so  that  trails  of  considerable  vertical 
extent  may  form.  Sometimes,  though  not  very  frequently,  the 
ice  crystals  in  the  trails  melt  into  small  water  droplets;  the 
trails  are  then  greyish,  in  contrast  with  their  usual  white 
appearance,  and  may  give  rise  to  the  formation  of  a  rainbow. 

The  trails  curve  irregularly  or  slant  as  a  result 
of  wind  shear  and  of  the  variation  in  size  of  the  constituent 
particles;  consequently,  Cirrus  filaments  near  the  horizon  do 
not  appear  parallel  to  it. 

Halo  phenomena  may  occur;  circular  haloes  almost 
never  show  a  complete  ring,  owing  to  the  narrowness  of  the 
Cirrus  clouds. 


2.  Cirrocumulus 


(1)  Definition 

Thin,  white  patch,  sheet  or  layer  of  cloud  without 
shading,  composed  of  very  small  elements  in  the  form  of  grains, 
ripples,  etc.,  merged  or  separate,  and  more  or  less  regularly 
arranged;  most  of  the  elements  have  an  apparent  width  of  less 
than  one  degree. 

(2)  Physical  Constitution 

Cirrocumulus  is  composed  almost  exclusively  of  ice 
crystals;  strongly  supercooled  water  droplets  may  occur  but 
are  usually  rapidly  replaced  by  ice  crystals. 

A  corona  or  irisation  may  sometimes  be  observed. 
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3.  Cirrostratus 


(1)  Definition 

Transparent,  whitish  cloud  veil  of  fibrous  (hair¬ 
like)  or  smooth  appearance,  totally  or  partly  covering  the 
sky,  and  generally  producing  halo  phenomena. 

(2)  Physical  Constitution. 

Cirrostratus  is  composed  mainly  of  ice  crystals. 

The  smallness  of  these  crystals,  their  sparseness  and  the  fact 
that  Cirrostratus  has  at  most  only  a  moderate  depth,  account 
for  the  transparency  of  this  cloud  through  which  the  outline 
of  the  sun  is  visible,  at  least  when  the  latter  is  not  too 
close  to  the  horizon. 

In  certain  types  of  Cirrostratus,  some  of  the  ice 
crystals  are  large  enough  to  acquire  an  appreciable  terminal 
velocity,  so  that  trailing  filaments  are  formed,  which  give 
these  Cirrostratus  clouds  a  fibrous  appearance. 

Halo  phenomena  are  often  observed  in  thin  Cirro- 
stratus;  sometimes  the  veil  of  Cirrostratus  is  so  thin  that  a 
halo  provides  the  only  indication  of  its  presence. 

4.  Altocumulus 


(1)  Definition 

White  or  gr^y,  or  both  white  and  grey,  patch,  sheet 
or  layer  of  cloud,  generally  with  shading,  composed  of  laminae, 
rounded  masses,  volls,  etc.,  which  are  sometimes  partly  fibrous 
or  diffuse  and  which  may  or  may  not  be  merged;  most  of  the 
regularly  arranged  small  elements  usually  have  an  apparent 
width  between  one  and  five  degrees. 

(2)  Physical  Constitution 

Altocumulus  is,  at  least  in  the  main,  almost  in¬ 
variably  composed  of  water  droplets.  This  is  evident  from  the 
fairly  low  transparency  of  the  microscopic  elements,  and  from 
the  fact  that  the  latter  show  sharp  outlines  when  separate. 
Nevertheless,  when  the  temperature  is  very  low,  ice  crystals 
may  form.  If  the  droplets  then  evaporate,  the  cloud  becomes 
entirely  an  ice  cloud  and  its  macroscopic  elements  cease  to 
present  sharp  outlines.  The  formation  of  ice  crystals  may 
take  place  in  all  species  of  Altocumulus;  it  occurs  most 
frequently  in  Altocumulus  castellanus  and  flaccus. 

A  corona  or  irisation  is  often  observed  in  thin 
parts  of  Altocumulus.  Parhelia  or  luminous  pillars  are  some¬ 
times  seen  in  Altocumulus,  indicating  the  presence  of  tubular¬ 
shaped  ice  crystals. 
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5.  Altostratus 


(1)  Definition 

Greyish  or  bluish  cloud  sheet  or  layer  of  striated, 
fibrous  or  uniform  appearance,  totally  or  partly  covering  the 
sky,  and  having  parts  thin  enough  to  reveal  the  sun  at  least 
vaguely,  as  through  ground  glass.  Altostratus  does  not  show 
halo  phenomena. 

(2)  Physical  Constitution 

Altostratus  nearly  always  appears  as  a  layer  of 
great  horizontal  extent  [several  tens  or  hundreds  of  kilo¬ 
meters  (several  tens  or  hundreds  of  miles)]  and  fairly  con¬ 
siderable  vertical  extent  [several  hundreds  or  thousands  of 
meters  (several  hundreds  or  thousands  of  feet)].  It  is  com¬ 
posed  of  water  droplets  and  ice  crystals.  In  the  most  com¬ 
plete  case,  three  superposed  parts  may  be  distinguished, 
namely: 

a.  An  upper  part,  composed  wholly  or  mainly  of 
ice  crystals, 

b.  A  middle  part,  composed  of  a  mixture  of  ice 
crystals,  snow  crystals  or  snowflakes  and  supercooled  water 
droplets, 

c.  A  lower  part,  composed  wholly  or  mainly  of 
ordinary  or  supercooled  water  droplets  or  drops. 

In  some  cases,  the  cloud  may  consist  of  only  two 
parts,  either: 

-  an  upper  part  like  (a)  and  a  lower  part  like  (c)  or 

-  an  upper  part  like  (b)  and  a  lower  part  like  (c). 

Less  frequently,  the  entire  cloud  may  also  be  like 

(a)  or  like  (b)  alone. 

The  constituent  particles  in  the  lower  part  of 
Altostratus  are  so  numerous  that  the  outline  of  the  sun  or 
moon  is  always  dimmed  and  the  surface  observer  never  sees 
halo  phenomena.  In  the  thickest  parts,  the  position  of  the 
luminary  may  be  completely  concealed. 

Raindrops  or  snowflakes  are  often  present  in  Alto- 
stratus  and  below  its  base.  When  precipitation  reaches  the 
ground,  it  is  generally  of  the  "continuous"  type  and  in  the 
form  of  rain,  snow  or  ice  pellets. 

6.  Nimbostratus 


(1)  Definition 

Grey  cloud  layer,  often  dark,  the  appearance  of 
which  is  rendered  diffuse  by  more  or  less  continuously  falling 
rain  or  snow,  which  in  most  cases  reaches  the  ground.  It  is 
thick  enough  throughout  to  blot  out  the  sun. 
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Low,  ragged  clouds  frequently  occur  below  the  layer, 
with  which  they  may  or  may  not  merge. 

(2)  Physical  Constitution 

Nimbo stratus  generally  covers  a  wide  area  and  is  of 
great  vertical  extent.  It  is  composed  of  water  droplets 
(sometimes  supercooled)  and  raindrops,  of  snow  crystals  and 
snowflakes,  or  of  a  mixuure  of  these  liquid  and  solid  particles. 
The  high  concentration  of  particles  and  the  great  vertical 
extent  of  the  cloud  prevent  direct  sunlight  from  being  ob¬ 
served  through  it.  The  cloud  produces  rain,  snow  or  ice  pel¬ 
lets  which,  however,  do  not  necessarily  reach  the  ground. 

7.  Stratocumulus 


(1)  Definition 

Grey  or  whitish,  or  both  grey  and  whitish,  patch, 
sheet  or  layer  of  cloud  which  almost  always  has  dark  parts, 
composed  cf  tessellations,  rounded  masses,  rolls,  etc.,  which 
are  non- fibrous  (except  for  virga)  and  which  may  or  may  not 
be  merged;  most  of  the  regularly  arranged  small  elements  have 
an  apparent  width  of  more  than  five  degrees. 

(2)  Physical  Constitution 

Stratocumulus  is  composed  of  water  droplets,  some¬ 
times  accompanied  by  raindrops  or  snow  pellets  and,  more 
rarely,  by  snow  crystals  and  snowflakes.  Any  ice  crystals 
present  are  usually  too  sparse  to  give  the  cloud  a  fibrous 
appearance;  during  extremely  cold  weather,  however,  Strato¬ 
cumulus  may  produce  abundant  ice  crystal  virga  which  may  be 
accompanied  by  a  halo.  When  Stratocumulus  is  not  very  thick, 
a  corona  or  irisation  is  sometimes  observed. 

8.  Stratus 


(1)  Definition 

Generally,  grey  cloud  layer  with  a  fairly  uniform 
base,  which  may  give  drizzle,  ice  prisms  or  snow  grains. 

When  the  sun  is  visible  through  the  cloud,  its  outline  is 
clearly  discernible.  Stratus  does  not  produce  halo  pheno¬ 
mena  except,  possibly,  at  very  low  temperatures. 

Sometimes  Stratus  appears  in  the  form  of  ragged 

patches. 

(2)  Physical  Constitution 

Stratus  is  usually  composed  of  small  water  droplets; 
this  cloud  may,  when  very  thin,  produce  a  corona  round  the  sun 
or  moon.  At  very  low  temperatures,  Stratus  may  consist  of 
small  ice  particles.  The  ice  cloud  is  usually  thin  and  may, 
on  rare  occasions,  produce  halo  phenomena. 
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Stratus,  when  dense  or  thick,  Oxten  contains 
drizzle  droplets  and  sometimes  ice  prisms  or  snow  grains,  it 
may  then  have  a  dark  or  even  a  threatening  appearance.  Stratus 
with  a  low  optical  thickness,  when  observed  at  more  than  90 
degrees  from  the  sun,  often  shows  a  more  or  less  smoky, 
greyish  tint  like  that  of  fog. 

9.  Cumulus 


(1)  Definition 

Detached  clouds,  generally  dense  and  with  sharp 
outlines,  developing  vertically  in  the  form  of  rising  mounds, 
domes  or  towers,  or  which  the  bulging  upper  part  often  re¬ 
sembles  a  cauliflower.  The  sunlit  parts  of  these  clouds  are 
mostly  brilliant  white;  their  base  is  relatively  dark  and 
nearly  horizontal. 

Sometimes  Cumulus  is  ragged. 

(2)  Physical  Constitution 

Cumulus  is  composed  mainly  of  water  droplets.  When 
of  great  vertical  extent,  Cumulus  may  release  precipitation 
in  the  form  of  rain  showers. 

Ice  crystals  may  form  in  those  parts  of  a  cumulus 
in  which  the  temperature  is  well  below  0°C;  they  grow  at  the 
expense  of  evaporating  supercooled  water  droplets,  thereby 
transforming  the  cloud  into  Cumulonimbus.  In  cold  weather, 
when  che  temperature  in  the  entire  cloud  is  well  below  0°C 
(32°F),  this  process  leads  to  the  degeneration  of  the  cloud 
into  diffuse  trails  of  snow. 

10.  Cumulonimbus 


(1)  Definition 

Heavy  and  dense  cloud,  with  a  considerable  vertical 
extent,  in  the  form  of  a  mountain  or  high  towers.  At  least 
part  of  its  upper  portion  is  usually  smooth,  or  fibrous  or 
striated,  and  nearly  always  flattened;  this  part  often  spreads 
out  in  the  shape  of  an  anvil  cr  vast  plume. 

Under  the  base  of  this  cloud  which  often  is  very 
dark,  there  are  frequently  low  ragged  clouds  either  merged 
with  it  or  not,  and  precipitation  sometimes  in  the  form  of 
virga. 

(2)  Physical  Constitution 

Cumulonimbus  is  composed  of  water  droplets  and, 
especially  in  its  upper  portion,  of  ice  crystals.  It  also 
contains  large  raindrops  and,  often,  snowflakes,  snow  pellets, 
ice  pellets  or  hailstones.  The  water  droplets  and  raindrops 
may  be  substantially  supercooled. 
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All.  other  atmospheric  phenomena  have  been  grouped  under  the  heading 
of  meteors.  The  World  Meteorological  Organization  defines  a  meteor  as 
a  phenomenon,  other  than  a  cloud,  observed  in  the  atmosphere  or  on  the 
surface  of  the  earth,  which  consists  of.  a  precipitation,  a  suspension 
or  a  deposit  of  aqueous  or  non-aqueous  liquid  or  solid  particles,  or 
a  phenomenon  of  the  nature  of  an  optical  or  electrical  manifestation. 
Meteors  are  divided  into  four  groups,  hydrometeors,  lithometeors, 
photometeors,  and  electrometeors,  according  to  the  nature  of  their  con¬ 
stituent  particles  or  the  physical  processes  involved  in  their  occurrence. 
The  following  definitions  are  taken  verbatim  from  the  International 
Cloud  Atlas  of  1956: 

(1)  Hydrometeors 

A  hydrometeor  Is  a  meteor  consisting  of  an  ensemble  of 
liquid  or  solid  water  particles,  falling  through  or  suspended 
in  the  atmosphere,  blown  by  the  wind  from  the  earth's  surface, 
or  deposited  on  objects  on  the  ground  or  in  the  free  air. 

Hydrometeors  consisting  of  ensembles  of  falling  particles 
(rain,  drizzle,  snow,  snow  pellets,  snow  grains,  ice  pellets, 
ice  prisms),  originate  mostly  in  clouds;  the  particles  may 
reach  the  earth's  surface  (precipitation)  or  they  may  completely 
evaporate  during  their  fall  (virga) 

(a)  Rain.  Precipitation  of  liquid  water  particles, 
either  in  the  form  of  drops  of  more  than  0.5  mm  (0.02  in.) 
diameter  or  of  smaller  widely  scattered  drops.  Freezing  rain 
is  defined  as  rain  the  drops  of  which  freeze  on  impact  with 
the  ground  or  with  objects  on  the  earth's  surface  or  with 
aircraft  in  flight. 

(b)  Drizzle.  Fairly  uniform  precipitation  composed 
exclusively  of  fine  drops  of  water  [diameter  less  than  0.5  mm 
(0.02  in.)],  very  close  to  one  another.  Freezing  drizzle  is 
defined  as  drizzle,  the  drops  of  which  freeze  on  impact  with 
the  ground  or  with  objects  on  the  earth's  surface  or  with  air¬ 
craft  in  flight. 

(c) .  Snow.  Precipitation  of  ice  crystals,  most  of  which 
are  branched  (sometimes  star-shaped). 

(d)  Snow  Pellets.  Precipitation  of  white  and  opaque 
grains  of  ice.  These  grains  are  spherical  or  sometimes  conical; 
their  diameter  is  about  2-5  mm  (0. 1-0.2  in.). 
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(e)  Snow  Grains.  Precipitation  of  very  small  white  and 
opaque  grains  of  ice.  These  grains  are  fairly  flat  or  elon¬ 
gated;  their  diameter  is  generally  less  than  1  ran  (0.04  in.). 

(f)  Ice  Pellets.  Precipitation  of  transparent  or  trans¬ 
lucent  pellets  of  ice,  which  are  spherical  or  irregular, 
rarely  conical,  and  which  have  a  diameter  of  5  raa  (0.2  in.) 

or  less. 

(g)  Hail.  Precipitation  of  small  balls  or  pieces  of 
ice  (hailstones)  with  a  diameter  ranging  from  5  to  50  mm 
(0.2  to  2.0  in.)  or  sometimes  more,  falling  either  separately 
or  agglomerated  into  irregular  lumps. 

(h)  Ice  Prisms.  A  fall  of  unbranched  ice  crystals,  in 
the  form  of  needles,  columns,  or  plates,  often  so  tiny  that 
they  seem  to  be  suspended  in  the  air.  These  crystals  may 
fall  from  a  cloud  or  from  a  cloudless  sky. 

Hydrometeors  consisting  of  ensembles  of  particles  sus¬ 
pended  in  the  air  (fog,  mist)  are  similar  to  clouds.  They 
are,  however,  considered  as  meteors,  owing  to  their  occurrence 
at  or  near  the  earth's  surface. 

(i)  Fog.  A  suspension  of  very  small  water  droplets 
in  the  air,  generally  reducing  the  horizontal  visibility  at 
the  earth's  surface  to  less  than  1  km  (5/8  mile). 

(j)  Ice  Fog.  A  suspension  of  numerous  ice  crystals  in 
the  air,  reducing  the  visibility  at  the  earth's  surface. 

(k)  Mist.  A  suspension  in  the  air  of  microscopic  water 
droplets  or  wet  hygroscopic  particles,  reducing  the  visibility 
at  the  earth ' s  surface . 

The  hydrometeors  consisting  of  ensembles  of  particles 
blown  by  the  wind  from  the  earth's  surface  (drifting  snow; 
blowing  snow  and  spray)  are  generally  confined  to  the  lowest 
layers  of  the  atmosphere. 

(l)  Drifting  Snow.  An  ensemble  of  snow  particles 
raised  by  the  wind  to  small  heights  above  the  ground.  The 
visibility  is  not  sensibly  diminished  at  eye  level.  Eye 
level  is  defined  as  1.8  m  (6  feet)  above  the  ground. 

(m)  Blowing  Snow.  An  ensemble  of  snow  particles  raised 
by  the  wind  to  moderate  or  great  heights  above  the  ground. 

The  horizontal  visibility  at  eye  level  is  generally  very  poor. 

(n)  Spray.  An  ensemble  of  water  droplets  torn  by  the 
wind  from  the  surface  of  an  extensive  body  of  water,  generally 
from  the  crests  of  waves,  and  carried  up  a  short  distance 
into  the  air. 

The  hydrometeors  in  the  form  of  deposits  (dew,  hoar-frost 
rime  and  glaze)  occur  either  in  the  form  of  aggregate'  of 
particles,  more  or  less  individually  discernible  in  spite  of 
the  fact  that  they  are  often  fused  together  (hoar-frost, 
rime),  or  as  homogeneous  and  smooth  layers  in  which  no 
particle  structure  can  be  discerned  (glaze). 

(o)  Dew.  A  deposit  of  water  drops  on  objects  at  or 
near  the  ground,  produced  by  the  condensation  of  water  vapor 
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from  the  surrounding  clear  air.  White  dew  is  a  deposit  of 
white  frozen  dew  drops. 

(p)  Hoar-Frost.  A  deposit  of  ice  having  a  crystalline 
appearance,  generally  assuming  the  form  of  scales,  needles, 
feathers  or  fans. 

(q)  Rime.  A  deposit  of  ice,  composed  of  grains  more  or 
less  separated  by  trapped  air,  sometimes  adorned  with  crystal¬ 
line  branches. 

(r)  Rime  (clear  ice).  A  generally  homogeneous  and  trans¬ 
parent  deposit  of  ice  formed  by  the  freezing  of  supercooled 
drizzle  droplets  or  raindrops  on  objects  the  surface  tempera¬ 
ture  of  which  is  below  or  slightly  above  0eC  (32°F). 

(3)  Spout.  A  phenomenon  consisting  of  an  often  violent 
whirlwind,  revealed  by  the  presence  of  a  cloud  column  or 
inverted  cloud  core  (funnel  cloud),  protruding  from  the  base 
of  a  Cumulonimbus,  and  of  a  "bush"  composed  of  water  droplets 
raised  from  the  surface  of  the  sea  or  of  dust,  sand  or  litter, 
raised  from  the  ground. 

(2)  Lithometeors 

A  lithometeor  is  a  meteor  consisting  of  an  ensemble  of 
particles  most  of  which  are  solid  and  non-aqueous.  The 
particles  are  more  or  less  suspended  in  the  air,  or  lifted 
by  the  wind  from  the  ground. 

The  lithometeors  which  have  more  or  less  the  character 
of  suspensions  in  the  atmosphere  are  haze,  dust  haze  and 
smoke;  they  consist  of  very  small  dust  particles,  of  sea-salt 
particles  or  of  combustion  products  (e.g.,  from  forest  fires). 

(a)  Haze.  A  suspension  in  the  air  of  extremely  small, 
dry  particles  invisible  to  the  naked  eye  and  sufficiently 
numerous  to  give  the  air  an  opalescent  appearance. 

(b)  Dust  Haze.  A  suspension  in  the  air  of  dust  or 
small  sand  particles,  raised  from  the  ground  prior  to  the 
t  me  of  observation  by  a  dust  storm  or  sandstorm. 

(c)  Smoke.  A  suspension  in  the  air  of  small  particles 
produced  by  combustion. 

The  lithometeors  resulting  from  the  action  of  the  wind 
are  drifting  and  blowing  dust  or  sand,  dust  storm  or  sand¬ 
storm,  and  dust  whirl  or  sand  whirl. 

(d)  Drifting  dust  or  drifting  sand.  Dust  or  sand,  raised 
by  the  wind  to  small  heights  above  the  ground.  The  visibility 
is  not  sensibly  diminished  at  eye  level. 

(e)  Blowing  Dust  or  Blowing  Sand.  Dust  or  sand,  raised 
by  the  wind  to  moderate  heights  above  the  ground.  The  hori¬ 
zontal  visibility  at  eye  level  is  sensibly  reduced. 

(f)  Dust  Storm  or  Sandstorm.  An  ensemble  of  particles 
of  dust  or  sand  energetically  lifted  to  great  heights  by  a 
strong  and  turbulent  wind. 

(g)  Dust  Whirl  or  Sand  Whirl.  An  ensemble  of  particles 
of  dust  or  sand,  sometimes  accompanied  by  small  litter,  raised 
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from  the  ground  in  the  form  of  a  whirling  column  of  varying 
height  with  a  small  diameter  and  an  approximately  vertical  axis. 

(3)  Photometeors. 

A  photometeor  is  a  luminous  phenomenon  produced  by  the 
reflection,  refraction,  diffraction  or  interference  of  light 
from  the  sun  or  the  moon. 

Photometeors  are  observed  in  more  or  less  clear  air 
(mirage,  shimmer,  scintillation,  green  flash,  twilight  colors), 
or  on  inside  clouds  (halo  phenomena,  coronae.  lrisation, 
glory)  and  on  or  inside  certain  hydrometeors  or  lithometeors 
(^lory,  rainbow,  fog  bow,  Bishop's  ring,  crepuscular  rays). 

(4)  Electrometeor 8 

An  electrometeor  is  a  visible  or  audible  manifestation 
of  atmospheric  electricity. 

Electrometeors  either  correspond  to  discontinuous  elec¬ 
trical  discharges  (lightning,  thunder)  or  occur  as  more  or 
less  continuous  phenomena  (Saint  Elmo's  fire,  polar  aurora). 

2.3  THE  THERMODYNAMICS  OF  THE  FORMATION  OF  ATMOSPHERIC  AEROSOLS 

The  processes  involved  in  the  formation  of  a  cloud  are  governed  by 
the  laws  of  thermodynamics  and  are  functions  of  changes  in  the  pressure, 
temperature  and  volume  relationships  of  the  atmosphere.  While  the  com¬ 
position  of  the  atmosphere  consists  of  a  mixture  of  gases  and  solid  and/ 
or  liquid  particles,  it  is  the  minor  constituent  water  vapor  that  plays 
the  dominant  role  on  the  formation  of  clouds.  The  fundamental  laws 
governing  the  behavior  of  real  gases  are  based  on  the  kinetic  theory  of 
an  ideal  gas  and  the  empirical  relationships  derived  by  the  early  inves¬ 
tigators  of  gaseous  behavior. 

One  of  the  first  fundamental  experiments  in  gaseous  behavior  was 
performed  by  Robert  Boyle  in  1660.  He  found  that  if  the  temperature 
and  mass  of  a  gas  were  held  constant  while  measuring  its  volume  at 
several  different  pressures,  the  product  of  the  pressure  and  volume 
remained  approximately  constant.  That  is,  for  a  fixed  mass  of  a  real 
gas  kept  at  a  constant  temperature 
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pV  =  constant 


(1) 


where  p  and  V  represent  the  pressure  and  volume.  Since  the  behavior  of 
real  gases  follows  this  pressure-volume  relationship  only  approximately, 
it  was  mathematically  convenient  to  define  an  ideal  gas  where  for  a 
fixed  mass  at  a  given  temperature 

pV  -  constant  (2) 

This  is  referred  to  as  Boyle's  Law. 

In  1802  Gay-Lussac  presented  his  findings  relating  the  changes  in 
the  volume  of  a  gas  to  changes  in  the  temperature  while  the  pressure  is 
held  constant.  Mathematically,  if  VQ  represents  the  volume  at  a  reference 
temperature,  tQ, 

V  *  Vq  [l+k(t-to)]  (3) 

where  V  is  the  volume  at  temperature  t  and  k  is  defined  as  the  volume 
coefficient  of  expansion.  It  is  found  that  if  k  was  measured  at  0°C 
for  a  number  of  gases  at  different  pressures,  its  value  remained  very 
nearly  constant  for  all  the  gases  and  approximately  equal  to  1/273  per 
centigrade  degree.  Thus,  at  a  reference  temperature  of  0°C,  Equation 
(3)  becomes 

V  «  Vo(l+kot)  (4) 

where  kQ  -  1/273  per  centigrade  degree. 

As  in  Boyle's  relationship,  it  was  mathematically  convenient  to 
extend  the  definition  of  an  ideal  gas  where  kQ  *  1/273  per  centigrade 
degree  and 

V  -  Vo(l+kot)  <5) 

This  relation  is  referred  to  as  Gay-Lussac's  Law. 
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We  can  combine  the  laws  of  Boyle  and  Gay-Lussac  to  obtain  an 
equation  of  state  for  an  ideal  gas  expressed  as  a  relationship  among 
pressure,  volume,  and  temperature.  Consider  a  given  mass  of  an  ideal 
gas.  If  the  temperature  is  varied  from  an  initial  value  tQ  ■  0°C 
to' a  value  t  at  constant  pressure  p^,  then  according  to  Gay-Lussac's 
law,  the  volume  at  temperature  t,  V  can  be  expressed  as 

Vt  -  VQ(l+kot)  (6) 

Now,  according  to  Boyle's  law,  if  we  vary  the  pressure  to  a  value  p 
while  holding  the  temperature  constant  at  t,  the  new  volume  V  can  be 
expressed  as 

PV  -  PQVt  (7) 

Substituting  the  expression  for  V  from  Equation  (6)  gives 

Pv  -  PQV0(l+kot)  (8) 

or 

py  -  PQV0ko(t+l/ko)  (9) 

If  we  assume  t  -  0°C  and  p  >  1  atm,  then  it  can  be  easily  shown 
o  o 

that 

p  V  k  -  nR  (10) 

o  o  o 

where  n  is  the  number  of  gram  molecules  contained  in  the  volume  Vq  and 
R  is  defined  as  the  universal  gas  constant. 

Now,  since  l/kQ  -  2738C,  the  term  in  parentheses  represents  a 
temperature  of  t  +  273  degrees  centigrade  which  defines  the  Kelvin 
temperature  scale.  With  T  ■  t  +  273  where  T  is  expressed  in  degrees 
Kelvin,  Equation  (9)  can  now  be  written  as 

pV  -  nRJ  (11) 
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and  ts  known  as  the  equation  of  state  of  an  ideaL  gas. 

In  deriving  the  equation  of  atate  for  atmospheric  air,  it  is 
preferable  to  express  Equation  (11)  in  other  forms.  Let  m  represent 
the  mass  and  M  the  gram-molecular  weight  of  the  volume  V,  then  by 
definition  the  density  is  given  by 


„  m 

p  ■  — 

M  V 

(12) 

m 

n  M 

(13) 

and 

the  specific  volume 

a  can  be  written  as 

1  V 

a  *  —  =  — 

p  m 

(14) 

The 

molecular  volume  v 

is  defined  as 

V 

v  =  — 
n 

(15) 

Therefore,  substituting 

these  relationships  into  Equation 

(11)  gives 

PV  -  ~  RT 

(16) 

or 

p  ■  °  fT 

(17) 

or 

p“  ■  ST 

(18) 

or 

pv  ■  RT 

(19) 

It  is  Impossible  to  measure  or  calculate  the  molecular  weight  of 
atmospheric  air  because  in  reality,  an  air  molecule  does  not  exist  since 
air  consists  of  a  mixture  of  gases;  however,  for  purposes  of  this  study, 
air  can  be  considered  to  be  composed  of  the  following  two  gases:  (1) 
water  vapor  and  (2)  pure  dry  air  (includes  all  atmospheric  gases  except 
water  vapor).  Now,  we  can  assume  that  the  total  pressure  exerted  by 
the  atmosphere  is  governed  by  Dalton's  Law  which  states  that  the  total 
pressure  exerted  by  a  mixture  of  gases  is  equal  to  the  sum  of  the 
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pressures  exerted  by  each  constituent  of  the  mixture.  The  individual 
pressures  exerted  by  each  constituent  are  called  partial  pressures. 
Dalton  discovered  that  the  partial  pressure  of  each  constituent  is  the 
same  as  would  be  the  actual  pressure  of  that  constituent  if  it  alone 
occupied  the  entire  volume  at  the  same  temperature  as  the  mixture. 

We  can  now  write  the  equation  of  state  for  water  vapor  as 

P  ■  p  “  T  (20) 

w 

Where  the  subscript  w  refers  to  water  vapor  and  pw  is  the  partial  pres¬ 
sure  of  the  water  vapor,  p^  is  the  vapor  density  or  absolute  humidity 

(mass  of  water  vapor  per  unit  volume)  and  M  is  the  molecular  weight  of 

w 

water. 

Similarly,  for  dry  air,  the  equation  of  state  can  be  written  as 


Where  the  subscript  d  refers  to  dry  air  and  pd  is  the  sum  of  partial 
pressures  of  the  constituent  gases,  is  the  sum  of  the  individual 
densities,  and  is  the  mean  molecular  weight  of  the  mixture. 

Now,  according  to  Dalton's  Law,  the  total  atmospheric  pressure  p 
is  given  by 

P  *  P„  +  Pd  (22) 

When  the  expressions  for  p^  and  p^  from  Equations  (20)  and  (21)  are 
substituted  in  Equation  (22),  we  obtain 


which  can  be  written 


„  R  rp 

Pw  M  T 
w 


pd  M  ,T 
d 


(23) 
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»-5»d  '1  +  &  <5T>1  1  W> 

d  Hd  w 

Now,  Substituting  the  relation  P  "  P„  +  Pd  into  Equation  (24),  gives 


which  can  be  written 

p-‘,M:(1-r>  [l  +  <ir>  <r>i  T  <26> 

d  d  w 

and  is  one  representation  of  the  equation  of  state  for  atmospheric  air. 
While  exact  measurements  indicate  that  the  equation  of  state  for  atmos¬ 
pheric  air  is  only  approximately  valid,  it  is  sufficiently  accurate  for 
most  meteorological  calculations.  And,  therefore,  the  real  gases  of 
the  atmosphere  may  be  treated  as  ideal  gases  for  most  computational 
purposes. 

We  know  from  basic  thermodynamics  that  a  substance  (system)  may  be 
caused  to  change  state  by  the  performance  of  work  or  the  transfer  of 
heat.  The  exact  process  involved  in  the  change  is  dependent  on  the 
definition  of  the  system  and  its  surroundings.  The  expansion  or 
compression  of  a  gas  is  accomplished  by  the  performance  of  work  and 
whenever  a  substance  such  as  1^0  changes  phase  (melts,  freezes,  evapor¬ 
ates,  condenses  or  sublimes),  the  change  of  state  is  brought  about  by 
a  transfer  of  heat.  The  kinetic  theory  of  gases  states  that  the 
internal  (kinetic)  energy  of  a  gas  is  directly  proportional  to  the  mean 
velocity  of  the  molecules  and  the  velocity  of  these  molecules  can  be 
increased  (decreased)  by  either  heating  (cooling)  or  by  compressing 
(expanding)  the  gas.  In  addition,  the  first  law  of  thermodynamics  states 
that  a  change  in  the  internal  energy  of  a  gas  can  be  accomplished  by  the 
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transfer  of  heat  or  by  performing  work  on  the  gas,  or  by  a  combination 
of  the  two  processes.  We  can  therefore  use  the  first  law  of  thermo¬ 
dynamics  to  describe  the  processes  involved  in  the  formation  of  a  cloud. 
If  the  formation  of  the  cloud  occurs  sufficiently  far  above  the  ground 
so  that  there  is  no  solid  or  liquid  surface  to  add  or  remove  heat,  we  can 
assume  that  the  formation  of  the  cloud  is  essentially  an  adiabatic 
process  which  is  determined  primarily  by  the  pressure  changes  experienced 
by  a  parcel  of  air  in  vertical  motion. 

If  we  represent  the  change  in  heat  content  by  dQ  and  the  change  in 
internal  energy  by  dU,  both  positive  for  an  increase  and  negative  for  a 
decrease,  and  4dW  represents  the  work  done  on  the  gas,  then  the  first 
law  of  thermodynamics  can  be  written  as 

dQ  -  dU  +  dW  (27) 

It  can  easily  be  shown  that  for  a  system  where  work  is  performed 
by  means  of  expansion  or  compression  only  (as  in  the  atmosphere),  then 

dW  =  pdV  (28) 

and 

dQ  =  dU  +  pdV  (29) 

In  the  atmosphere  we  are  dealing  with  fixed  masses  of  air,  and  if 
we  divide  Equation  (29)  by  the  mass  of  the  system  we  obtain  another  form 
of  the  first  law: 

dq  ■  du  +  pda  (30) 

where  dq  is  the  change  in  heat  content  per  unit  mass,  du  is  the  specific 
internal  energy,  and  da  is  the  change  in  specific  volume. 

Now,  we  know  by  definition  that  an  ideal  gas  satisfies  the  perfect 
ga3  law.  It  can  also  be  shown  by  the  kinetic-molecular  theory  of  gases 
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chat  the  internal  energy  of  an  ideal  gas  is  a  function  of  temperature 
only.  While  no  real  gas  obeys  these  two  conditions  exactly,  the  degree 
“O  which  these  conditions  are  satisfied  has  been  detormined  experimentally 
for  most  real  gases.  The  results  of  these  experiments  Indicate  that 
for  most  meteorological  computations  the  gases  of  the  atmosphere  can 
be  assumed  to  be  ideal  and  that  their  internal  energy  is  a  function  of 
temperature  only. 

The  temperature  of  a  gas  can  be  changed  by  varying  the  pressure 
while  the  volume  is  held  constant,  by  varying  the  volume  while  the 
pressure  is  held  constant,  or  by  simultaneously  varying  both  the 
pressure  and  volume.  In  each  situation  the  amount  of  heat  necessary  for 
the  change  is  different;  however,  only  the  first  two  cases  are  of 
practical  interest  and  by  convention  we  denote  the  specific  heat 
capacity  of  a  gas  at  constant  volume  as  cv  and  the  specific  heat  capacity 
at  constant  pressure  as  cp.  For  a  constant  volume  process,  i.e.,  dV  »  0, 
it  has  been  shown  experimentally  that 

dq  *»  du  *  cydT  (31) 

The  value  of  cv  determined  experimentally  has  been  found  to  be  constant 
and  Equation  (3i)  is  accepted  as  the  definition  of  the  change  in  internal 
energy  and  may  be  used  whether  the  volume  is  held  constant  or  not.  Now, 
if  the  experiment  is  repeated  while  holding  the  pressure  constant  (i.e. 
the  same  amount  of  the  same  gas  used  in  the  derivation  of  Equation  (31) 
is  heated  to  the  same  temperature  by  varying  the  volume  while  the  pressure 
is  held  constant) ,  we  obtain  the  following  expression 

c  dT  *  c  dT  +  pda  (32) 

p  v 

To  obtain  a  relationship  between  the  two  specific  heats,  Equation 
(32)  may  be  written  in  the  form 
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(cp-Cv)dT  -  pda 


(33) 


and  by  substituting  the  equation  of  state  Equation  (18)  in  its  differen¬ 
tial  form,  i.e.,  by  substituting 

pda  +  ado  -  ^  dT  (34) 


We  are  now  in  a  position  to  express  the  first  law  of  thermodynamics 
in  terms  of  the  easily  measured  atmospheric  variables  of  temperature 
and  pressure.  Starting  with  the  first  law  in  its  general  form 

dq  •»  c^dT  +  pda  (36) 

we  substitue  the  value  for  a  from  Equation  (18),  for  pda  from  Equation 
(34)  and  the  relationship  of  Equation  (35)  to  obtain 

d,  -  cpdT  -  |  T  f  (37) 


Equation  (37)  represents  the  first  law  of  thermodynamics  for  air 
in  terms  of  temperature  and  pressure  only. 

Since  an  adiabatic  process,  i.e.,  one  in  which  dq  =  0,  is  the 
dominant  mechanism  in  the  formation  of  atmospheric  clouds  Equation  (37) 
can  be  simplified  accordingly  to  obtain  several  important  relationships 
that  describe  the  formation  process.  A  complete  mathematical  derivation 
is  quite  lengthy  and  will  not  be  presented  here.  The  reader  is  referred 
to  basic  and  advanced  textbooks  on  meteorology,  such  as  Hess  [Ref.  21] 
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or  Byers  [Ref.  10]  fora  comprehensive  treatment  of  the  formation  of  an 
atmospheric  cloud.  Briefly,  the  results  of  a  comprehensive  mathematical 
treatment  show  that  a  cloud  is  formed  by  the  lifting  of  a  parcel  of 
moist  air  which  is  cooled  by  adiabatic  expansion  at  the  lower  pressures 
to  a  temperature  (dev  point)  where  some  of  the  water  vapor  begins  to 
condense  into  a  fog  «r  cloud  of  water  droplets. 

Looking  at  Equation  (26),  we  can  now  define  several  of  the  more 
important  meteorological  variables  that  are  used  to  describe  the  thermo¬ 
dynamic  conditions  n  a  given  region  of  the  atmosphere.  The  following 
definitions  are  derived  from  measurements  taken  in  air  over  a  plane 
surface  of  pure  watsc  where  the  water  and  air  space  above  it  are  at  the 
same  temperature. 

1.  Vapor  oresstre,  p^:  The  vapor  pressure  is  simply  the  partial 
pressure  exerted  by  the  water  vapor  in  the  air  space.  When  the  vapor 
pressure  equals  the  yapor  tension  of  the  water  surface,  the  air  is  said 
to  be  saturated  witi  water  va.'or.  If  the  vapor  pressure  is  less  than 
the  vapor  tension,  the  air  is  unsaturated.  For  the  state  of  saturation 
the  vapor  pressure  is  usually  referred  to  as  the  saturation  vapor 
pressure  and  it  has  been  determined  empirically  that  it  is  a  function 
of  temperature  independent  of  pressure. 

2.  Absolute  humidity,  pw:  The  absolute  humidity  is  defined  as 
the  water  vapor  density,  pw* 

3.  Mixing  ratio,  io:  The  nixing  ratio  is  defined  as  the  ratio  of 
the  mass  of  water  vapor  present  in  the  total  air  mixture  to  the  mass  of 
the  dry  air  alone.  Thus 


m 
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Since,  from  Dalton's  Law,  all  the  gases  occupy  the  same  volume  and 


p  *  m/V,  Equation  (38)  can  be  expressed  as 

P. 


0)  ■ 


w 


(39) 


4.  Specific  humidity,  q:  The  specific  humidity  is  defined  as  the 
ratio  of  the  mass  of  water  vapor  present  in  the  total  air  mixture  to 
the  total  mixture.  Thus 
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where  the  denominator  represents  the  total  density  of  the  air,  p. 

It  can  be  easily  shown  that  the  specific  humidity  is  related  to 
the  mixing  ratio  by 


(42) 


5.  Relative  humidity,  h:  The  relative  humidity  is  defined  as 
the  ratio  of  the  mixing  ratio  of  a  sample  of  air  taken  at  some  tempera¬ 
ture  other  than  saturation  to  the  saturation  mixing  ratio  of  the  air  at 
the  same  remperature  t.  Thus 


where  the  subscript  refers  to  the  value  at  saturation, 
humidity  can  also  be  expressed  as 


(43) 

The  relative 
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(44) 


and 
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(45) 
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6.  Dew-point  temperature,  T^:  The  dew-point  temperature  is 
defined  as  the  temperature  at  which  saturation  is  achieved  for  air  cooled 
at  a  constant  pressure  and  mixing  ratio. 

By  substituting  the  relationships  for  the  mixing  ratio  and  the 
specific  humidity  Equation  (26)  can  be  written  as 


d  M 

p  *  p^— <i-q>  [i  +  u(jp)  ]  t 

d  w 


(46) 


Now,  the  ratio  Ma/Mw  is  a  dimensionless  constant  which  can  be  easily 
calculated  (Mj  *  28.966,  =  18)  to  be  1.609.  By  eliminating  w  from 

Equation  (46),  we  get 


p  -  p  |kl+.609q)T 
Md 


(47) 


which  is  known  as  the  equation  of  state  for  air  and  is  usually  expressed 
as 


P  - 


(48) 


where 

T*  -  <l+.609q)T  (49) 

and  is  called  the  virtual  temperature. 


2.4  CONDENSATION  IN  PARTICLE-FREE  AIR 

Even  though  scientists  had  been  seriously  investigating  the  proces¬ 
ses  involved  in  the  formation  of  clouds  since  early  in  the  nineteenth 
century,  they  had  evidently  ignored  the  role  played  by  nuclei  until 
late  in  that  century,  when  the  results  of  John  Aitken's  [Ref.  2]  research 
showed  that  some  sort  of  nuclei  are  involved  in  the  condensation  process. 
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Ait ken's  first  experiments  consisted  of  blowing  steam  into  a  receptacle 
containing  air  which  had  been  filtered  through  cotton  wool.  He  found 
that  a  more  dense  cloud  formed  in  unfiltered  air  than  in  filtered  air 
and  therefore  concluded  that  the  condensation  of  water  vapor  occurs  on 
some  sort  of  nuclei  which  he  referred  to  as  "dust"  particles.  He 
then  constructed  an  expansion  chamber  fitted  with  a  microscope  for 
counting  the  number  of  particles  per  unit  volume  of  air  which  has  come 
to  be  known  as  Aitken's  "dust"  counter.  The  results  of  his  experiments 
with  the  "dust"  counter  indicated  that  below  saturation,  condensation 
occurs  first  on  those  "dust"  (nuclei)  particles  that  are  hygroscopic, 
i.e.,  the  nuclei  that  have  an  affinity  for  water  vapor.  For  hydrophobic, 
i.e.,  non-hygroscopic  nuclei,  some  degree  of  supersaturation  is  required 
before  condensation  can  occur.  Since  saturation  is  defined  as  100 
percent  relative  humidity,  supersaturation  is  defined  as  the  percent 
relative  humidity  above  100,  i.e.,  a  relative  humidity  of  102  percent 
corresponds  to  a  supersaturation  of  2  percent.  Aitken  also  discovered 
that  if  the  air  was  free  of  nuclei,  it  was  still  possible  to  produce 
condensation  but  only  for  an  appreciable  degree  of  supersaturation. 

These  results  indicated  that  condensation  can  occur  in  nuclei- free  air 
by  the  joining  of  sufficient  numbers  of  water  molecules  to  form  an 
initial  nucleus.  However,  the  results  of  Aitken's  work  were  open  to 
argument  because  of  the  possibility  that  the  steam  itself  contained 
condensation  nuclei  and  that  the  filtering  process  removed  only  those 
particles  which  exceeded  a  certain  size.  It  was  for  this  reason, 
together  with  the  knowledge  that  the  data  from  other  researchers  were 
very  incomplete  and  contradictory,  that  C.T.R.  Wilson  [Ref.  58]  in  1897 
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initiated  his  own  research  effort  to  determine  what  degree  of  supersat- 
uration  can  be  obtained  in  moist  air  without  condensation  taking  place. 
Wilson's  experimental  apparatus  was  similar  to  Aitken's  in  that  it  too 
was  an  expansion  chamber  and  was  actually  the  prototype  of  what  was 
later  to  become  known  as  the  Wilson  Cloud  Chamber.  Wilson's  approach 
to  obtain  an  unequivocal  relation  between  degree  of  saturation  and  con¬ 
densation  in  nuclei- free  air  was  to  repeatedly  expand  the  same  sample 
of  air  while  protecting  it  from  contamination.  He  found  that  the  first 
expansion  always  produced  a  fog  if  the  air  had  not  been  allowed  to 
settle  as  completely  as  possible  before  contracting  to  its  original 
volume  so  that  the  nuclei  could  be  removed  by  being  carried  down  into 
the  water  surface  below  by  the  condensation  droplets.  For  each  successive 
expansion  and  removal  of  nuclei,  the  resulting  fogs  became  progressively 
thinner  and  thinner,  the  droplets  becoming  fewer  and  larger  until  the 
fo"  passed  into  a  fine  rain.  Generally,  thereafter,  only  one  more 
expansion  was  needed  to  remove  the  remainder  of  the  nuclei  and  for  any 
greater  expansion,  no  visible  effect  was  produced  until  an  expansion 
ratio  corresponding  to  a  relative  humidity  of  420  percent  (supersatur¬ 
ation  of  320  percent)  was  achieved  at  a  final  temperature  of  -5.8°C. 

The  initial  temperature  had  been  20°C.  With  larger  expansions,  a 
shower  of  drops  was  invariably  produced  until  at  an  expansion  ratio 
corresponding  to  a  supersaturation  of  690  percent,  a  transition  from 
rain-like  to  cloud-like  condensation  occurred.  The  temperature  of  the 
air  at  this  point  of  change  was  -15.8°C.  Wilson  found  that  the  nuclei 
responsible  for  rain-like  condensation  could  not  be  removed  by  filtering 
the  air  through  cotton  wool  nor  did  it  produce  any  effect  upon  the 


expansion  required  for  condensation  to  begin.  Wilson  then  irradiated 
the  filtered  air  in  his  chamber  with  X  rays  and  found  that  no  conden¬ 
sation  occurred  until  exactly  the  same  minimum  supersaturation  levels 
as  in  his  original  experiments  were  reached.  However,  he  did  observe 
that  a  much  denser  cloud  had  formed.  Then,  as  Wilson  increased  the 
supersaturation  above  this  upper  limit,  the  condensation  increased 
with  enormous  rapidity  with  the  result  that  there  was  no  indications  of 
any  limit  to  the  number  of  droplets  that  could  be  produced.  Wilson  cal¬ 
culated  the  size  of  the  smallest  drops  which  would  be  able  to  grow  in 
a  highly  supersaturated  vapor  using  a  relationship  between  the  equilibrium 

(saturation)  vapor  pressure  p  over  the  surface  of  a  droplet  of  pure 

rs 

water  of  radius  r  and  equilibrium  vapor  pressure  pws  over  a  plane  surface 
of  pure  water  at  the  same  temperature  derived  by  Lord  Kelvin  to  be 
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where  r  is  the  surface  tension,  R'  is  the  individual  gas  constant 
(R'»  R/Mw)  and  T  is  the  temperature.  Wilson  found  that  between  the 
limits  of  supersaturation  already  mentioned,  the  value  for  r  varied 
from  6.4xl0~9  to  8.7xlo"8cm,  The  diameter  of  a  water  molecule  is 
4x10  8cm. 

Wilson's  experimental  results  and  calculations  have  been  confirmed 
by  later  researchers  and  the  only  conclusion  that  can  be  drawn  from 
Wilson's  work  is  that  if  in  nuclei-free  and  ion-free  air  the  relative 
humidity  increases  to  800  percent,  water  droplets  occur  spontaneously  when 
water  vapor  molecules  combine  to  form  small  aggregates  of  water  molecules 
which  then  act  as  condensation  nuclei. 
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2.5  THE  NUCLEI  OF  ATMOSPHERIC  CONDENSATION  AND  SUBLIMATION 

While  Wilson's  experiments  in  the  laboratory  demonstrated  that 
enormously  large  supersaturations  are  required  to  produce  cloud-like 
condensation  in  nuclei- free  air,  atmospheric  clouds  are  not  formed  in 
this  manner  because  in  the  atmosphere  both  hygroscopic  and  nonhygroscopic 
particles  and  small  ions  are  always  present  to  act  as  condensation  and 
sublimation  nuclei  when  the  air  is  only  slightly  supersaturated  or  in 
many  occurrences,  even  when  the  relative  humidity  is  less  than  100  per¬ 
cent.  The  atmosphere  contains  a  large  variety  of  particles  in  a  wide 
range  of  sizes  and  concentrations  and  while  scientists  agree  that  some 
sort  of  nuclei  other  than  water  molecules  are  involved  in  the  formation 
of  an  actual  cloud-droplet,  identification  of  the  nature  and  physical 
composition  of  the  actual  nucleus  has  proven  very  difficult  and  is  the 
subject  of  considerable  controversy. 

In  1913,  Wigand  [Ref.  57]  measured  the  number  of  nuclei  in  a 
sample  of  air  with  an  Aitken  counter  before  and  after  he  added  non¬ 
hygroscopic  particles  such  as  coal  dust.  He  found  no  increase  in  the 
nulcei  count  as  a  result  of  the  Increased  dust  content  and  therefore 
concluded  that  nonhygroscopic  dust  particles  are  not  suitable  conden¬ 
sation  nuclei.  Wigand' s  results  and  conclusions  were  verified  in  1926 
by  Boylan  [Ref.  7].  In  1936,  Junge  [Ref.  26]  performed  similar  experi¬ 
ments  with  a  variety  of  nonhygroscopic  particles  such  as  coal  dust,  dirt, 
emery  powder,  ground-up  clay,  household  dust,  and  pulverized  paraffin 
oil.  Using  an  Aitken  counter,  he,  in  contrast  to  Wigand  and  Boylan, 
measured  an  increase  in  nuclei  count  as  a  result  of  increased  dust 
content,  but  only  when  there  was  a  large  concentration  of  particles 
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smaller  than  20  microns  diameter.  Junge  concluded  that  the  ability  of 
nonhygroscopic  particles  to  act  as  condensation  nuclei  depends  on  the 
size  of  the  individual  particle  and  is  generally  independent  of  its 
chemical  composition. 

Many  researchers  have  investigated  the  behavior  of  various  hygros¬ 
copic  particles  both  theoretically  and  experimentally  to  determine  to 
what  extent  they  serve  as  condensation  nuclei.  A  theoretical  expression 
showing  how  the  radius  of  a  droplet  of  solution  varies  with  relative 
humidity  was  derived  by  Kohler  [Ref.  32]  and  is  based  on  Lord  Kelvin's 
derivation  for  a  pure  water  droplet  [Equation  (50)].  A  modified 
version  of  KShler's  expression  first  proposed  in  1936  by  Wright  [Ref.  63] 
was  the  accepted  form  used  for  most  calculations  until  a  slightly  more 
useful  form  was  introduced  by  Mason  [Ref.  40].  Because  scientists 
originally  believed  that  sea-salt  particles  were  the  main  source  of 
cloud-droplet  nuclei  (and  many  still  do),  a  large  amount  of  theoretical 
and  experimental  data  has  been  generated  on  sodium  chloride  (NaCl) 
crystals  and  solutions.  In  1926,  Owens  [Ref.  45]  found  that  NaCl 
crystals  became  deliquescent  at  a  relative  humidity  of  79  percent. 
Experimental  confirmation  of  the  general  theory  and  Owen's  work  was 
published  in  1949  by  a  french  physicist,  H.  Dessens  [Ref.  14].  Dessens 
collected  NaCl  crystals  on  fine  spider  threads  and  measured  the  diameters 
of  the  droplets  formed  about  the  crystals  while  varying  the  relative 
humidity.  He  found  that  the  crystals  became  droplets  of  saturated 
solution  at  a  relative  humidity  of  78  percent.  Again  in  1952,  Junge 
[Ref.  27]  repeated  Dessens*  spider-thread  approach  for  artificially 
produced  nuclei  and  obtained  the  same  results.  He  also  conducted 
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experiments  on  particles  consisting  of  a  mixture  of  insoluble  and 
soluble  matter.  From  his  results  and  Wright's  [Ref.  64]  observation 

t 

that  visibility  is  a  function  of  humidity  for  relative  hvmidites  above 
70  percent,  Junge  suggested  that  in  addition  to  hygroscopic  and  non- 
hygroscopic  nuclei,  the  atmosphere  also  contains  condensation  particles 
consisting  of  insoluble  matter  with  a  thin  covering  of  some  hygroscopic 
substance.  He  called  these  particles  "mixed  nuclei"  in  that  at  relative 
humidities  above  70  percent,  they  react  as  hygroscopic  nuclei,  while 
at  humidities  below  70  percent,  they  behave  as  nonhygroscopic  particles. 

The  chemical  analysis  and  electron  microscopy  of  rime,  frost,  rain 
water,  fog  and  cloud  droplets  made  by  many  other  investigators  tend  to 
verify  Junge' s  theory.  KShler  [Ref.  33]  reported  an  almost  constant 
chloride  content  of  rime  deposited  around  Haldee  Observatory  in  Sweden. 

The  results  of  Woodcock's  [Ref.  60,  61]  research,  which  in  addition  to  lab¬ 
oratory  research,  included  data  collected  from  a  lighthouse  in  Florida, 
and  from  an  aircraft  over  the  mainland  and  sea  around  Hawaii,  indicated 
that  sea-salt  particles  are  present  in  the  atmosphere  up  to  3,000  meters 
and  the  concentration  and  sise  distribution  of  the  particles  did  not 
change  appreciably  for  sea  air  that  had  penetrated  inland  up  to  110  km. 

The  first  attempt  to  identify  actual  condensation  nuclei  with  the  elec¬ 
tron  microscope  was  made  by  Kuroiwa  [Ref.  36]  who  collected  fog  and 
cloud  droplets  on  the  coast  and  in  the  mountains  of  Japan.  Kuroiwa 
found  that  some  droplet  residues  contained  NaCl  crystals  with  traces  of 
MgCl2>  and  MgSO^  while  many  more  contained  nonhygroscopic  particles  such 
as  carbon  granules  and  a  few  consisted  of  "mixed  nuclei."  He  therefore 
concluded  that  combus'  on  products  are  the  most  abundant  source  of 
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condensation  nuclei  with  sea  salt  making  a  secondary  but  numerous 
contribution.  A  similar  study  was  performed  by  Ogiwara  and  Okita 
[Ref.  44]  who  compared  the  residues  of  fog  and  cloud  droplets  with 
those  of  haze  nuclei,  sea-salt  nuclei,  and  artificially  produced  combustion 
products.  Their  results  indicated  that  combustion  products  were  the 
main  source  of  condensation  nuclei  while  the  influence  of  sea- salt 
appeared  to  be  negligible.  The  results  of  a  later  Investigation  by 
Yamamoto  and  Ohtake  [Ref.  40]  agreed  more  with  Kuroiwa's  findings. 

Junge  [Ref.  28]  compared  the  electron-diffraction  patterns  of  conden¬ 
sation  nuclei  gathered  from  actual  fogs  with  the  patterns  of  known 
chemical  compounds.  The  results  of  his  work  indicated  that  particles 
of  amnonium  sulphate  SO^  may  be  an  important  source  of  condensation 

nuclei.  He  also  performed  a  chemical  analysis  of  actual  fog  droplets 
and  found  that  droplets  of  a  radius  greater  than  one  micron  exhibited 

a  high  cloride  content  while  the  smaller  particles  contained  a  higher 
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proportion  of  NH^  and  SO^  .  Two  similar  methods  for  measuring  the 
chemical  composition  of  airborne  particles  from  an  aircraft  were  pro¬ 
posed  by  Seely  [Ref.  50]  and  Lodge  [Ref.  38].  Both  methods  are  based 
on  the  chemical  reaction  between  the  airborne  particles  and  specific 
chemical  reagents  and  both  produced  the  same  results  as  other  methods. 

It  is  generally  believed  that  the  sublimation  process  requires 
some  kind  of  solid  particle  to  act  as  a  nucleus  because  droplets  of 
pure  solution  probably  could  not  induce  sublimation.  Cloud  droplets 
in  the  atmosphere  have  been  observed  in  the  liguid  state  at  temperatures 


33 


well  below  OX  and  experimental  resultB  indicate  ice  crystals  form  more 
commonly  by  the  freezing  of  these  supercooled  droplets  than  by  direct 
sublimation. 

While  the  exact  nature  and  chemical  composition  of  the  particles 
which  serve  as  condensation  and  sublimation  nuclei  is  still  the  subject 
of  considerable  controversy,  scientists  generally  agree  that  atmospheric 
nuclei  can  originate  in  the  following  ways:  (a)  by  condensation  and 
sublimation  of  the  vapors  produced  from  industrial  and  domestic  combus¬ 
tion  processes;  (b)  by  condensation  and  sublimation  of  the  vapors 
produced  from  gaseous  reactions  in  the  atmosphere;  (c)  by  the  disinte¬ 
gration  of  matter,  such  as  the  formation  of  mineral  dust  and  sea  spray; 
and  (d)  by  the  coagulation  of  nuclei.  The  major  unanswered  question  is 
which  process  is  the  main  source  of  condensation  and  sublimation  nuclei. 
2.6  THE  GROWTH  AND  SIZE  DISTRIBUTION  OF  CLOUD  DROPLETS 

As  we  have  discussed  earlier,  the  meteorlogical  variables  used  for 
describing  the  thermodynamic  conditions  of  the  atmosphere  are  based 
measurements  taken  in  a  closed  volume  of  air  over  a  plane  surface  of 
pure  water,  and  while  these  measurements  are  sufficiently  accurate  for 
most  meteorlogical  purposes  such  as  weather  forecasting,  a  study  of  the 
processes  Involved  in  droplet  growth  must  include  the  effects  of  a 
spherical  surface  of  a  solution  of  impure  water.  Wilson  demonstrated 
that  in  nuclei  and  ion-free  air,  condensation  nuclei  formed  from  samll 
aggregates  of  water  molecules  occur  spontaneously  at  extremely  high 
supersaturations.  These  nuclei  will  remain  in  the  liquid  state  and  grow 
into  water  droplets  only  if  they  attain  a  critical  size  determined  by 
the  relationship  expressed  in  Equation  (SO).  We  rearrange  Equation  (50) 


36 


in  a  slightly  different  form,  i.e.. 
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to  emphasize  the  relationship  that  must  be  satisfied  in  order  that  a 
condensation  nuclei  survive  and  not  evaporate.  This  expression  for  the 
equilibrium  vapor  pressure  over  a  curved  liquid  surface  of  pure  water 
was  originally  derived  by  Lord  Kelvin  [Ref.  29].  An  atmospheric  droplet 
is  never  composed  entirely  of  pure  water  and  an  expression  similar  to 
Lord  Kelvin's  for  the  equilibrium  vapor  pressure  over  a  curved  liquid 
surface  of  pure  solution  was  first  derived  by  Kohler  [Ref.  32]  in  1921. 
While  several  variations  of  Kdhler's  expression  have  been  proposed, 
only  Wright's  [Ref.  63]  and  Mason's  [Ref.  40]  versions  have  found  general 
acceptance  by  most  scientists.  Wright's  expression  for  the  equilibrium 
vapor  pressure  over  a  droplet  of  solution  may  be  written 


exp(- 


rs 


ws 


2tM 

Jr-)  -  imM 

p  RTr  w 

w _ 

V3irr3p  M 
w 


(52) 


where  m  is  the  mass  and  M  is  the  molecular  weight  of  the  solute 
(nucleus)  and  i  is  known  as  the  van't  Hoff  factor  which  is  a  function 
of  the  chemical  composition  and  degree  of  dissociation  of  the  solute. 
Mason's  expression,  which  is  more  comprehensive  and  lengthy,  will  not  be 
presented  here.  A  complete  and  vigorous  derivation  is  given  in  Mason 
[Ref.  40], 

If  a  nucleus  droplet  attains  or  exceeds  *-he  critical  radius  necessary 
to  ensure  that  it  will  not  evaporate  during  the  early  stages  of  its 
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formation,  then  the  droplet  will  continue  to  grow  if  a  sufficient  degree 
ot  supersatutation  is  maintained.  The  growth  from  a  nucleus  droplet 
to  a  cloud  droplet  is  accomplished  by  the  diffusion  of  water  vapor  to, 
and  its  subsequent  condensation  upon,  the  surface  of  the  droplet,  and 
by  collision  and  coalescence  with  other  droplets. 

From  Fick's  Law  of  diffusion  the  flux,  F,  of  water  vapor  in  a  radial 
direction,  R,  toward  the  droplet,  is  given  by 


F 


(53) 


and  the  growth  of  the  droolet  of  diffusion  is  therefore 

,  dP  dp 

* DA  dr " 47rR  D  dir  (54) 

where  ~  is  the  time  rate  increase  of  droplet  mass,  A  =  4tcr2  is  the  sur¬ 
face  area  of  the  droplet,  D  is  the  coefficient  of  diffusion  of  water 
dP  w 

vapor  in  air  and  is  the  radial  gradient  of  vapor  density  which  is 

assumed  to  be  symmetrical  around  the  drop.  If  we  assume  that  steady- 

dm 

state  conditions  exist,  then  —  is  a  constant  and  Equation  (54)  can  be 
integrated  with  respect  to  the  radial  distance  from  the  surface  of  the 
droplet  where  the  vapor  density  is  to  infinity  where  it  is  P^. 
Therefore, 
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which  gives 
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Now,  we  can  also  express  the  time  rate  change  of  droplet  mass  as 
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dm  .  ,  dr 

dl  "  ^P»r2  ~ 


where  p  is  the  density  of  the  droplet. 


dt 


(57) 


Equating  Equation  (56)  and  Equation  (57)  gives 


rdr  ■  r  <‘V‘>r)dt 


(58) 


Integrating  the  resulting  expression,  assuming  that  as  the  radius  of 
the  droplet  varies  from  its  critical  value  r^  to  r,  the  time  goes  from 
t  '  0  to  t  =  t  gives 


r2  ■  r2  +  ~  (p  -p  )t 
c  p  '  w  r 

This  expression  can  be  written  in  terms  of  vapor  pressure  as 
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where  R  is  the  universal  gas  constant,  pr  and  are  respectively  the 
vapor  pressure  and  temperature  at  the  surface  of  the  droplet  and  and 
T  are  the  ambient  vapor  pressure  and  temperature.  We  see  that  a  neces¬ 
sary  condition  for  droplet  growth  (condensation)  is  that  the  term  in 
parentheses  be  positive.  A  negative  term  signifies  a  state  of  evapor¬ 
ation.  And  since  the  temperature  decreases  with  increasing  distance 
from  the  droplet,  the  necessary  condition  for  droplet  growth  is  that 
the  ambient  vapor  pressure  be  greater  than  ie  pressure  at  the  surface 
of  the  droplet. 

In  addition  to  the  other  factors  which  determine  the  growth  rate 
of  a  droplet,  the  condensation  of  water  vapor  to  or  the  evaporation 
from  a  droplet  involves  the  addition  or  removal  of  heat  from  the 
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surrounding  air  and  is  given  by 


(61) 


dQ  da 
dt  "  L  dt 
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Where  ^  is  the  rate  of  latent  heat  transfer,  L  is  the  latent  heat  of 

Jin 

condensation  or  evaporation,  K  is  the  thermal  conductivity  of  air,  ~ 

qr 

is  the  temperature  gradient  and  the  choice  of  sign  depends  on  whether 
heat  is  being  added  or  removed  from  the  air.  This  expression  may  be 
integrated  in  a  similar  manner  to  Equation  (54)  and  for  the  condensation 
process  (negative  sign)  we  obtain 

L  ^  -  4ffKr(Tr-T)  (62) 

and  as  before,  by  substitution  of  Equation  (57)  we  obtain 

r  ft  ■  TTt  <V»  <«> 


which  can  be  integrated  in  a  similar  manner  as  Equation  (58). 

Now  we  can  see  from  Equation  (58)  that  for  a  given  gradient  of 
djr 

vapor  density,  is  inversely  proportional  to  r  and  from  Equation  (62) 
we  see  that  for  a  given  increment  of  radius  the  rate  of  transfer  of 
latent  heat  increases  as  the  radius  of  the  droplet  increases,  i.e., 

^  *  r.  Therefore,  since  is  assumed  constant,  an  increase  in  the 
radius  of  *V  droplet  produces  a  decrease  in  the  vapor  gradient  until 
the  excess  water  vapor  is  exhausted  and  an  equilibrium  condition  is 
reached.  It  is  for  these  reasons  that  most  researchers  agree  that  the 
condensation  process  alone  is  unlikely  to  produce  droplets  of  radii 
'>30  microns.  Actual  measurements  of  droplet  sizes  have  shown  that 
while  large  droplets  with  radii  up  to  100  microns  (borderline  between 
cloud  drops  and  rain  drops)  exist,  most  cloud  and  fog  droplets  have 
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radii  of  less  Chan  30  microns  with  peak  concentrations  somewhere  between 
1  and  10  microns.  Condensation  on  very  large  nuclei  may  produce  a  few 
of  the  larger  droplets;  however,  growth  by  condensation  for  the  larger 
droplets  requires  several  hours.  It  is  believed  that  the  growth  to 
larger  droplets  and  raindrops  is  accomplished  by  the  collision  and 
coalescence  of  droplets  moving  relative  to  each  other.  Theoretical 
calculations  show  that  cloud  and  fog  droplets  grow  more  rapidly  by  con¬ 
densation  during  the  early  stages  of  development  and  the  droplets  must 
attain  a  radius  of  at  least  20  microns  before  the  coalescence  process 
begins  to  predominate.  Experimental  verification  has  been  obtained  by 
several  researchers  and  an  excellent  discussion  of  their  work  is  presented 
in  Mason  [Ref.  40] . 

Much  effort  has  been  devoted  to  developing  experimental  techniques 
and  equipment  for  measuring  the  size  and  size  distribution  of  atmospheric 
and  laboratory  simulated  fog  and  cloud  droplets.  The  development, 
however,  of  reliable,  accurate,  and  rapid  labor  saving  instruments  to 
measure  droplet  size  and  distribution,  especially  in  the  1-5  microns 
range  is  an  urgent  requirement  for  cloud  physics  research.  Notably 
lacking,  also,  are  suitable  techniques  and  equipment  for  use  in  airborne 
research  such  as  from  an  aircraft.  The  newest  instrument,  a  laser 
disdrometer,  in  which  holography  is  used  for  measuring  size  and  dis¬ 
tribution  of  droplets  within  a  given  volume  is  limited  in  its  application 
and  has  not  been  adapted  to  airborne  use.  In  addition,  the  illuminated 
volume  is  small  which  requires  a  large  number  of  holographs  to  obtain  a 
representative  sample  of  the  size  distribution  and  an  ever  tedious 
drop-by-drop  measurement  is  the  only  reliable  technique  presently 
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Available  for  sorting  the  droplets.  The  most  reliable  data  on  actual 
droplet  distributions  were  obtained  primarily  by  simple  methods  where 
the  droplets  are  caught  on  a  specially  prepared  surface  and  the  diameter 
of  the  droplets  or  their  impressions  are  measured  under  the  microscope. 

A  summary  of  some  of  the  more  important  parameters  characterising  the 
size  distribution  of  droplets  in  different  types  of  clouds  along  with 
the  appropriate  reference  is  shown  in  Table  1.  The  parameters  used  in 
Table  1  are  defined  as  follows: 

^min  *s  Minimum  diameter  of  the  droplets 

d  is  the  mean  diameter  (the  sum  of  all  droplet  diameters 
x“  divided  by  the  total  number  of  droplets) 

d^  is  the  mode  diameter  (the  diameter  corresponding  to  the 
maximum  number  of  droplets) 

d  is  the  predominant  volume  diameter  (the  diameter  corres- 
v  ponding  to  the  maximum  volume  of  water) 

d  is  the  maximum  diameter  of  the  droplets 
max 

Much  effort  has  been  devoted  to  formulating  an  analytical  expression 
that  would  provide  a  best  fit  for  the  many  different  empirical  distri¬ 
butions  obtained  from  the  numerous  droplet  size  measurement  programs. 

Best  [Ref.  5]  examined  the  experimental  data  of  several  researchers 
[Refs.  15,  16,  17,  18]  and  found  that  most  of  the  observed  droplet  size 

distributions  could  be  represented  by  the  expression 

b 

1  -  F  •  exp  [(~p)  ]  (64) 

where  F  is  the  fraction  of  liquid  water  in  the  air  composed  of  drops  of 
diameter  less  than  x,  a  is  a  numerical  factor  dependent  on  the  liquid 
water  content  and  b  is  a  constant  with  a  mean  value  of  3.27.  One  of  the 
first  universally  accepted  distributions  was  formulated  by  Junge  [Ref.  28] 
and  is  given  by  the  expression 
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Table  1 

Characteristics  of  Natural  Fog  and  Cloud  Droplets 
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N(r)  ■  cr 

Where  N(r)  is  the  number  of  particles  per  cubic  centimeter  of  size  r 
and  c  and  p  are  numerical  constants.  A  more  recent  and  widely  accepted 
distribution  proposed  by  Deirmendjian  [Ref.  13]  is  a  generalization  of 
the  distribution  first  proposed  by  Khrgian  and  Mazin  which  may  be 
written 

n(r)  ■  ar°e  (66) 

where  n(r)  is  the  volume  concentration  at  the  radius  r  and  a,  a,  b  and 
y  are  positive  constants.  While  many  other  distributions  have  been 
proposed,  the  three  distributions  quoted  here  represent  the  most  widely 
accepted  forms.  In  general,  the  attempts  of  the  many  researchers  to 
analytically  describe  the  droplet  size  distribution  have  resulted  in 
the  conclusion  that  atmospheric  aerosol  distributions  exist  in  con¬ 
siderable  variety. 


SECTION  III 


OPTICAL  PROPERTIES  OF  ATMOSPHERIC  AEROSOLS 

3.1  INTRODUCTION 

The  propagation  of  electromagnetic  energy  through  the  atmosphere 
has  been  studied  for  many  decades,  and  the  quantity  of  available  litera¬ 
ture  on  the  subject  is  immense.  The  subject  has  received  the  attention 
of  people  in  astronomy,  communications,  radar,  meteorology,  physics,  j 
physical  chemistry,  optics,  spectroscopy  and  many  other  fields  of  endeavor. 
In  some  cases,  this  high  degree  of  documentation  has  been  accompanied 
by  ever-increasing  precision.  Yet  in  the  case  of  atmospheric  aerosols 
there  is  still  controversy  as  to  whether  salt  crystals,  sulphuric  acid 
or  nitrous  oxide  plays  the  dominant  role  in  the  nulceation  of  clouds 
and  fogs. 

As  we  discussed  in  the  previous  chapter  the  atmosphere  is  composed 
of  a  mixture  of  gases  and  water  vapor  in  which  are  susnended  a  variety 
of  hygroscopic  and  hydrophobic  substances  some  of  which  may  act  as  con¬ 
densation  nuclei  in  the  formation  of  a  cloud.  These  chemical  and  particu¬ 
late  constituents  of  the  atmosphere  attenuate  electromagnetic  radiation 
resulting  in  a  loss  of  intensity  greater  than  that  experienced  in  free 
space.  The  attenuation  of  the  propagating  beam  Is  caused  by  two  Independent 
processes :  (1>  absorption,  whereby  energy  is  removed  from  the  beam  by 

converting  the  radiation  to  some  other  form  of  energy,  and  (2)  scattering, 
whereby  some  of  the  beam  energy  Is  redirected  from  its  origins]  direction 
of  propagation.  The  gases  and  water  vapor  cause  the  radiation  to  be 
absorbed  and  the  scattering  of  radiation  is  caused  by  non  ini f oral ties 
which  may  be  due  either  to  the  presence  of  particulate  matter,  to 
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macroscopic  variation  in  the  index  of  refraction  (atmospheric  turbulence), 
or  to  fluctuations  on  a  microscopic  scale  of  the  number  of  molecules  per 
unit  volume  (molecular  scattering).  The  total  attenuation,  then,  is 
equal  to  the  sum  of  the  absorbed  and  scattered  radiation  and  is  a  function 
of  the  wavelength  of  the  propagating  energy. 

When  a  collimated  beam  of  monochromatic  electromagnetic  radiation 
of  intensity  1  propagates  through  the  atmosphere,  a  small  increment  dl 
of  the  intensity  will  be  attenuated  in  passing  through  an  infinitesimal 
thickness,  dx,  so  we  may  write 

dl  -  -oldx  (67) 

where  °  is  a  measure  of  the  loss  per  unit  length  and  is  called  the 
extinction  (attenuation)  coefficient.  For  a  wave  propagating  over  a 
finite  distance,  x.  Equation  (67)  can  be  integrated  to  give 

I  -  loe‘OX  (68) 

where  Iq  is  the  incident  Intensity  and  I  is  the  intensity  at  the  distance 
x. 

If  we  assume  that  the  attenuation  of  the  propagating  beam  is  due 
to  the  two  Independent  processes  of  absorption  and  scattering,  then  the 

extinction  coefficient  is  the  sum  of  an  absorption  coefficient,  o  , 

ft 

and  a  scattering  coefficient,  o^,  and  Equation  (68)  becomes 

I  ■  I  e'<a.  +  %>x  (69) 

o 

Equations  (68)  and  (69)  are  often  referred  to  as  Bouguer's  Law 
named  after  the  French  scientist  who  in  1760  formulated  the  law  by 
experiment . 

Both  the  coefficients  o#  and  o#  are  functions  of  the  wavelength 
of  the  propagating  wave,  and  a  medium  is  arbitrarily  called  transparent 
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if  Che  extinction  distance  (defined  to  be  1/a)  is  large  compared  to  the 

path  length  of  interest.  In  the  hard  X-ray  region  — .01  to  .1  Angstroms 
0 

(A) —  the  attenuation  is  due  both  to  molecular  absorption  and  molecular 
scattering  (because  of  the  short  wavelengths  the  air  molecules  themselves 
act  as  scatterers)  and  the  extinction  distance  is  greater  than  30  meters. 
As  the  wavelength  increases,  the  extinction  distance  decreases  until 

O 

in  the  2  to  20A  region  it  is  less  than  10  centimeters.  As  the  wave¬ 
length  continues  to  increase  from  the  X-ray  region  to  the  far  ultra- 
0 

violet  (around  100A)  the  extinction  distance,  due  mostly  to  absorption, 
is  still  less  than  10  centimeters.  At  2100A  a  window  appears  which  for 
pure  oxygen  and  nitrogen  extends  out  into  the  far  infrared  (IR) . 

The  extinction  distance  would  be  greater  than  300  meters.  However, 
water  vapor  produces  strong  absorption  bands  in  the  regions  near  l.S, 

1.9,  2.7  and  6.3  microns  (u)  while  carbon  dioxide  produces  strong 
absorption  bands  at  2.7  and  4.3y.  The  presence  of  both  of  these  gases 
in  the  atmosphere  reduces  the  extinction  distance  to  less  than  one 
meter  in  these  regions.  As  the  wavelength  increases  through  the  far  IR 
into  the  millimeter  wavelength  region,  the  attenuation  is  due  primarily 
to  the  scattering  from  particulate  matter  and  decreases  steadily 
until  below  X-bsnd  the  atmospheric  attenuation  is  essentially  negligible. 

In  this  Investigation  we  will  restrict  ourselves  to  wavelengths 
within  the  regions  of  major  atmospheric  windows  where  the  extinction 
length  is  much  greater  than  the  path  lengths  of  interest  and  where  two  of 
the  scattering  mechanisms,  atmospheric  turbulence  and  molecular  scat¬ 
tering,  are  essentially  negligible.  In  atmospheric  turbulence,  the 
sise  of  the  disturbing  elements  is  several  orders  of  magnitude  larger  then 
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macroscopic  variation  in  the  index  of  refraction  (atmospheric  turbulence), 
or  to  fluctuations  on  a  microscopic  scale  of  the  number  of  molecules  per 
unit  volume  (mole:ul.ar  scattering).  The  total  attenuation,  then,  is 
equal  to  the  sum  of  the  absorbed  and  scattered  radiation  and  is  a  function 
of  the  wavelength  of  the  propagating  energy. 

When  a  collimated  beam  of  monochromatic  electromagnetic  radiation 
of  intensity  I  propagates  through  the  atmosphere,  a  small  increment  dl 
of  the  intensity  will  be  attenuated  in  passing  through  an  infinitesimal 
thickness,  ax,  so  we  may  write 

dl  *  -oldx  (87) 

where  cr  is  a  measure  of  the  loss  per  unit  length  and  is  called  the 
extinction  (attenuation)  coefficient.  For  a  wave  propagating  over  a 
finite  distance,  :c,  Equation  (67)  can  be  integrated  to  give 

I  »  I  e  *  v'  (68) 

o 

where  Iq  is  the  incident  intensity  and  I  is  the  intensity  at  the  distance 
x . 

If  we  assume  that  the  attenuation  of  the  propagating  beam  is  due 

to  the  two  Independent  processes  of  absorption  and  scattering,  then  the 

extinction  coefficient  is  the  sum  of  an  absorption  coefficient,  o^, 

and  a  scattering  coefficient,  a  ,  and  Equation  (68)  becomes 

s 

T  T  "  "b  0  )x  f  (LQ\ 

I  ■  I  e  a  s  (b9) 

o 

Equations  (68)  and  (69)  are  often  referred  to  a3  Bouguer's  Law 
named  after  the  French  scientist  who  in  1760  formulated  the  law  by 
experiment. 

Both  the  coefficients  c  and  a  are  functions  of  the  wavelength 

3  S 

of  the  propagating  wave,  and  a  medium  Is  arbitrarily  called  transparent 
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Scattering  is  probably  responsible  for  the  many  attempts  to  formulate  a 
similar  relationship  between  the  scattering  coefficient  and  wavelength 
for  atmospheric  aerosols.  However,  the  derivation  of  a  theory  describing 
the  scattering  mechanism  in  aerosols  such  as  clouds  and  fogs  is  compli¬ 
cated  by  the  many  variables  such  as  particle  composition,  concentration, 
size  distribution,  etc.  Therefore,  a  straightforward  relationship  has 
not  been  found  as  in  the  case  of  Rayleigh  scattering. 

In  1908  a  German  meteorologist,  Gustav  Mle  [Ref.  42],  derived  in 
terms  of  electromagnetic  wave  theory  a  solution  for  the  electric  and 
magnetic  vectors  describing  the  radiation  scattered  by  a  homogeneous 
sphere  of  any  diameter  situated  in  a  homogeneous  medium.  While  his 
work  was  originally  in  the  field  of  colloidal  chemistry,  it  is  applic¬ 
able  to  many  other  problems  and  Stratton  and  Houghton  [Ref.  52]  were  the 
first  researchers  to  apply  the  Mie  Theory  of  Scattering  to  the  problem 
of  optical  transmission  in  atmospheric  aerosols. 

3.2  THE  MIE  THEORY  OF  SCATTERING 

In  1908  Gustav  Mie  [Ref.  42]  derived  the  first  rigorous  solution  for 
the  diffraction  of  a  linearly  polarized  monochromatic  plane  wave  by  a 
single  homogeneous  sphere  of  any  description  situated  in  a  homogeneous 
medium.  Mie  approached  the  problem  in  terms  of  the  electromagnetic  wave 
theory  of  radiation  in  which  a  periodic  wave  incident  on  a  homogeneous 
sphere  of  arbitrary  size  and  composition  induces  a  forced  oscillation  of 
the  free  and  bound  charges.  These  forced  oscillations  themselves  generate 
a  secondary  field  both  inside  and  outside  the  sphere  where  the  resultant 
field  outside  the  sphere  is  the  sum  of  the  incident  and  secondary  fields. 
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In  1909  D«bye  [Ref.  12]  published  a  similar  solution  of  the  same  prob¬ 
lem  where  he  used  a  pair  of  potential  functions  that  have  come  to  be 
known  as  Debye’s  potentials.  While  Mie's  original  work  was  concerned 
with  the  passage  of  light  through  a  colloidal  suspension  of  metallic 
particles,  the  Mle  Theory  of  Scattering  is  of  great  practical  value. 

Since  then  many  researchers  have  applied  the  theory  to  a  variety 
of  problems.  Stratton  and  Houghton  [Ref.  52]  appear  to  be  the  first 
writers  to  apply  the  theory  to  the  transmission  of  light  through 
clouds  and  fogs. 

A  thorough  treatment  of  Mie's  theory  is  beyond  the  scope  of  this 
paper  and  we  shall  therefore  restrict  ourselves  to  those  parts  of  the 
theory  that  are  directly  applicable  to  our  investigation.  The  reader 
is  referred  to  Stratton  [Ref.  53]  or  Born  and  Wolf  [Ref.  6]  for  a 
vigorous  mathematical  treatment  and  to  van  de  Hulst  [Ref.  54]  for  a  clear 
and  comprehensive  discussion. 

Mie  treated  the  scattering  of  light  by  a  homogeneous  sphere  in 
terms  of  Maxwell's  field  relations  and  the  corresponding  wave  equations. 
We  know  that  an  electromagnetic  field  which  is  a  function  of  both  space 
and  time  cannot,  in  general,  be  derived  solely  from  a  scalar  function. 

For  a  solution  to  this  problem  the  integration  of  Maxwell's  equations 
can  be  reduced  to  deriving  a  vector  and  a  scalar  potential  that  both 
satisfy  the  same  differential  equation  (commonly  referred  to  as  the 
inhomogeneous  wave  equation)  as  that  obeyed  by  the  fields.  On  the 
other  hand  Hertz  [Ref.  20]  has  shown  that  an  electromagnetic  field  can 
be  expressed  in  terms  of  a  single  vector  function.  We  shall  show  that 
an  electromagnetic  field  within  an  isotropic,  homogeneous,  charge-free 
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conducting  medium  can  be  resolved  into  two  partial  fields,  each  derivable 
from  separate  Hertz  vectors  each  of  which  satisfy  the  vector  wave 
equation. 

We  begin  by  expressing  Maxwell's  time  varying  field  relations  for 
a  homogeneous.  Isotropic,  charge-free,  conducting  medium  in  the  following 
form 


7  X  H  -  (a  +  e— )E 

(70) 

5  x  *  -  -*n 

(71) 

$4  -  o 

(72) 

— ► 

7-E  ■  0 

(73) 

where 

H  »  magnetic  field  intensity 
E  ■  electric  field  intensity 
a  =>  conductivity  of  the  medium 
e  «  permittivity  of  the  medium 
p  *  permeability  of  the  medium 

Since  the  if  field  is  solenoidal  it  can  be  represented  by  the  curl 
of  another  vector  /t,  i.e. 

H  -  $  X  t  (74) 

and  it  can  be  shown  that  the  corresponding  £  field  is  given  by 

E  -  4*  -  ||  (75) 

Now  since  the  E  field  is  also  solenoidal  an  alternate  solution  can 
be  derived  in  terms  of  a  vector  potential  A  and  a  scalar  potential  . 

I  -  -  -  $  X  1*  (76) 

£ 
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(77) 


Now,  if  we  let 


-*■  *  3A  a  * 

"If  “IA 


A  *  (a  +  e-^)n 


(78) 


where  n  is  another  vector  potential  then  Equations  (74)  and  (75)  become 

3 


H  -  (a  +  e~)  Vxl 


E  -  -  7^  -  (o  +  e~)  || 


(79) 


(80) 


Now,  since  the  choice  of  the  scalar  potential  ip  is  arbitrary  as  long  as 
it  satisfies  the  wave  equation,  we  can  define  it  as 


\p  =  -V»JI 


(81) 


Equation  (80)  now  becomes 


I  -  7^4  -  (0  +  c~)  || 


(82) 


Upon  substitution  of  Equations  (79)  and  (82)  into  Equation  (70)  we 
obtain  the  vector  equation 


ulxl  -  $v4  +  po||  +  pe—  -  0 


at 


(83) 


3t^ 


every  solution  of  which  determines  an  electromagnetic  field  given  by 


H  ■  (0  +  e—)  $  X  f 


(84) 


I  -  $  X  $  X  it  (85) 

An  alternate  solution  satisfying  Equations  (76)  and  (77)  can  be 

■fit  * 

constructed  in  terms  of  A  and  ip  where 


jr* 

'3t 


(86) 
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i 


(87) 


\p  »  -v*n 

Consequently, 

H  ■  ^  X  I  X  f  (88) 

i  x  if  <89) 

The  general  solution,  then,  for  a  time-varying  electromagnetic 
field  in  a  homogeneous,  isotropic,  charge-free,  conducting  medium  can 
be  expressed  as  the  sum  of  two  partial  fields,  each  derivable  from 
separate  Hertz  vectors  (?  and  II*) . 

Thus, 

I  -  ?  X  ?  X  «  -  y?  X  |2  (90) 

h  «  5x  ?  x  I*  +  h  (o  +  £— -)n  (91) 

3t 

Since  the  scattered  fields  represent  outgoing  spherical  waves, 

Mie  chose  to  represent  Maxwell's  equations,  the  incident  radiation, 
and  the  scattered  radiation  in  terms  of  the  spherical  coordinates,  r, 

9,  and  <p  defined  by 

x  ■  r  cos9  (92) 

y  ■  r  sin9cos<J>  (93) 

z  ■  r  sin8sin<j>  (94) 

where  x,  y,  and  z  represent  a  rectangular  coordinate  system  with  origin 
at  the  center  of  the  sphere.  The  z-direction  is  taken  as  the  direction 
of  propagation  of  the  incident  wave  with  the  x-direction  in  the  direction 
of  its  electric  field  vector.  In  our  derivation  of  Mie's  formulas  we 
modify  Mie's  geometry  slightly  and  redefine  the  angle  6  to  be  the  angle 
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between  the  direction  of  propagation  of  the  incident  wave  and  the 

scattered  wave.  With  this  definition  the  angle  6  is  universally 

referred  to  as  the  scattering  angle  and  the  corresponding  curvilinear 

coordinate  system  is  defined  as 

x  *  r  sin0cos<|>  (95) 

y  ■  r  sinGsin^  (96) 

z  -  r  cos6  (97) 

where  the  plane  of  scattering  is  defined  as  the  plane  containing  the 

axis  of  propagation  and  the  radius  vector  from  the  center  of  the  sphere 

to  a  point  of  observation.  If  we  define  the  plane  of  oscillation  as  • 

the  plane  containing  the  electric  field  vector  of  the  incident  wave 

and  the  axis  of  propagation,  then  the  angle  <f>  represents  the  angle 

between  this  plane  and  the  plane  of  scattering  (See' Fig.  1)'. 

.  -► 

In  spherical  coordinates,  the  curl  of  any  vector  A  is  given  by 


$ 


X  A 


t 


rsin9 


9 (sin0 A, )  9A0 


.  3A  '  9(rA.)  _ 

+  ( 2  [ — i - £  _  —  *  t.  }  1  a 

‘  plsine  -3*  9r  J»  fl 


.  9(rA)  9A 

*  \ 


(98) 
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Geometry  for  Mie  Scattering  Problem 
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(102) 


1  3<rHfl>  3Hr 

7t-5T-  -  if’  ■  <"  + 

,  3(rsin0E.)  3(rE.) 

_A_  f. - _i - 5^-]  -  -IttuH 

r2sin0  30  3*  r 


,  3(rE  )  3E 

v'-ir-  -  if  i  ■ 


*  — 1 ' 

3E  3(rsin0E  ) 

Tib'lf  -  — IT-4"  ‘  _1“uHe  <104> 

,  3<rEe>  3Er 

i>-ir-  -  ifi  ■  (105) 

,  3(r2H  )  .  3(sin0H.)  .  3H 

_L - E_  +  — i - —  +  - i  .  o  (io6) 

^2  3r  rsin0  30  rsin0  3 <(> 

,  3(r2E  )  .  3 (sin0Efl)  .  3E 

JL _ L_  +  — 1 - §_  +  —I - i  .  o  (107) 

o  3r  rsln@  30  rsin@  34 

r* 

Mle  represented  the  solution  of  these  equations  as  a  superposition 
of  two  linearly  independent  fields  such  that 

Er  *  0,  Hr  -  0  (108) 


1 

3(r2Hr) 

3(sin0HQ) 

+  1 

3H. 

r2 

3r 

rsin6 

30 

rsin0 

3  4 

1 

3(r2Ef) 

+  1 

3 (sin0E0) 

+  1 

!!i 

3r 

rsln@ 

30 

rsinO 

34 

Er  -  0,  Hr  *  0 

Substituting  Equation  (108)  into  Equations  (101)  and  (102)  gives 

,  3(rHj 


(o  +  iwe)E0  -  -  7  — ^ 


1  3<rV 

(°  +  iwe>%  "  x  57” 


(109) 


(110) 


(111) 


Upon  substituting  these  equations  in  Equations  (104)  and  (105)  we  obtain 

I  0  ^  -  3  <rHe>  4-  k2rH.  (112) 


aine  34  3r2  e 

3E  32 (rH . ) 

-(o  +  iwe)  -TT-  •  — ' — +  k2rH 
36  0r2  4 


(113) 
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where  k2  ■  w2ue  -iuiyo 

Similarly,  for  the  relations  of  Equation  (109)  a  solution  to 

Maxwell's  equations  can  be  represented  by 

l  3<rV 
“"‘"e  ’  7  S?- 

1  3<rEfl> 

^  "  -  7  Tr 


3H  32(rEa) 

.  JXSL  .  -  6  +  k2rE 

sine  3$  ar2  k  Ee 


(114) 


(115) 


(116) 


3H  32(rE.) 

iww  -r~  - - *-  +  k2rE 

39  a-2  ♦ 


(117) 


It  is  not  difficult  to  see  that  the  relationshipa  expressed  in 
Equations  (108)  and  (109)  define  a  complementary  pair  of  Hertz  potentials 
ft  and  ft*  such  that 


nr  +  o,  n0  -  *  o 


(118) 


nr  *  °*  ne  *  %  -  0 


(119) 


For  II  0,  n  ■  n  ■  0  we  epply  the  vector  expansions  of  Equation  (98) 
r  ®  ♦ 

and  (99)  to  Equations  (84)  end  (85)  to  obtain 


jn 

»<»*-  ir> 


r2slnO 


r2sin2e  3#2 


(120) 


Be  7  3r39 


(121) 


"  rsin9  3r34 


(122) 


(123) 


a  ,  i°  t  i.y£)  !jl 
0  rsine  a  $ 


(124) 


_  (p  +_>ic) 

♦  r  39 

If  we  now  substitute  from  these  equations  into  Equations  (112)  and 
(113)  we  obtain 


(125) 


.  a2n  . 
i_  t — L  +  — - — 

3<^  ar2  r2sin0 


an 

a(sin9-rr^)  .  a2n 

- + - i - -  +  k2nj  -  0  (126) 

30  r2sin20  3*2  r 


an 


.  a2n 

Is  l — "  + 

30  ar2  r2sine 


r  .  i  3<8ln&-alT> 

ae 


r'  .  a2n 

+ - - - -  +  k2n  ]  -  0  (127) 

r2sin20  3$2  r 


The  vanishing  of  the  6  and  $  derivatives  may  be  satisfied  by  equating 
the  expression  within  the  brackets  to  zero,  which  gives 

a2n 


an 

a(sine-^> 


ae 


i  »2n 

- i - L  +  k2n  «  o 

r28in20  3$2  r 


(128) 


3r2  r2sin0 

Comparing  this  expression  with  Equation  (120)  we  see  that  may  be 
written  as 

a2n, 

e  -  — -  +  k2n 
r  dr2  r 

In  a  similar  manner  for  n*  #  0»  n*  »  fl*  ■  0  we  apply  the  vector 

r  s  s 

expansions  of  Equations  (98)  and  (99)  to  Equations  (88)  and  (89)  to  obtain 


(129) 


Er  •  0 


(130) 


.  an 

_  1mm  r 

S  "  "  rsine  ae 


(131) 


57 


(132) 


H  -  - 
r  o 


r*-sin0 


311 
_£ 
r  30 
* 
an 

9<sin97#> 

30 


3  2 II 


r2sin20  3<t>2 


H„  - 


.  32n* 

1  r 


0  r  3r3G 


(133) 


(134) 


D?-n 


H,  - 


<j>  rsinO  3r3(|) 

and  upon  substitution  from  these  expressions  into  Equations  (116) 
(117)  we  obtain 


(135) 

and 


32JI 


?n 


3  (sine- 


30  ' 


o  * 

32n 


3r2  r2sin0 
Consequently, 


36 


- ± - L  +  k2n* 

r2sin20  dp2  C 


H  - 

r 


s2n' 


3  r2 


+  k2n 


Now  let  nr  ■  rip,  then  Equation  (128)  becomes 


1  a2(r»)  |  1 _ 

r  3r2  r2sin0 


+  — J—  +  k2*  -  0 

r2sin20  3<{>2 


A 

and  for  n  *  rip  Equation  (136)  becomes 

2 i  9<8in0lV> 


1  d'(rip  )  ¥  1 

r  3r2  r2sin6 


iii*+ 


r2sin20  dp2 

It  is  obvious  that  Equation  (138)  and  (139)  are  only 


k2<|>  -  0 


the  scalar  wave  equations 

V2ip  +  k2\p  ■  0 

expressed  in  spherical  coordinates. 


(136) 


(137) 


(138) 


(139) 


(140) 
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Both  scalar  potentials  and  \p*  are  solutions  of  Equation  (140) 
and  the  complete  solution  of  the  electromagnetic  field  is  obtained  by 
adding  the  partial  fields  derived  from  i|>  and  $*.  Thus 


E  -  +  k2ri{> 

r  3r  2 

(141) 

-  1  32riJ)  imp  Sriji 

0  r  3r36  rsin0  3$ 

(142) 

E  .  1  ifsi  +  m! 

4>  rsin0  3r34>  r  30 

(143) 

H  •  +  ^2,.^* 

r  3r2 

(144) 

..  B  1  32rijj*  (o  +  ime)  3riJ; 

0  r  3r30  rsin0  3<|> 

(145) 

..  1  32riji  (o  +  icoe)  3ri|» 

$  rsin0  3r3<Js  r  30 

(146) 

Thus  the  problem  of  integrating  Maxwell's  equations  has  now  been 
reduced  to  the  problem  of  finding  two  mutually  independent  solutions  of 
the  scalar  wave  equation  evaluated  so  as  to  satisfy  the  boundary 
conditions. 

It  is  easily  verified  that  these  relations  satisfy  Maxwell's  equations. 

Many  solutions  of  the  scalar  wave  equation  in  spherical  coordinates 
are  available  [Refs.  6,  35,  53,  54]  and  the  discussion  presented  here  will  be 
based  partially  on  Stratton's  [Ref.  53]  treatment.  Equation  (138)  is  the 
scalar  wave  equation  in  spherical  coordinates,  and  since  it  is  separable 
in  the  spherical  coordinates  r,  9  and  integral  solutions  exist  of 
the  form 

i|>  “  R(r)  T(9)  ?(♦)  (147) 
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5 

j 

- - - 

,  ---- - 

£ 

> 

\ 

i 

$ 

where 

1 

R(r)  *  a  function  of  r  only 

\ 

* 

i 

T(8)  ■  a  function  of  6  only 

,  i 

P($)  *  a  function  of  only 

Substituting  this  expression  into  Equation  (138) ,  we  obtain 

»  j 

TP  difiaa.  +  “  d<3J1°6^>  +  RI  d2p  +  k2RTP  -  0 

r  dr2  r2sin28  r28in28  d$2 

(148) 

• 

Dividing  by  i|»  *  RTP  and  multiplying  by  r2  gives 

r 

r  d2 (rR)  1  d(sln  ¥>  .,  1  d2P  +  ^  .  Q 

R  dr2  sinST  d@  sin28P  dtp2 

(149) 

' 

Now,  setting  all  terms  that  are  a  function  only  of  r  equal  to  a  constant, 

l 

a2,  we  get 

t 

!■/ 

*. 

;  +  k2r2,  „2 

R  dr2 

(150) 

and 

r  d2.LrM  +  (k2r2  .  a2)R  .  o 
dr2 

(151) 

i 

* 

If  we  let  a2  *  n(n  +1)  n  ■  an  integer 

(152) 

, 

x  -  kr 

(153) 

■j 

then  dx  -  kdr 

(154) 

* 

and  upon  substituting  these  relations  in  Equation  (151)  we  obtain 

::  * 

x2  +  2x  ~  +  [x2  -  n(n  +  1)]  R  ■  0 

dx2  dx 

(155) 

i 

•*l/2 

Now,  if  we  make  the  substitution  R  •  x  y,  Equation  (155)  becomes 

\ 

x2  +  [x2  -  n(n  +  1)  -  t]  y  “  0 

dx2  dx  4 

(156) 

.* 

% 

\ 

| 

i 

* 

\ 
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£ 

or 

x2  +  x  &  +  [x2  -  n(n  +  y)2]  y  =  0  (157) 

dx2  dX  2 

which  is  known  in  the  literature  as  Bessel's  differential  equation  of 
order  n  +  1/2. 


The  general  solution  is  of  the  form 


y  "  UnJn+l/2(x>  +  V„Vl/2<x> 


(158) 


-1/2. 


Now,  since  R  ■  x  y  and  x  ■  kr,  the  general  integral  of  Equation 
(151)  becomes 


/kr 


n+l/2 


(kr)  + 


/kr 


N 


n+l/2 


(kr) 


(159) 


where 


^n+l/2  ^r)  "  ®esse^  function  the  first  kind 

N  . ,  .  (kr)  *  Bessel  function  of  the  second  kind 
n+l/2 

sometimes  called  the  Neumann  function 
We  can  find  an  alternate  solution  in  terms  of  spherical  Bessel 
functions  [Ref.  35]  such  that 


where 


R  =  u  j  (kr)  +  v  n  (kr) 
n  n  n  m 


Vkr> 


Jn+l/2 


(kr) 


(160) 


nn(kr) 


'  2ilr  Nn+l/2(kr) 


Now,  by  substituting  Equation  (150)  into  Equation  (149)  and  upon 
multiplying  the  resultant  expression  by  sin20,  we  obtain 


eln8  J(sln95?> 
T  de 


+  ^  +  a2sin20 

P  d*2 


0 


(161) 
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The  relationship  expressed  by  this  equation  can  be  satisfied  only  if 
the  second  term  is  constant,  i.e. 


1  d2P  2 

where  m2  is  a  constant.  A  single-valued  solution  of  the  form 


(162) 


P  *  a  cos(m$)  +  b  sin(m4>) 
m  m 


(163) 


exists  over  the  range  0  <  <  2JI  for  m  an  integer  such  that  n  >  |m| 

where  a2  ■  n(n  +  1) . 

Substituting  Equation  (162)  into  Equation  (161)  and  multiplying  by 
T/sin20  gives 


sin20 


(164) 


which  is  known  as  the  associated  Legendre  equation.  The  general 


solution  is 


T  =  f  P  m(cos6)  +  g  Q  m(cos0) 
n  n  ®n^n 


(165) 


where  P^  (cos9)  *  associated  Legendre  function  of  the  first  kind 

Qnm(cos0)  -  associated  Legendre  function  of  the  second  kind 

Qnm(cos0)  has  a  singularity  at  the  poles  9  ■  0  and  it;  consequently, 
since  the  fields  exist  at  these  poles,  an  associated  Legendre  function 
of  the  second  kind  does  not  apply  to  our  problem.  Therefore,  g^  must 
equal  zero  and  Equation  (165)  reduces  to 


T  -  f  P  (cos6) 
n  n 


(166) 


which  vanishes  for  m  >  n. 


Therefore,  the  general  solution  of  the  scalar  wave  equation  can  be 


written  as 


"  n 

^(r ,  6,$)  -  l  l  [a  Cos(m<{0  +  b  sin(m<f>)] 
n=0  m=-n  <li 


[Cnjn(kr)  +  dnNn(kr)]  [P^cosS)] 


(167) 


where  a  ,  b  ,  c  and  d  are  arbitrary  constants, 
m  m  n  n 


We  are  now  prepared  to  derive  Mie's  formulae  assuming  the  incident 
radiation  is  linearly  polarized  with  the  x-axis  in  the  plane  of  the 
electric  field  vector  and  the  direction  of  propagation  along  the  positive 
z  axis.  The  origin  is  taken  at  the  center  of  the  sphere  of  radius  a. 

We  designate  the  medium  surrounding  the  sphere  as  medium  1  defined  by 
0],  £},  pi  and  kj;  and  the  sphere  is  designated  as  medium  2  defined  by 
02 ,  £2 >  M2  and  k2 .  Thus,  the  incident  linearly  polarized  radiation  is 
defined  as 

E  -  E  e1(a)t  +  kiz)  (168) 


E  ei(mt  +  kjz) 

0 

(168) 

H  ei(wt  +  k!z) 

0 

(169) 

E  -  H  -  H  -  0 

Z  X  z 

(170) 

We  can  immediately  calculate  the  relationship  between  E  and  H 

x  y 


from  Maxwell's  equations,  i.e. 


3H 

if 


(171) 


Thus, 


For  convenience,  let 


ikl 

Hy  iwui  Ex 

o«o  +  luc 
g  ■  —  iu>W 


(172) 

(173) 

(174) 
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where 


k2  •  o0  *  (i)2yo  -  itupa 


(175) 


Then  may  be  expressed  as 


H  ■  “2 —  E  e 
y  01  o 


i(uit  +  k^z) 


(176) 


Now,  It  is  necessary  to  express  E  and  H  in  terms  of  spherical 

x  y 


coordinates.  To  transform  any  vector 


A'Aa  +  A  a  +  A  a 
xx  y  y  z  z 


(177) 


in  the  Cartesian  coordinate  system  to  the  spherical  coordinate  system. 


A  ■  A  a  +  A.aa  +  A.a. 
r  r  6  0  d> 


(178) 


we  represent  the  position  vector  of  any  point  as 


p  ■  xa  +  ya  +  za 
x  *  y  z 


(179) 


substituting  Equations  (95)  through  (97)  into  Equation  (179)  gives 


p  ■  rsin8cos$a  +  rsin0sin$a  +  rcosBa 
x  y  z 


(180) 


In  the  spherical  coordinate  system  the  unit  vectors  are  defined  as 


a  *  — r —  ■  sin0cos$a  +  sin6sin$a  +  cosda 
r  1 3p  |  x  y  z 


(181) 


a.  »  — r —  •  coa0coa$a  +  cos0sin+a  -  sin0a 

9  iHi  *  y 


(182) 


•*34  •+  •+ 

a.  "  "TT“  ■  ’  •!*»♦«  ♦  coa^a 

♦  ||L|  x  y 


(183) 


Then  solving  these  equations  simultaneously  for  the  unit  vectors  in  the 
Cartesian  systeai  gives 


I 


a  »  sin0cos#a  +  cos0cos#aQ  -  sin^a^ 
X  r  o  <t 


(184) 


a  *  sin0sin#a  +  cos0sin#an  +  cos#a, 
y  r  U  9 


a  ■  cos0a  -  sin0aA 
z  r  vJ 


(185) 

(186) 


Substituting  Equations  (184)  through  (186)  Into  Equation  (177)  we  see 
that  any  vector  A  expressed  in  the  Cartesian  coordinate  system  can  be 
transformed  to  the  spherical  coordinate  system  according  to 

X  •  (Axsin0cos#  +  Ay8in0sin<t>  +  Azcos0)  af 

+  (A^cosOcos#  +  AyCos0sin4>  -  Azsin0)  aQ 


+  (-A  slnq>  +  A  cos#)  a. 
x  y  # 


(187) 


Thus,  the  components  of  the  incident  field  can  be  expressed  in  the 
spherical  coordinate  system  as 


.  E  +  k,rco»9) 

r  o 

(188) 

£  -  E  ei("e  +  k'rcos6,co,9co8* 

0  0 

(189) 

_  i(wt  +  kjrcosO)  .  . 

E  •  -E  e  1  sin# 

9  o 

(190) 

u  _  **1  i(ut  +  k;rcos8)  .  .  .  , 

H  »  E  -s —  e  *  siadsin# 

r  o  Bi 

(191) 

H  .  E  .1(ot  +  k>rco,9,M.9.1o» 

0  OD] 

(192) 

H  ■  E  -jr—  «1<“t  +  k't”,9)cc» 

#  o  8j 

(193) 

Since  these  equations  represent  the  total  field,  it  is  necessary 
to  use  Equation  (141)  through  (146)  to  evaluate  the  scalar  potentials 
#  and  for  the  incident  wave.  For  convenience  in  computation  we  choose 
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to  work  with  Equation  (141),  i.e. 

E  ■  +  k2ri|> 

r  3r2 

Substituting  for  E^  from  Equation  (188)  gives 

<  .1(“t  +  klrco,e).ln9co.*  -  list  +  k*r„ 
0  Sr*  ' 

Now,  applying  Bauer's  formula  [Ref.  55] 


(141) 


(194) 


ircos6 


l  in(2n  +  1)  J  (r)  P  (cos8) 
n  n 

n*o 


(195) 


we  obtain 


eik1rcos6  _  j  in(2n  +  1)  Jn(k!r)  Pn(cos0) 


(196) 


n»o 


If  we  differentiate  this  expression  with  respect  to  6,  i.e. 


elk,rco.e  .  |_  j  tn(2B  +  „  j_(kir)  p_(c„.e, 


36 


(197) 


n»o 


we  obtain 


eikirco808lne  l  in(2n  +  1)  J^r)  Pn(1)(cose)  (198) 

1  n»l 

where  the  associated  Legendre  functions  of  the  first  kind  are  defined  in 
terms  of  the  Legendre  functions  of  the  first  kind  by 

d*P  (cos8) 


P  *(cos8)  ■  sin®8  - — 

n 


(199) 


d(cos8) 


Thus,  it  can  be  readily  verified  that 


dP  (cos8) 

“35 - "n<0<"-9> 


(200) 


Now, Er  can  be  expressed  as 


N 
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Er  .  */<•'  +  klrcos9,sln9cos* 


E/"'  l  ln'‘(2n  +  1)  jn(k1r)  ?n(l)(cos8) 


(201) 


Assuming  a  similar  form  for  \ l>,  i.e., 


*  ■  V1"'  Kn  Vk>r)  '.“’t"*") 


(202) 


we  obtain  upon  substituting  it  and  Equation  (201)  into  Equation  (194) 


the  relation 

32rj  (lqr)  ,11-1,.  .  .* 

r - - - +  [k?r2  -  i j  (kjr) 

3r2  1  n  n 

which  is  nothing  but  Bessel's  equation  where 

„  _  iH*l (2n  +  1) 
n  1  n(n  +  1) 


(203) 


(204) 


Therefore,  for  the  incident  wave  the  scalar  electric  potential  is  given 


j.  iwt 


and  similarly,  for  the  magnetic  potential 


E  eiwtsin$ 


*  t  e  sinf  *  m  n  m 

*i  -  -Hr —  ‘  er  JH)  J„<k>r>  p„  ’  <«••> 

To  determine  the  scalar  field  potentials  for  the  scattered  wave  and 
for  the  wave  within  the  sphere,  it  is  necessary  to  determine  the  arbitrary 
constants  of  Equation  (167)  in  such  a  way  as  to  satisfy  the  boundary 
conditions  and  the  conditions  at  the  origin  (center  of  the  sphere)  and 
at  infinity.  The  conditions  to  be  satisfied,  then,  are 


a.  The  field  within  the  sphere  must  be  finite  at  the  origin. 

Since  the  spherical  Bessel  function  of  the  second  kind,  n  (kr) ,  has 

n  ’ 

a  singularity  at  r  -  0  and  becomes  infinite  there,  d  »  0  for  the  region 

n 

within  the  sphere. 

b.  The  scattered  field  must  vanish  at  infinity.  The  function 

nn<2V>  -  Jn(kr)  -  lnn<kr)  -  JgL  h£>/2<1*)  (207) 

is  known  as  the  Hankel  function  of  type  2  (or  Bessel  function  of  the  third 

kind)  and  has  the  property  that  it  vanishes  at  infinity.  Therefore, 

c  -  1  and  d  -  -i  in  the  solution  of  the  scattered  wave, 
n  n 

c.  The  required  boundary  conditions  are  that  the  tangential  com¬ 
ponents  of  the  electric  field  vector  and  the  magnetic  field  vector  must 
be  continuous  across  the  surface  of  the  sphere,  i.e. ,  at  r  *  a 


Eei 

m  E 

a02 

(208) 

V 

E(J>2 

(209) 

He, 

“H02 

(210) 

“  %2 

(211) 

If  we  assume  that  the  form  used  to  represent  the  incident  wave 
dictates  the  form  of  the  complete  solution,  then  the  scattered  wave  and 
the  field  within  the  sphere  can  be  represented  by  the  following  scalar 
potentials: 

a ;  Scattered  wave  outside  of  sphere 
„  imt 

E  e  COS41  00  /-x  \ 

*  -  -2— ; -  7  a  h(2;(k,r)PU;(cos0)  (212) 

rs  k’  L.  n  n  1  n 

1  n*l 
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(213) 


_  i<»!t  .  . 

E  e  sin4>  00  ikj 


*  “o  r  (2)  (i) 

K - in -  b„TThn  (klr)P„  (cos9) 

n-i 


b.  Wave  within  the  sphere 


lu)t 

E  e  cos4> 
o 


l  wn  j  (k2r)P^U(cos0) 
n-1 


_  i<*)t  .  . 

*  E  e  sinij»  *  ik2  v 

*  «-£— jJ -  l  x  —  j  (k2r)Pr;(cos0) 

k2  n.i  n  s2  n  n 


(214) 


(215) 


Now  since  we  know  that  the  total  field  outside  the  sphere  is  equal 
to  the  sum  of  the  Incident  field  and  the  scattered  field,  we  apply  the 
boundary  conditions  to  Equations  (142),  (143),  (145)  and  (146)  to  obtain 


8r(ij»  +  *  )  Sr* 
is  o 


3r(*  +  **)  3ri|>* 

i  a  _  o 


(216) 


(217) 


c»i  £r <<*»i  +  *8)1  -  02(r^0) 


(218) 


8j  (r(vi  *8)1  »  B2(r*o) 


(219) 


which  are  valid  only  at  r  ■  a.  Upon  substituting  Equations  (205)  and 
(206)  and  Equations  (212)  through  (215)  into  these  boundary  conditions  where 
the  derivatives  of  Bessel's  and  Henkel's  functions  are  denoted  by  primes: 

‘“"'mtH}  rr  *;;«'•>  ♦  *«  iK’r<k**>  •  -» fe  <220) 


l"  h  J«<ki*)  *  b.  «r  l^)<k>*>  ■  *.  r2  ua) 
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“2 


in’‘  n(n  7 1)  kf  Vkla)  +  an  kf  in<k2a)  (222) 

*"  Vk>'>  +  b„ihi2)'kla>  -  W1**’  <*”> 

Upon  eliminating  w  from  the  first  and  third  equations  and  x  from  the 
n  n 

remaining  pair  we  obtain  the  following  expressions  for  a  and  b  : 

n  n 


Ln-i  2n+l  aiVkia)jn(k2a)  "  °2 Jn(k2a) j^C^a) 

n(n  +  X)  ^J^aJh^kja)  -  ^'(k^h^^a) 


(224) 


b  xn-l  2n  4-  1  eiVM^ayO  -  B2Jn(k2a)Ji;(k1a)  ^ 
n  n(n  +  D  e2jn(k2a)h<2) (kja)  -  Bj j;(k2a)h<2) (kja) 


Similarly, 

.  <"-1  2„  +  t  -  k2„1Jr[Q11a)h<2>;ic1a) . 

Wn  n(n  +  1)  ....  *,  (2),.  s  v  .  ,,  \u(2)/.  .  “  * 

k1a1jn(k2a)hn  (k^)  -  k2a2 jn(k2a)hn  (kia) 

2n  +1  -  S.VMJhf'Jk.a) 

n  n(n  "  D  3i j'(k2a)h£2) (kl3)  -  82jn(k2a)h<2) (kl3) 

We  have  now  determined  the  unknown  constants  of  Mie's  problem  and 
we  can  finally  calculate  the  field  at  any  point  inside  or  outside  the 
sphere.  For  our  problem  we  are  only  interested  in  the  behavior  of  the 
scattered  wave  and  on  substituting  Equation  (212)  and  (213)  into  Equations 
(141)  through  (146)  we  can  derive  expressions  for  the  components  of  the 
field  vectors  of  the  scattered  wave.  In  deriving  the  radial  components, 

Er  and  H^,  it  is  necessary  to  evaluate  the  first  and  second  derivatives 
of  Equation  (207)  and  while  it  is  a  tedious  job,  the  use  of  the  following 
recurrence  formulas  readily  facilitate  the  calculations: 
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dJ  (y) 

V  *"  V 


X"" 

dy 

-  -  J  (y)  - 
y  x  * 

J,+, (y) 

(228) 

dJ  (y) 

A 

dy 

*  —  J  ( y ) 

y  x  ' 

+  J  (y) 

X-l  J 

(229) 

where  (y)  *  Bessel  function  of  the  first,  second  or  third  kind  of  order 
x  (may  be  fractional,  integral  or  complex)  and  argument  y 
Thus,  we  obtain  the  following  expressions  for  the  radial  components 
of  the  field: 


_  lust  , 

E  e  cos<4>  " 
o 


E  = -  l  n(n  +  l)a  h(2)(k,r)P(l)(cose)  (230) 

r  2  n  n  A  n 


kirz  n=l 


H 


ltot  .  . 

iEoe  sin<(>  ™ 


l  n(n  +  l)b  h(2)(k1r)P(1)(cos0)  (231) 

2  n  n  n 


Bir^  n=l 

In  deriving  the  tangential  components  of  the  field  it  is  convenient 
to  make  the  substitution 

/irkir 


This  gives 


„  ioit  , 
E  e  cos<J) 
o 


c  tusif/  —  ,  -  \  \  » 

E  ■ - — -  [  J  a  (kir)P'  '(cos0)sin0 

0  kir  n  n  1  n 


(232) 


-  ib  z(2)(kir)P(l)(cos0) 

n  n  1  7  n  v  sin0J 


(233) 


o  iwt  .  , 

E  e  sin$  «  , 

E  -  -  - -  [  l  a  z(2)(kir)P(  }(cos0)  — 

<p  k*r  n  n  1  n  sin 


sine 


-  ib  z^  (kjr)P^1^  (cos0)sin0] 


n  n 


(234) 
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H  -  - 
8 


_  i(i)t  .  . 

E  e  sin$  00 
o 


^  a  ^ 


sin0 


+  ib  z^  (kjr)P^  (cos6)sin6] 
n  n  n 


(235) 


r,  v.  , 

E  e  cos$  ®  .  .  . 

H$  "  - 1  ^  Vn  J(klr)Pn  (cosQ)  sine 

n*l 

+  iVf)'k‘r>Pn1>(“s9)'iI561  <236) 

The  primes  denote  the  derivative  of  a  function  with  respect  to  its 
argument.  We  have  now  completed  the  solution  of  Mie's  problem.  While 
the  solution  was  derived  for  a  single  sphere,  the  results  apply  to  any 
optically  thin  medium  containing  many  spherical  particles  of  the  same 
composition  and  of  the  same  or  different  ulzes  if  the  distance  between 
particles  is  large  compared  to  the  wavelength  of  the  incident  radiation. 
We  can  then  assume  that  there  is  no  multiple  scattering  and,  therefore, 
the  total  intensity  scattered  is  equal  to  the  sum  of  the  intensities 
scattered  by  each  particle. 

In  summary,  we  have  derived  the  formulas  for  the  electric  and 
magnetic  field  components  which  describe  the  scattered  field  arising 
from  a  linearly  polarized  monochromatic  plane  wave  incident  upon  a 
spherical  particle  of  arbitrary  composition  and  size.  We  see  from 
Equations  (230)  and  (231)  that  in  the  far  field  (i.e.,  at  distances 
very  much  greater  than  the  wavelength  of  the  incident  radiation)  the 
radial  components  of  the  scattered  wave  vanish  as  the  inverse  square 
of  the  distance  from  the  origin,  while  the  tangential  components 
diminish  as  the  inverse  of  the  distance.  Thus  in  the  far  field  the 
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radial  components  may  be  neglected  In  comparison  with  the  tangential 
components  and  the  fields  can,  therefore,  be  assumed  to  be  transverse 
to  the  direction  of  propagation.  Furthermore,  the  scattered  radiation 
from  the  sphere  is,  in  general,  elliptically  polarized  except  in  the 
planes  of  scattering  defined  by  <f>  ■  0  and  <j>  ■  -j,  i.e.  the  xz  and  the  yz 
planes.  We  see  from  Equations  (233)  and  (234)  that  E^  vanishes  for 
<j>  *  0  and,  alternately,  E0  =  0  for  4>  *  4;  consequently,  the  scattered 
radiation  is  linearly  polarized  for  all  scattering  angles,  0,  when  the 
plane  of  scattering  is  parallel  or  perpendicular  to  the  electric  field 
vector.  This  polarization  characteristic  provides  the  basis  for  the  cloud/ 
target  discrimination  technique  described  in  the  next  chapter. 

3.3  APPLICATION  OF  MIE  THEORY  TO  ATMOSPHERIC  AEROSOLS 

Many  investigators  have  considered  the  consequences  of  Mie's  theory 
in  terms  of  the  physical  parameters,  such  as  the  wavelength  of  the  incident 
radiation,  the  index  of  refraction  of  the  sphere  and  the  surrounding 
medium,  the  composition,  density,  size  and  size  distribution  of  the 
scattering  particles,  etc.  In  addition,  there  Is  an  immense  amount  of 
literature  on  the  application  of  the  theory  to  other  fields  of  research. 

We,  therefore,  refer  the  reader  interested  in  a  complete  coverage  of 
the  theory  and  all  of  its  aspects  to  the  publications  already  referenced. 
For  purposes  of  our  research  we  are  primarily  interested  in  the  extinc¬ 
tion  coefficient,  o,  and  the  angular  volume  scattering  coefficient,  8(0) 
defined  as  the  ratio  of  the  flux  scattered  into  a  unit  solid  angle  by  a 
unit  volume  of  a  homogeneous  medium  containing  particles  of  different 
radii  in  the  direction  0,  to  the  flux  incident  on  the  unit  cross  section 
of  this  volume. 
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Both  o  and  8(6)  are  defined  mathematically  for  the  far  field  where 

(2) 

the  spherical  Bessel  function  h^  (kr)  converges  very  rapidly  for  large 
values  of  its  argument  and  can,  therefore,  be  approximated  by 

,n+l 


Substituting  this  expression  into  Equations  (212)  and  (213)  gives 


(237) 


A  s 

S 


i(o)t-kir) 
iE  e  1  cos$ 
o _ _ 

*  ' 


l  in+1a„pf ^ (cos6) 
n**l 


n  n 


(238) 


l<u)t-kir)  ^ 

*  _ 0 _ _ 

=  Mir  nml 


00 

l  i^b^^cose)  (239) 


and  accordingly  the  tangential  field  components  become 
iE  e1(“t_klr) 


cos4> 


Ee  = 


kjr 


{  l  i1^1  [a  (cos0)sin9  +  bP^(cos0) 

1  ",  n  n  n  n  si.no 

n*l 

(240) 


E$  * 


iE  ei(a)t“kir)sin<(> 
o _ _ 

kjr 


CO 

{[  i1^1  [anPn1>(CO80>  rd*  +  bnP^l)(cos0)sin6]} 


n-1 


sin6  n  n 


(241) 


He  *  ‘ 


E  e1(a)t-klr)sin* 
0 


FT -  {  I  i”*1  U  P(l)(cos0)  ~~  +  b  P(l) <cos0)sin0] } 

8ir  1  1  n  n  '  sin8  n  n  1 

n  x  /a/ a\ 


(242) 


„  i(wt-kir)  , 
E  e  1  co8$ 

o 


H? - -  f  l,  l0+‘  tanP‘l)t«»6)»i«e  +  Vi"<co,e)  ihe’l 

1  n*l 

(243) 

Now,  Equations  (240)  and  (243)  can  be  rewritten  as 

U  el(»t-k,r)coM 


,(i) 


Ee  " 


{A„<e)} 


(244) 
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iE  ei(tot“klr) 


cos$ 


kir 


{Ax<e>} 


(245) 


He  =  ' 


_  i(ii)t-kir)  .  . 
E^e  1  sin<j> 

eir 


{ajl(6)  } 


(246) 


E  ei(aJt-klr)8in^ 

H*  =  - - 317 -  fAii<e» 


(247) 


where  A()(0)  and  AA(0)  are  called  amplitude  functions  and  are  equal  to 


A||(6)  =  l  i 0+1  [b  E  (cos0)  -  a  c  (cos0)J  (248) 

"  ,  n  n  n  n 

n-1 


A, (6)  -  l  i1*1  fa  E  (cos0)  -  be  (cos0)] 
,  n  n  n  n 

n-1 


where  we  let 


dP(1)(cos0)  dP(1)(coa0) 

^n(cos9)  "  -2"d0 - Sin6  d(cos0)  ■ 


(249) 


(250) 


Sn(c°89)  -  p‘‘>(cose)  ^ 


(251) 


We  introduce  the  intensity  functions,  i(|  and  iA,  defined  by  Mie  to  be 

i„  "  l A,,(e)  |2  (252) 

iA  -  |Aa(0)|2  (253) 


where  il(  is  proportional  to  the  angular  component  of  the  scattered 
field  whose  electric  vector  is  in  the  plane  of  scattering  and  iA  is 
proportional  to  the  component  whose  electric  vector  is  perpendicular  to 
the  plane  of  scattering. 

We  now  state  without  proof  the  expression  for  the  angular  volume 
scattering  coefficient  8(0)  for  a  medium  containing  spherical  particles 
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of  different  radii 


6(6)  -  —  l  N(a)Aa(ij,  +  iA)  (254) 

8ir2  a 

and  the  extinction  coefficient  assuming  negligible  absorption,  i.e., 

0  ■  a,  is  obtained  by  integrating  this  expression  over  all  directions; 
thus 

2  IT  IT 

o  -  J  /  6(8)sin0d6d$  (255) 

o  o 


7T 

-  2ir  /  6(0)sin8d0  (256) 
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SECTION  IV 


THE  DERIVATION  OF  THE  CLOUD/TARGET 
POLARIZATION  DISCRIMINATION  TECHNIQUE 

4.1  QUALITATIVE  DESCRIPTION  OF  THE  DISCRIMINATION  TECHNIQUE 

The  general  solution  for  the  scattered  fields  (Includes  refracted, 
reflected,  and  diffracted  fields)  arising  from  an  electromagnetic  wave 
Incident  on  a  smooth  (plane)  Interface  between  two  media  Is  well-known 
and  in  the  preceeding  chapter  we  discussed  an  exact  solution  for  one 
such  problem.  On  the  other  hand,  neither  a  general  nor  an  exact  solu¬ 
tion  has, been  determined  as  yet  for  similar  problems  involving  a  rough 
interface  between  media  where  "rough  Interface"  is  defined  as  a  surface 
that  will  diffusely  scatter  the  radiation  of  an  incident  plane  wave. 
Correspondingly,  then,  a  smooth  interface  is  defined  as  a  surface  that 
will  reflect  the  incident  radiation  specularly.  The  degree  of  roughness 
(or  smoothness)  is  dependent  on  the  wavelength  and  angle  of  incidence 
of  the  incident  radiation,  and  can  be  described  in  terms  of  a  widely  ac¬ 
cepted  roughness  criterion  attributed  to  Lord  Rayleigh.  Briefly 
stated,  Rayleigh's  roughness  criterion  arbitrarily  considers  a  surface 
to  be  smooth  if  either  the  ratio 

X  -  0  (257) 

or  the  angle  of  incidence,  y,  measured  from  the  normal  to  the  surface 
approaches  0,  i.e., 

Y  -  0  (258) 

where  *  is  the  wavelength  of  the  incident  radiation  and  t  is  the  height 
of  the  surface  irregularities  (See  Fig.  2). 
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While  the  general  problem  of  describing  the  amplitude  and  phase 
relationships  of  the  field  components  reflected  from  rough  surfaces  has 
received  considerable  attention  from  several  writers,  very  little  effort 
has  been  devoted  to  determining  whether  or  not  the  polarization  of  the 
scattered  wave  is  the  same  as  that  of  the  incident  wave.  A  thorough 
treatment  of  rough- surface  scattering  is  beyond  the  scope  of  this  paper 
and  the  reader  is  referred  to  Beckmann  and  Spizzichino' s  book  [Ref.  4] 
for  a  clear  and  comprehensive  discussion.  They  have  shown  that,  in  gen¬ 
eral,  the  scattered  wave  arising  from  an  arbitrarly  polarized  plane  wave, 
incident  on  an  arbitrarily  inclined  rough  surface  of  arbitrary  conduc¬ 
tivity  whose  radii  of  curvature  is  large  compared  to  the  wavelength,  will 
be  depolarized. 

We  have  now  discussed  the  resultant  phenomenon  arising  from  the 
interaction  of  a  linearly  polarized,  monochromatic  electromagnetic 
plane  wave  incident  on  (1)  a  homogeneous  medium  containing  a  size  dis¬ 
tribution  of  spherical  particles  the  radii  of  which  are  on  the  same  order 
of  magnitude  as  the  wavelength  of  the  Incident  radiation,  and  (2)  an 
arbitrarily  Inclined  rough  surface  of  arbitrary  conductivity  as  de¬ 
fined  by  Rayleigh's  roughness  criterion.  For  the  first  case  the  Hie 
theory  of  scattering  shows  that  in  the  far  field  (r»X)  the  scattered 
radiation  in  planes  either  parallel  or  perpendicular  to  the  electric 
field  vector  of  the  incident  wave  is  linearly  polarised  for  all  scatter¬ 
ing  angles.  And  for  the  second  case  Beckmann  and  Spizzichino  have 
shown  that,  in  general,  the  scattered  radiation  will  be  depolarized  if 
the  radius  of  curvature  of  the  surface  variations  is  comparable  co  the 
wavelength  of  the  incident  radiation. 
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As  we  discussed  in  Seccion  I  the  primary  limitation  to  the  use  of 
the  laser  in  many  atmospheric  applications  is  the  severe  path  loss  due  to 
absorption  by  atmospheric  gases  and/or  scattering  by  cloud  and  fog  par¬ 
ticles  «  Furthermore,  in  detection  and  ranging  applications  the  cloud 
and  fog  particles  may  scatter  the  transmitted  radiation  toward  the  re¬ 
ceiver,  thus  producing  a  false  return  which  must  be  distinguished  from 
the  target  return.  It  is  this  latter  situation  that  we  are  most  con¬ 
cerned  with,  i.e.,  the  problem  of  distinguishing  between  a  cloud  reuurn 
and  a  target  return.  Since  the  majority  of  useful  lasers  emit  radiation 
in  the  visible  or  near  infrared  portion  of  the  electromagnetic  spectrum, 
atmospheric  clouds  or  fogs  actually  represent  Mie  scattering  media  and 
most  targets — earth  terrain  including  man-made  irregularities,  , airborne 
objects  such  as  aircraft,  etc., —  can  be  considered  to  have  rough  sur¬ 
faces  at  these  wavelengths.  It  is  obvious,  then,  that  one  possible 
method  of  distinguishing  a  cloud  return  from  a  target  return  could  be 
based  on  the  difference  in  the  polarization  characteristics  of  each 
return. 

A  qualitative  description  of  a  possible  cloud/target  polarization 
discrimination  technique  is  illustrated  in  Fig.  3.  The  output  of  the 
laser  is  plane  polarized  with  the  plane  of  vibration  of  its  electric 
vector  (represented  by  the  dots)  perpendicular  to  the  page  and  a  direc¬ 
tion  of  propagation  as  illustrated.  A  polarizer  is  placed  at  the  re¬ 
ceiver  input  that  will  transmit  only  those  components  whose  electric 
vectors  (represented  by  the  arrows)  vibrate  at  right  angles  to  the 
direction  of  the  electric  vector  of  the  transmitted  wave  and  will  elimi¬ 
nate  those  that  vibrate  parallel  to  this  direction.  When  used  in  this 
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Qualitative  Illustration  of  the  Cloud/Target  Polarization  Discrimination  Technique 


manner  the  polarizer  is  usually  called  an  analyzer.  We  see  then  that, 
according  to  the  Mie  theory  of  scattering,  those  portions  of  the  trans¬ 
mitted  wave  that  are  scattered  by  the  spherical  cloud  particles  are  of 
the  same  polarization  as  the  transmitted  wave.  The  radiation  re¬ 
flected  from  a  diffuse  target,  however,  is  depolarized  (represented  as 
a  combination  of  a  dot  and  an  arrow)  according  to  the  rough  surface 
theory  of  scattering.  We  assume  that  any  depolarization  effects  due  to 
multiple  scattering  are  negligible.  Therefore,  the  receiver  will  "see" 
only  those  components  of  the  radiation  reflected  from  the  target  whose 
electric  vectors  vibrate  at  right  angles  to  the  electric  vector  of  the 
transmitted  wave  while  all  returns  due  to  scattering  from  cloud  parti¬ 
cles  will  be  eliminated  by  the  analyzer.  We  use  the  term  "eliminated" 
to  arbitrarily  represent  the  attenuation  of  optical  energy  by  any  one 
of  the  processes  of  absorption,  reflection,  refraction,  or  scattering. 

4.2  MATHEMATICAL  DERIVATION  OF  THE  DISCRIMINATION  TECHNIQUE 

The  purpose  of  this  section  is  to  derive  a  mathematical  expression 
analogous  to  the  radar  range  equation  for  a  ladar  system.  In  addition, 
we  will  derive  an  expression  for  the  radiation  scattered  by  a  cloud  in 
the  direction  of  the  ladar  receiver.  Finally,  we  will  evaluate  the  ef¬ 
fectiveness  of  the  discrimination  technique  in  terms  of  the  signal-to- 
noise  (target  return-to-cloud  return)  ratio.  We  chose  the  signal-to- 
noise  ratios  as  our  criterion  because  it  is  the  basic  system  parameter 
which  determines  the  minimum  detectable  signal  necessary  to  achieve  a 
specified  probability  of  detection  without  exceeding  a  specified  prob¬ 
ability  of  false  alarm.  Strictly  speaking,  the  return  from  a  cloud 
represents  a  range  ambiguity  rather  than  a  noise  factor  limiting  receiver 
sensitivity;  however,  if  we  define  "noise"  as  unwanted  electromagnetic 
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energy  which  interferes  with  the  ability  of  the  receiver  to  detect  and 
recognize  the  target  return,  then  we  can  consider  the  return  from  a 
cloud  as  noise. 

We  know  that  the  energy  of  electromagnetic  radiation  (i.e„,  radiant 
energy)  emitted  or  reflected  per  unit  time  by  an  object  is  commonly  re¬ 
ferred  to  as  radiant  flux  or  radiant  power.  The  physical  quantity  ra¬ 
diant  flux  is  considered  as  a  fundamental  concept  £rum  which  other  phy¬ 
sical  concepts  called  ''radiometric  quantities"  have  been  derived.  An 
analogous  set  of  concepts  referred  to  as  "photometric  quantities"  has 
been  derived  for  the  visible  part  of  the  spectrum.  This  system  is  psy¬ 
chophysical  in  nature  being  based  on  the  effect  of  light  on  the  eye. 

Since  the  radiometric  quantities  refer  to  radiation  in  general,  we  will 
use  this  system  in  our  development.  An  excellent  treatise  on  both  radio- 
metric  and  photometric  concepts  is  presented  by  Middleton  [Ref.  41].  The 
definitions  of  the  radiometric  quantities  we  will  use  are  given  in 
Table  2.  Note  that  the  irradiance,  E,  has  the  dimensions  of  power  den¬ 
sity  and  for  a  point  source  alternate  expressions  can  be  derived  that 


are  valid  in  the  far  field. 
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4.2.1  DERIVATION  OF  THE  LADAR  RANGE  EQUATION 

If  the  linear  dimensions  of  the  source  are  small  compared  to  the 
distance  between  the  transmitter  and  the  target,  the  energy  may  be  assumed 
to  be  radiating  from  a  point  source.  Furthermore,  assuming  that  the 
target  is  far  enough  away  so  that  the  transmitted  rays  can  be  considered 
parallel  to  the  line  of  incidence  between  the  transmitter  and  the  target, 
the  radiant  intensity,  I,  at  the  target  is 

I  -  Ioe“0r  (263) 

Let  ds  in  Fig.  4  represent  an  infinitesimal  area  of  the  target  surface 
whose  normal  makes  an  angle  y  with  the  line  segment  r  to  the  transmitter 
and  an  angle  6  with  the  line  segment  l  ■  v^r2  +  d2  to  the  receiver.  The 
solid  angle,  dm,  subtended  by  ds  at  the  source  is 

4w  -  (264) 

r2 

From  Table  2,  the  radiant  intensity  of  the  source  can  be  written  as 

I  -  £  (265) 

Therefore,  the  radiant  flux  incident  on  the  area  ds  is 

dP  =  Idw  =  -C°"Y-d-S  (266) 

r2 

and  on  substituting  Equation  (263)  into  Equation  (266)  we  obtain 

I  e  °rcosyds 

dP  =  -2 -  (267) 

r2 

The  target  surface  has  a  reflective  efficiency  represented  by  its  re¬ 
flection  coefficient  p.  The  radiant  flux,  dP^,  reflected  from  the  sur- 
surface  element  ds  is 
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(268) 


dPfc  *  PdP  * 


t  -or  , 
pi  e  cosyds 
o 

r2 


Assuming  the  target  surface  to  be  a  diffuse  scatterer,  Lambert's  cosine 
law  states  that  the  total  radiant  flux  reflected  from  a  surface  element 
into  a  hemisphere  is 


dPt  »  irBds 


Thus, 


B  « 


t  -or 
pi  e  cosy 
o 


irr* 


(269) 


(270) 


In  addition,  according  to  Lambert's  law  the  radiant  intensity,  dl^,  re¬ 
flected  by  the  surface  element  ds  in  the  direction  of  the  receiver  is 
defined  as 


dl^  «  Bcos6ds 


and  on  combining  Equations  (270)  and  (271)  we  obtain 


-or  r 

pi  e  cosycos6 

dl  - - 2  -  ....  dg 

trr2 


(271) 


(272) 


Assuming  the  infinitesimal  surface  area  ds  approximates  a  point  source, 

the  irradiance,  dE  ,  of  ds  incident  on  the  receiving  aperture  is, 

K 


according  to  Equation  (262), 


dl 


dE. 


R 


R  *2 

where  the  radiant  intensity,  dIR,  at  the  receiver  is  given  by 

-a£ 


(273) 


dlR  -  dlte 


(274) 
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Thus, 


Pi  d  °^r+^c08Ycos<5 


ds 


(275) 


rr2*2 

and  the  receiver  power  (radiant  flux)  based  on  the  infinitesimal  target 
surface  area  is 


dP, 


R 


>IQARe  0 (^^cosY Cos'S 
7rr2£2 


ds 


(276) 


(277) 


or 


pi  A_e"a  (r^>cosYCOsd 

«»_ .  -O - d8 

R  rr2  (i^+d2  ) 


(278) 


where  is  the  area  of  the  receiving  aperture. 

Equation  (278)  is  but  one  form  of  the  general  expression  for  the 
ladar  range  equation  and  if  the  separation,  d,  between  the  transmitter 
and  receiver  is  very  small  compared  to  the  distance  to  the  target, 
Equation  (278)  reduces  to 


pi  e  ^0rcos2YA_ 
dP_  =  —2 - 2  d8 


R 


irr’ 


or 


Pi  e”2orco82YAp 

P  a  - - - 3  /  d8 


(279) 


(280) 


irr’ 


Since,  in  general,  the  output  power,  P^,  and  the  beamwidth,  8^,  of 


a  laser  source  is  known,  the  radiant  intensity,  IQ,  is 

4P„ 


T 


(281) 
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Thus,  we  see  that  the  fourth  power  range  relationship  of  radar  systems 
is  also  valid  for  ladar  systems.  Furthermore,  in  cases  where  the  target 
is  larger  than  the  transmitter  beam,  the  target  return  varies  as  the 
range  squared. 

4.2.2  DERIVATION  OF  THE  CLOUD  RETURN  EQUATION 
In  deriving  an  expression  for  the  radiant  flux  at  the  receiver 
(which  is  now  assumed  to  be  located  at  the  transmitter)  due  to  scatter¬ 
ing  from  the  cloud,  we  consider  an  infinitesimal  element  of  volume  of  the 
cloud  (Fig.  S)  where 

dv  ■  dAdx  (286) 


■  dux2dx 


(287) 


89 


Figure  5.  Geometry  for  Derivation  of  Cloud  Return 
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The  irradiance,  E,  of  the  source  incident  on  this  infinitesimal 


volume  of  cloud  particles  is 


E  - 


1  e 
o 


-ox 


(288) 


If  we  can  assume  the  radiation  scattered  from  the  element  of  volume, 

dv,  as  coming  from  a  point  source  of  which  the  intensity  is  proportional 

to  the  number  of  particles,  then  the  radiant  intensity,  1  ,  for  any 

c 


scattering  angle,  6,  is  given  by 


dl  ■  B(0)Edv 
c 


(289) 


In  the  direction  of  the  receiver,  i.e.,  for  0  ■  tr,  the  radiant  in¬ 
tensity  is 

dl  ■  0(ir)Edv  (290) 

c 


8(ir)I  e  °* 

- - 2 - dv  (291) 

x2 

* 

Accordingly,  then,  the  irradiance,  dER,  at  the  receiver  is 


-  dIR 

dE  -  ~ 

R  X2 

(292) 

where 

,,  -OX 

'  dIc* 

(293) 

Thus, 

.  8(1.)  I.*'2** 

dE_  •  - dv 

R  4 

XH 

(294) 

and  the 

total  power  received  from  the  element  of  volume  is 

K  •  Y‘« 

(295) 

.  a(’>Yo*'J0X  dv 

X4 

(296) 
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On  substituting  Equation  (287)  into  Equation  (296)  we  obtain 


e“2ox 

dpR  -  eoOA^d*  — —  dx 


(297) 


®R  e~2oX 

8(7t)ArIo2it(1  -  cos  -|)  — dx 

X 


(298) 


where  —  is  the  half-angle  of  the  receiver's  field  of  view.  There 
fore,  the  radiant  flux  (power)  at  the  receiver  is 


.  0  r  -2ox 

Pr  -  BO^I^sd  -  cos  -|)  /  - - dx 


(299) 


r  x‘ 
o 


The  complicated  integration  process  suggested  by  Equation  (299)  is 
beyond  the  scope  of  this  paper.  An  excellent  coverage  of  Integrals  of 
this  type  (usually  called  exponential  integrals  depending  on  the  limits 
of  integration)  is  given  in  Appendix  I  of  Kourganoff's  book  [Ref.  34]. 
The  omission  of  the  solution  will  not  affect  the  results  of  our  deri¬ 
vation.  For  convenience,  then,  we  let 


r  -2ox 

E2(2o)  •  /  - dx 


(300) 


r  x* 
o 


Thus,  Equation  (299)  becomes 


PR  -  8(*)ArIo2w(1  -  cos  ~)E2(2o) 


(301) 


4.2.3  SIGNAL -TO- HOI SB  ANALYSIS 


In  the  derivation  of  Equations  (280)  and  (301)  nothing  was  said  about 
the  state  of  polarisation  of  the  transmitted,  reflected  or  scattered 
waves.  We  assumed  the  output  radiation  of  the  laser  to  be 
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unpolarized  where  p  is  the  reflection  coefficient  for  unpolarized  radia¬ 
tion  and  8 (it)  is  the  volume  backscatter  coefficient  for  unpolarized  ra¬ 
diation.  We  begin  our  evaluation  of  the  polarization  discrimination 
technique  by  repeating  Equations  (280)  and  (301) : 

.  f  -2or  5 
pA  I  e  cos^y 

PR - 2 - 7 -  /  ds  (280) 

irr4  s 

PR  -  6(tt)ArIo2tt(1  -  cos  -|)E2(2o)  (301) 

If  we  now  insert  a  linear  polarizer  at  the  transmitter  and  an  analyzer 
at  the  receiver  we  obtain 


.  T  -2or  p 

n2  pl¥oe  cos  Y 

2  4 
irrH 


/  ds 


jr 

2 


V">Vp2it(1 


cos  — |)E2(2o) 


(302) 


(303) 


where  the  subscript  "1"  refers  to  the  fraction  of  the  plane-polarized 
transmitted  wave  which  is  reflected  with  its  polarization  vector  orth¬ 
ogonal  to  that  of  the  transmitted  radiation.  Similarily,  the  subscript 
"  II"  refers  to  radiation  with  its  polarization  vector  parallel  to  that 
of  the  transmitted  wave;  p  is  the  transmission  efficiency  of  the  polar¬ 
izers  (assumed  to  be  the  same  for  both).  The  1/2  represents  the  re¬ 
duction  in  transmitter  output  due  to  the  polarization  process. 

We  now  define  the  signal- to-noise  ratio  for  the  situation  where  the 
output  of  the  transmitter  is  unpolarized  as  the  ratio  of  Equation 


(280)  to  (301),  i.e.. 


(304) 
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where  the  subscript  "+"  denotes  that  the  transmitted  radiation  is  un¬ 
polarized.  Similarly,  the  signal-to-noise  ratio  for  the  system  whose 
transmitter  output  ie  plane  polarized  is 


Upon  careful  scrutiny  of  Equation  (306)  we  observe  that  while  the  value 
of  the  term  in  parenthesis  can  never  be  greater  than  unity  (and  is  usually 
much  less),  the  value  of  the  term  within  the  brackets  can  be  very  large 
if  it  should  turn  out  that  6„(v)  can  be  reduced  to  a  very  small  value 
by  the  use  of  a  perfect  analyzer  at  the  receiver.  We  see,  then, 
that  the  signal-to-noise  ratio  for  the  "cross-polarized"  system  will 
be  much  greater  than  that  of  the  unpolarized  system  where  the  degree  of 


improvement  is  given  by  the  polarization  discrimination  coefficient. 


defined  as 


PDC  *  <w] 


(307) 


We  have  neglected  the  effects  of  other  noise  sources  in  our  devel¬ 
opment  because  they,  in  general,  limit  the  ability  of  the  receiver  to 
detect  the  wanted  signal  rather  than  its  ability  to  discriminate  be¬ 
tween  wanted  and  unwanted  signals.  However,  for  situations  where  the 
effects  of  other  noise  sources  are  not  negligible,  the  appropriate  ex¬ 
pressions  should  be  modified  accordingly. 
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SECTION  V 


THE  LABORATORY  SIMULATION  OF  THE  FORMATION  OF  ATMOSPHERIC  AEROSOLS 

5.1  INTRODUCTION 

The  practical  application  of  the  cloud/target  discrimination  tech¬ 
nique  developed  in  the  preceding  section  is  dependent  on  the  experimental 
verification  of  the  technique  both  in  the  laboratory  and  in  the  field. 

A  high  degree  of  confidence  in  any  conclusion  based  on  the  resultant 
data  from  such  an  effort  can  be  realized  only  if  quantitative  descrip¬ 
tions  can  be  obtained  of  both  the  target  and  the  cloud.  In  order  to 
satisfy  these  "conditions  for  confidence"  in  the  laboratory  it  is  nec¬ 
essary  to  create  an  artificial  cloud  that  best  represents  the  atmos¬ 
pheric  variety.  Even  though  a  laboratory  cloud  generation  system  cannot 
reproduce  all  of  the  natural  processes  involved  in  the  actual  formation 
and  evolution  of  an  atmospheric  cloud,  an  adequate  laboratory  model  can 
be  created  by  a  system  possessing  the  following  characteristics: 

a.  Individual  control  of  the  following  parameters: 

(1)  Temperature 

(2)  Humidity 

(3)  Volume 

(4)  Pressure 

(5)  Dispersion  of  Condensation  Nuclei 

b.  Environmental  conditions  well  controlled,  instrumented,  and 
repeatable. 

c.  Minimum  interference  from  walls  of  container. 

It  is  our  purpose,  then,  in  this  section  to  describe  the  design, 
fabrication,  and  instrumentation  of  a  laboratory  cloud  chamber  and  the 
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corresponding  characteristics  of  the  simulated  uloud  used  for  the  cloud/ 
target  polarization  discrimination  experiments. 

5.2  THE  CLOUD  CHAMBER 

The  cloud  is  generated  in  a  chamber  constructed  of  modular  sections 
which  can  be  coupled  together  to  form  a  convenient  length.  Each  module 
was  fabricated  from  1/4-inch-thick  sheet  aluminum  in  the  form  of  a 
cube  with  two  opposite  sides  open.  The  internal  dimensions  of  each 
module  are  1  meter  x  1  meter  x  1  meter.  Aluminum  edge  strips  were 
welded  to  the  outside  of  the  open  ends  to  provide  rigidity  and  to  serve 
as  gasket  surfaces  to  obtain  an  airtight  and  watertight  seal  when  coup¬ 
ling  one  module  to  another  or  for  attaching  a  door.  Two  1-meter  x  1- 
meter  doors  were  fabricated  from  1/4-inch-thick  sheet  aluminum.  View¬ 
ing  ports  10  centimeters  in  diameter  were  installed  in  the  center  of  each 
door.  Additional  viewing  ports  were  installed  in  the  center  of  one  side 
of  the  modules  (see  Fig.  6).  The  inside  of  each  module  is  painted  with 
nonreflecting  flat  black  (camouflage)  paint  and  the  outside 
(top,  bottom,  and  sides)  is  insulated  with  two  layers  of  2-inch-thick 
headboard.  The  doors  are  similarily  insulat'd.  The  modules  are  sup¬ 
ported  on  a  specially  built  1-inch-thick  plywood-top  table  and  can  be 
coupled  together  (with  the  doors  attached  to  the  end  modules)  to  form 
a  variable-length  cloud  chamber. 

Each  module  is  provided  with  an  independent  water  vapor  injection 
system  at  the  bottom  and  independent  cooling  coils  at  the  top.  The 
water  vapor  in  the  form  of  steam  is  generated  outside  the  chamber  in  a 
steam  boiler  (Model  MB-6L,  Automatic  Steam  Products  Corp.,  N.Y.,  N.Y.). 
The  steam  is  injected  into  each  module  via  a  3/8-inch  O.D.  copper  tube 


97 


VIEWING  PORTS 
(10cm  OIA) 


CtPCO 

QiUJZ 

=>ce 

=J<cr 

Oujo 

XCOlr- 


u 


3 


98 


Figure  6.  Overall  Dimensions  of  Clo 


through  an  electric  solenoid  valve  (Type  GL-181,  Series  500,  Atkomatic 
Valve  Co.,  Inc.,  Indianapolis,  Indiana)  into  a  needle  valve  control. 

The  needle  valve  control  is  mounted  along  the  center  line  and  2  inches 
above  the  chamber  floor  in  the  wall  opposite  the  side  containing  the 
viewing  port.  An  independent  set  of  serpentined  cooling  coils  made  from 
1/2  inch  I.D.  copper  tubing  is  mounted  inside  of  each  module.  The 
coils  are  attached  directly  to  the  top  section  which  thus  acts  as  a 
giant  cooling  fin.  To  avoid  heating  the  laboratory  air  and  to 
maintain  minimum  noise  and  vibration  levels,  the  compressors,  1-1/2 
horsepower  heavy  duty  air-conditioning  units  (Model  AH-16H,  Dunham-Bush 
Inc.,  West  Hartford,  Connecticut),  are  mounted  outside  the  laboratory. 

The  refrigerant  (Freon-12)  is  pumped  to  each  module  through  3/8-inch  O.D. 
copper  tubing,  and  through  a  thermostatic  expansion  valve  (Type  CFE-1-1/2  C, 
Sporlan  Valve  Co.,  St.  Louis,  Missouri)  mounted  at  the  top  of  the  out¬ 
side  wall  opposite  the  side  with  the  viewing  port.  Combination  filter- 
drier  units  (Catch-All  Series,  £>prolan  Valve  Co.,  St.  Louis,  Missouri) 
are  installed  in  the  liquid  line  ahead  of  the  expansion  valve  and  in 
the  suction  line  to  remove  moisture  and  foreign  material  from  the  sys¬ 
tem.  The  thermostats  are  set  to  maintain  the  temperature  in  the  modules 
between  -8°C  and  -5°C. 

The  original  design  concept  for  the  cloud  chamber  provided  for  the 
injection  of  liquid  nitrogen  and  compressed  air  at  varying  rates  conflu¬ 
ent  with  the  stedm  to  (a)  offset  the  steam  heatload  and  (b)  provide  some 
measure  of  variability  in  particle  size  distributions  that  can  be  pro¬ 
duced.  Additionally,  the  compressed  air  injection  system  was  to  deliver 
(relatively)  dry  air  over  the  cooling  coils  to  prevent  Ice  from  forming 
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on  che  colls.  We  observed,  however,  during  the  Initial  operation  of  the 
chamber,  that  the  compressors  handled  the  steam  heatload  adequately  and 
variable  particle  size  distributions  could  be  generated  by  varying  the 
steam  injection  rates  and  the  steady  state  temperature  of  the  chamber. 

We  observed,  also,  that  blowing  compressed  air  over  the  cooling  coils  pro¬ 
duced  undesirable  turbulence  within  the  cloud.  We,  therefore,  omitted 
the  liquid  nitrogen  and  compressed  air  injection  systems  in  the  final 
design. 

Remote  sensing  thermometers  are  used  to  monitor  the  room  temperature 
and  the  coil  and  air  temperature  of  each  modular  section.  A  probe  type 
absorption  hygrometer  based  on  electrical  changes  produced  by  the  absorp¬ 
tion  of  water  on  conducting  materials  (Model  W611A,  Minneapolis  Honeywell, 
Minneapolis,  Minnesota)  is  used  to  monitor  the  relative  humidity  within 
the  chamber  while  a  membrane  type  absorption  hygrometer  (Bacharach  Indus¬ 
trial  Instruments  Co.,  Pittsburgh,  Pennsylvania)  is  used  to  monitor  the 
relative  humidity  within  the  laboratory. 

Provision  has  been  made  in  the  design  for  the  insertion  of  trans¬ 
parent  baffle  plates  between  sections  for  the  creation  of  different  con¬ 
ditions  in  selected  portions  of  the  chamber  to  more  realistically  simu¬ 
late  actual  clouds.  It  would  then  be  possible  to  have  the  cloud  local¬ 
ized  near  the  detector,  near  the  target,  or  at  some  intermediate  point. 

A  cut-away  view  of  a  modular  section  is  shown  in  Fig.  7. 

5.3  CLOUD  GENERATION  TECHNIQUES 

As  we  discussed  in  Section  II,  an  atmospheric  cloud  is  the  result  of 
either  lowering  the  temperature  of  a  volume  of  air  to  a  point  where  its 
capacity  for  water  vapor  is  greatly  reduced  or  of  adding  water  vapor  to  the 
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Figure  7,  Cut-Away  View  of  Cloud  Chamber  Module 


air  to  the  extent  necessary  to  increase  the  content  to  the  maximum  capa¬ 
city  possible  at  the  prevailing  temperature.  A  combination  of  these 
processes  may  also  result  in  the  formation  of  a  cloud. 

Two  basically  different  techniques  are  used  to  produce  artificial 
clouds  in  the  laboratory:  the  mechanical  technique  and  the  thermodynamic 
technique.  In  the  mechanical  process,  water  under  pressure  is  forced 
through  a  tiny  orifice  and  mixed  with  high  pressure  air.  The  mixture  is 
then  expanded  to  form  a  fine  spray.  In  the  thermodynamic  process,  one  of 
the  parameters  (pressure  or  temperature)  of  the  cloud  generation  system 
is  varied,  resulting  in  the  condensation  of  water  vapor  contained  in  the 
saturated  atmosphere  within  the  system.  Controlled  cloud  production 
utilizes  variations  or  combinations  of  these  two  methods.  We  chose  to 
use  a  combination  of  the  two  techniques  in  our  chamber  facility. 

We  begin  the  cloud  generation  process  by  sealing  off  the  chamber 
from  the  laboratory  environment  (primarily  accomplished  by  covering 
viewing  ports  with  quartz  plates)  and  then  lowering  the  temperature  of 
the  air  within  the  chamber  from  approximately  25°C  (room  temperature) 
to  -8°C  thereby  increasing  the  relative  humidity  within  the  chamber 
from  about  23  per  cent  to  60  per  cent.  During  this  phase,  the  pressure 
within  the  chamber  is  assumed  constant  and  equal  to  the  atmospheric 
pressure.  We  now  inject  steam  under  pressure  through  the  individual 
orifices  in  each  section.  The  transmissivity  of  the  cloud  is  continu¬ 
ously  monitored  and  recorded  on  a  strip  chart  recorder  (Hewlett-Packard 
Model  680)  using  a  100  mw  He-Ne  (helium-neon)  gas  laser  (Model  125, 
Spectra-Physics,  Palo  Alto,  California)  modulated  with  a  KDP  (potassium 
dihydro6en  phosphate)  crystal  at  a  pulse  repetition  frequency  of  1  kHz 
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that  is  beamed  into  a  photomultiplier  tube  over  a  one-way  path  through 
the  length  of  the  cloud  chamber.  Steady  state  conditions  appear  to  occur 
when  the  chamber  temperature  is  between  -5°C  and  0°C  and  the  relative 
humidity  in  the  chamber  reaches  90  percent.  During  this  steady  state 
phase  we  have  been  able  to  hold  the  transmissivity  to  20  percent  for  two 
hours  with  only  five  to  ten  percent  fluctuations.  The  corresponding 
droplet  diameters  as  monitored  with  microscopic  photography  techniques 
indicate  particle  size  distributions  which  peak  (mode  diameter)  at  5 
to  6  microns  with  mean  droplet  diameters  ranging  from  8.9  to  9.1  microns 
in  good  agreement  with  the  distributions  of  Table  1  for  stratus  type 
clouds  (includes  fog  which  is  defined  as  a  stratus  cloud  that  touches 
the  ground).  For  the  temperature  range  between  0  and  +.5°C  the  droplet 
size  distribution  indicates  that  the  cloud  resembles  a  warm  temperature 
coastal  fog  (mode  diameter:  8-10  microns,  mean  diameter:  12  -  15 
microns).  The  transmissivity  of  the  cloud  is  greater  (40  -  60  percent) 
for  this  phase  indicating  a  less  dense  cloud  which  is  another  character¬ 
istic  of  a  coastal  fog.  Additionally,  the  transmissivity  fluctuates 
considerably  over  long  periods  of  time  and  is  impossible  to  hold  con¬ 
stant  for  longer  periods  than  15  minutes. 

The  original  design  concept  for  the  chamber  provided  for  the  injec¬ 
tion  of  condensation  nuclei  to  facilitate  cloud  formation.  However,  we 
have  had  phenomenal  results  using  the  nuclei  present  in  the  nat¬ 
ural  air  in  the  chamber  and  carried  by  the  steam.  The  existence  of 
nuclei  is  indicated  by  the  fact  that  the  relative  humidity  within  the 
chamber  never  Increases  above  92  percent.  If  the  process  of  producing 
a  cloud,  and  then  allowing  it  to  dissipate,  is  repeated  several  times, 
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the  resulting  clouds  become  successively  less  dense,  requiring  increas¬ 
ingly  longer  time  intervals  in  which  to  form. 

5.4  CLOUD  SAMPLING  TECHNIQUES 

It  is  surprising,  indeed,  that  while  much  effort  has  been  devoted 
to  developing  experimental  techniques  for  measuring  the  size  and  concen¬ 
tration  distributions  of  cloud  droplets,  we  still  do  not  have  an  entirely 
satisfactory  technique  for  accurately  and  reliably  measuring  the  size 
distribution  without  physically  disturbing  the  cloud  and  that  does  not 
require  extremely  tedious  drop-by-drop  measurements.  The  following  are 
the  most  commonly  used  methods  used  for  measuring  droplet  sizes  within 
a  cloud  or  fog: 

a.  An  optical  technique  based  on  the  diffraction  theory  of  Airy 
[Ref.  1]  which  depends  on  the  measurement  of  the  distribution  of  light 
intensity  (coronas)  produced  by  the  diffraction  of  a  beam  of  light. 

Wilson  [Ref.  59]  has  argued  mathematically  that  this  method  of  sampling 
favors  the  larger  diameter  droplets  (>20  microns)  and  is  useful  only 

in  homogeneous  fogs. 

b.  Kneusel  [Ref,  30]  built  an  apparatus  to  observe  the  rate  of 
fail  of  atmospheric  fog  particles  in  undisturbed  air.  Assuming  this 
velocity  of  fall  was  due  to  the  earth's  gravitational  force  of  attrac¬ 
tion,  he  calculated  the  corresponding  droplet  sizes.  His  distributions, 
however,  appear  to  favor  the  smaller  droplets. 

c.  An  electrostatic  precipitation  method  first  devised  by 
Panthenier  and  Brun  (Ref.  46]  uses  a  metal  cylinder  lined  with  a  special 
absorbent  paper  dyed  with  methylene  blue.  The  droplets  within  the  cylin¬ 
der  are  precipitated  onto  the  paper  by  an  electrostatic  field  and 


104 


produce  impressions,  the  diameters  of  which  are  proportional  to  the 
actual  particle  diameters. 

d.  Kohler  [Ref.  31]  attempted  to  collect  fog  droplets  on  fine 
wire  but  was  unable  to  correct  for  the  effects  jf  evaporation.  On  the 
other  hand,  Dessens  [Ref.  14]  and  later  Arnulf  and  Bricard  [Ref.  3]  ob¬ 
tained  results  (see  Table  1)  using  spiders'  webs. 

e.  Probably  the  most  reliable  technique  employed  by  the  majority 
of  researchers  involves  the  collection  of  the  droplets  on  suitably  pre¬ 
pared  surfaces  and  the  subsequent  photographing  of  the  droplets  or  their 
impressions  under  a  darkfield  microscope.  In  one  method  employed  by 
many  researchers,  a  microscope  slide  is  coated  with  a  greasy  substance 
(usually  petroleum  jelly)  while  in  an  equally  popular  method,  the  slide  is 
covered  with  oil  (mineral  oil,  kerosene,  etc.).  According  to  Houghton 
[Ref.  22]  droplets  caught  on  slides  coated  with  petroleum  jelly  flatten  out 
to  the  point  where  the  droplet  radius  measured  from  the  slide  is  al¬ 
most  twice  as  large  as  that  of  the  original  droplet;  those  caught  on 
slides  covered  with  oil,  however,  sink  and  remain  essentially  spherical 
[Ref.  8].  An  alternate  method  is  to  coat  a  slide  with  magnesium  oxide 

in  which  the  droplets  produce  impressions,  the  diameters  of  which  are 
proportional  to  the  actual  particle  diameters. 

f.  Laser  holographic  techniques  have  been  developed  that  provide 
a  three-dimensional  view  of  fog  droplets  in  a  specific  volume  of  space. 
However,  the  use  of  this  technique  is  limited  in  application,  due  to 
stringent  stability  and  alignment  requirements. 

g.  Aerosol  counters  are  coonerciai ly  available  that  determine  the 
site  distribution  and  concentration  of  aerosol  particles  by  measuring 
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the  intensity  of  the  light  scattered  by  a  single  particle  which  is  ill¬ 
uminated  by  a  visible  light  source. 

After  careful  consideration  of  these  methods  with  regard  to  the 
availability  of  resources  (money,  personnel,  equipment,  materials,  lab¬ 
oratory  rpace,  etc.)  and  the  reliability  and  accuracy  of  the  resulting 
data,  we  chose  to  measure  and  compare  the  size  distribution  by  two  sep¬ 
arate  methods:  microscopic- slide  photography  and  a  light-scattering 
aerosol  counter.  The  aerosol  counter  has  the  additional  capability  of 
providing  concentration  data. 

5.4.1  MICROSCOPIC  SAMPLING  TECHNIQUE 

The  microscopic  sampling  technique  consists  of  thoroughly  cleaning 
suitable  microscope  slides  and  then  applying  melted  petroleum  jelly  to  the 
slide  surface  to  form  a  very  thin  film.  A  drop  of  mineral  oil  or  kero¬ 
sene  is  next  applied  to  the  jelly  surface.  The  slide  is  cooled  to 
the  ambient  temperature  of  the  chamber  and  then  submerged  in  the  cloud 
for  a  finite  time.  The  droplets  sink  through  the  kerosene  and  barely 
rest  on  the-  petroleum  jelly,  the  buoyancy  of  the  kerosene  preventing  them 
from  flattening  out  on  or  penetrating  into  the  jelly.  The  validity  of 
this  technique  is  indicated  (when  viewing  a  sample  under  the  microscope) 
by  the  fact  that  while  the  droplets  do  not  appear  to  evaporate  they  do 
wander  and  collide  to  some  extent,  The  corresponding  degree  of  coales¬ 
cence  is,  however,  very  small  with  never  a  complete  merging  of  droplets; 
that  is,  occssionally  a  smaller  droplet  will  appear  to  be  "stuck"  to  s 
larger  droplet.  The  slide  is  removed  from  the  chamber  and  photographed 
under  a  darkfield  microscope  (Bausch  and  Lomb  Oynaxoom  Flat  Field  Bi¬ 
nocular  Body  Microscope,  Bausch  and  Lomb,  Rochester,  New  York).  <*e 
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have  observed,  using  the  petroleum  Jelly  kerosene  mixture,  that  the  captured 
droplets  do  not  begin  to  evaporate  or  coalesce  for  about  30  minutes  at 
room  temperature.  If  maintained  at  the  ambient  temperature  of  the 
chamber,  their  size  and  concentration  distributions  remain  essentially 
unchanged.  The  photomicrograph  is  then  magnified  under  a  measuring  mi¬ 
croscope  (Carl  Zeiss  Tool  Makers  Microscope)  and  the  diameter  of  each 
particle  is  measured  with  an  accuracy  of  1  micron  in  125  millimeters. 

The  measured  values  of  the  droplet  diameters  are  then  sorted  into  size 
groups  and  recorded  in  increments  of  0.1  micron. 

In  determining  the  size  distribution  and  the  corresponding  degree 
of  repeatability  of  our  artificial  cloud,  we  compared  the  data  collected 
by  two  different  sampling  techniques: 

a.  In  the  first  method  we  prepare  seven  slides  with  the  petroleum 
jelly  kerosene  mixture  in  the  manner  described  in  the  preceding  paragraph 
and  then  place  them  in  a  narrow  covered  tray  extending  across  the  cross 
sectional  width  of  the  chamber.  The  slides  are  allowed  to  cool  to  the 
ambient  temperature  of  the  chamber  before  being  exposed  to  the  cloud. 

When  steady  state  conditions  exist  within  the  cloud  for  a  particular 
transmissivity,  the  tray  is  uncovered  and  the  slides  are  exposed  to  the 
cloud  for  a  finite  time.  At  the  end  of  this  sampling  time  the  first 
slide  is  removed  from  the  chamber  and  immediately  photographed  under  the 
microscope.  The  remaining  slides  are  left  in  the  chamber  with  their 
surfaces  exposed  to  the  cloud  and  are  subsequently  removed  and  photo¬ 
graphed  in  turn. 

This  sampling  method  was  used  to  collect  data  for  clouds  of  differ¬ 
ent  densities  for  the  temperature  ranges  -5°  to  0°C  and  0°  to  +5°C. 
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The  density  of  the  individual  clouds  was  continuously  monitored  as  a 
function  of  the  transmissivity.  While  the  individual  size  distributions 
did  not  appear  to  vary  as  a  function  of  density  (as  indicated  by  the 
transmissivity)  within  either  of  the  temperature  intervals,  there  was  a 
definite  difference  in  size  distributions  and  corresponding  densities 
of  the  clouds  obtained  for  the  two  temperature  intervals.  For  the  tem¬ 
perature  interval  -5°  to  0°C,  the  mode  diameter,  d^,  was  5  microns  and 
the  corresponding  mean  diameter,  dm,  was  8.9  microns,  while  for  the  tem¬ 
perature  interval  0°  to  +5°C,  d,  was  8  microns  and  d  was  12.3  microns. 

d  m 

In  the  lower  temperature  interval  the  transmissivity  ranged  from  20  to 
100  percent  while  in  the  higher  temperature  interval  the  transmissivity 
ranged  from  40  to  100  percent.  We  conclude,  then,  from  the  increase  in 
droplet  size  coupled  with  a  decrease  in  density  (indicated  by  an  increase 
in  transmissivity)  that  the  density  of  the  clouds  in  the  temperature  in¬ 
terval  above  0°C  is  less  than  that  of  the  clouds  in  the  interval  below 
0°C.  The  resultant  size  distributions  for  arbitrary  density  (trans¬ 
missivity)  as  a  function  of  temperature  are  shown  in  Figs.  8  and  9. 

b.  In  the  second  sampling  technique  the  slide  surfaces  are  pre¬ 
pared  and  exposed  to  the  cloud  in  the  same  manner  as  previously  de¬ 
scribed.  But  at  the  end  of  the  sampling  interval  when  the  first  slide 
is  removed  for  photographing  the  remaining  slides  are  covered  and 
thusly  isolated  at  the  ambient  temperature  of  the  chamber  and  are  sub¬ 
sequently  removed  and  photographed. 

This  sampling  method  was  used  to  collect  data  as  a  function  of 
cloud  density  and  chamber  temperature  in  the  same  manner  as  in  the 
first  method.  Again,  the  results  indicated  a  temperature  dependence 
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while  «  variation  in  the  density  (transmissivity)  did  not  indicate  a 
change  in  size  distribution.  The  resultant  sir.v  distributions  are  shown 
in  Figs.  10  and  11.  For  the  temperature  interval  -5°  to  0°C,  d^  *  6  mi¬ 
crons  and  d  ■  9.2  microns:  for  the  interval  0°  to  +5°C,  d,  *  10  microns 

m  a 

and  d  ■  15  microns, 
m 

Examples  of  the  photomicrographs  from  which  the  temperature  depen¬ 
dent  size  distributions  were  derived  are  shown  in  Figs.  12  and  13. 

A  comparison  of  the  results  obtained  by  the  two  sampling  techniques 
show  that,  while  the  mode  diameters  differ  by  about  20  percent,  the  mean 
diameters  differ  by  only  3.5  percent.  We  assume  that  the  distributions 
obtained  by  the  first  method  more  realistically  represent  the  actual 
cloud  because  almost  twice  as  many  particles  were  collected  by  the  first 
method.  The  sampling  intervals  (based  on  the  first  slide)  were  equal 
for  the  two  methods  and  when  this  interval  was  doubled  using  the 
second  method  of  sampling  in  subsequent  experiments,  the  resultant  dis¬ 
tributions  corresponded  to  those  obtained  by  the  first  method. 

F ARTICLE  MONITOR-SAMPLING  TECHNIQUE 
The  particle  monitor-sampling  technique  consists  of  using  a  light¬ 
scattering  aerosol  counter  for  determining  the  size  distribution  and 
concentration  of  aerosol  particles.  A  light- scattering  aerosol  counter 
(Model  202,  Royco  Instruments,  Menlo  Park,  California)  determines  the 
size  distribution  and  number  concentration  of  aerosol  particles  by 
measuring  the  intensity  of  visible  light  scattered  by  a  single  particle 
in  the  range  of  scattering  angles  from  75°  to  105°.  The  air  (containing 
the  aerosol  particles)  being  monitored,  is  blown  through  a  volume 
illuminated  by  the  radiation  of  a  white  light  source  of  color  temperature 
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Figure  12.  Photomicrograph  (X400)  of  Laboratory 
Generated  Cloud  Droplets 
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Figure  13.  Photomicrograph  (X400)  of  Laboratory 
Generated  Coastal  Fog  Droplets 
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2800°  K.  The  detector  is  a  photomultiplier  tube.  Each  particle 
remains  within  the  illuminated  volume  for  a  short  time  only  and  the 
magnitude  of  its  electrical  signal  depends  on  its  size.  The  signals 
are  electronically  classified  according  to  their  magnitude  in  a  narrow 
interval.  The  number  of  the  impulses  in  each  interval  is  counted  and 
printed  out  in  definite  time  intervals. 

The  size  distributions  obtained  from  the  aerosol  counter  differed 
greatly  from  those  obtained  from  the  microscope-slide  photography  tech¬ 
nique.  In  fact,  the  particle  counter  scarcely  indicated  the  presence  of 
any  particles,  while  at  the  same  time  approximately  10,000  droplets  were 
collected  on  the  slides.  We  believe  the  reason  that  the  particle 
counter  does  not  give  a  true  representation  of  the  size  distribution  is 
because  the  droplets  evaporate  before  they  can  be  measured  aue  to 
the  large  temperature  difference  between  th»  droplet  and  the  interior 
of  the  counter.  The  validity  of  this  argumenc  is  indicated  by  the  fact 
that  the  counter  appears  to  operate  satisfactorily  when  counting  "steam 
particles,"  i.e.,  for  steam  injected  into  the  chamber  when  the  tempera¬ 
ture  of  the  chamber  was  approximately  25°C  (room  temperature)  before 
sceam  infection. 
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SECTION  VI 


EXPERIMENTAL  VERIFICATION  OF  THE  CLOUD/TARGET 
POLARIZATION  DISCRIMINATION  TECHNIQUE 

6.1  INTRODUCTION 

In  this  chapter  we  present  the  results  of  our  experimental  research 
program  to  determine  the  validity  of  the  cloud/target  polarization  dis¬ 
crimination  technique.  The  program  is  divided  into  a  laboratory  veri¬ 
fication  phase  and  a  field  evaluation  phase.  In  the  laboratory  verification 
phase  of  the  investigation  we  measured  the  signal  returns  from  a 
reference  target  and  various  clouds  generated  in  the  cloud  chamber 
facility  described  in  Section  V  as  a  function  of  cloud  transmissivity. 

The  polarization  properties  of  the  return  signals  were  then  compared 
using  the  plane  of  vibration  of  the  electric  field  vector  of  the  trans¬ 
mitted  wave  as  a  reference.  In  this  respect  we  define  the  following 
symbols  assuming  the  transmitted  wave  to  be  linearly  polarized: 

"+"  signifies  unpolarized  radiation 

"ll"  signifies  radiation  whose  electric  field  vector  vibrates 
parallel  to  that  of  the  tx-ansmitted  wave 
"l"  signifies  radiation  whose  electric  field  vector  vibrates 
orthogonally  to  that  of  the  transmitted  wave 
A  qualitative  verification  of  the  polarization  discrimination 
technique  is  depicted  in  the  series  of  photographs  shown  in  Figs.  15  and 
17.  These  photographs  present  a  visual  record  of  the  scattering  phenomena 
we  observed  when  the  output  of  a  linearly  polarized  He-Ne  (helium-neon) 
gas  laser  was  transmitted  through  a  laboratory  generated  cloud  and  was 
incident  on  a  diffusely  reflecting  target.  The  wavelength,  X,  of  the 
laser  radiation  was  in  the  visible  (red)  portion  of  the  electromagnetic 
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3  in. 


spectrum  ;-nd  equal  to  6328  Angstroms  (10**  Angstroms  =  1  micron).  « 
x  5  in.  piece  of  fine  grain  sandpaper  sprayed  with  white  paint  was  used 
as  the  diffuse  target.  A  33mm  Nikon  F  camera  was  used  to  photograph  the 
radiation  scattered  from  the  cloud  and  the  target.  An  analyzer  was  posi¬ 
tioned  in  front  of  the  camera  and  mounted  so  that  it  could  be  easily  removed 
and  replaced  in  its  original  position.  Schematic  diagrams  of  the  exper¬ 
imental  arrangements  used  in  obtaining  the  phonographic  illustrations  of 
Figs.  15  and  17  are  shown  in  Figs.  14  and  16,  respectively. 

The  experimental  arrangement  illustrated  in  Fig.. 14  was  used  to 
obtain  the  photographs  shown  in  Fig.  15.  For  this  situation  the  laser 
beam  was  transmitted  along  the  long  dimension  of  the  cloud  chamber  and 
impinged  upon  the  diffuse  target  which  was  placed  at  the  far  end  of  the 
chamber.  The  camera  was  positioned  coaxial  to  that  portion  of  the  trans¬ 
mitted  beam  being  propagated  through  the  cloud  chamber.  The  analyzer  was 
not  used  for  the  photograph  shown  in  Fig.  15a.  The  large  circle  in  the 
center  of  the  photograph  is  the  viewing  port.  The  smaller  red  circle 
within  the  viewing  port  represents  the  sum  of  (a)  the  radiation  reflected 
from  the  target  at  the  far  end  of  the  chamber  and  (b)  the  radiation  back- 
scattered  from  the  cloud  along  the  length  of  the  chamber.  The  diameter 
of  the  viewing  port  is  10  centimeters;  and  while  the  diameter  of  the  trans¬ 
mitted  laser  beam  is  only  one  centimeter,  the  diameter  of  the  red  circle 
is  approximately  six  centimeters.  The  apparent  increase  in  the  diameter  of 
the  transmitted  beam  was,  in  reality,  caused  by  a  halo  of  scattered  light 
surrounding  the  beam.  This  halo  of  light  was  primarily  composed  of  the 
radiation  scattered  out  of  the  transmitted  beam  by  the  cloud  particles. 

The  black  annular  ring  surrounding  the  red  circle  represents  that  portion 
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Figure  15.  Photographic  Illustration  of  the  Cloud/Target  Polarization  Dis¬ 
crimination  Technique,  (a)  Laser  Output  Linearly  Polarized,  No 
Analyzer  on  Camera;  (b)  Laser  Output  Linearly  Polarized,  Plane 
of  Polarization  of  Analyzer  (Located  at  Camera)  Orthogonal  to 
That  of  Transmitted  Beam. 


Figure  17.  Photographic  Illustration  of  the  Cloud/Target  Polarization  Dis¬ 
crimination  Technique,  (a)  Laser  Output  Linearly  Polarized,  No 
Analyzer  on  Camera;  (b)  Laser  Output  Linearly  Polarized,  Plane 
of  Polarization  of  Analyzer  (Located  at  Camera)  Orthogonal  to 
That  of  Transmitted  Beam. 

Beam  enters  at  left  side  of  each  photograph.  Target  holder 
visible  in  lower  part  of  each  photograph. 


of  the  cloud  (as  seen  through  the  viewing  port)  that  is  not  illuminated 
by  the  laser  beam.  The  photograph  shown  in  Fig.  15b  was  taken  with  the 
analyzer  positioned  in  front  of  the  camera  and  oriented  so  that  its 
plane  of  polarization  was  perpendicular  to  the  plane  of  polarization  of 
the  transmitted  beam.  The  small,  irregularly  shaped  red  spot  in  the 
center  of  the  photograph  represents  the  sum  of  (a)  those  components  of 
the  radiation  reflected  from  the  target  whose  plane  of  polarisation  was 
orthogonal  to  that  of  the  transmitted  wave,  (b)  a  negligible  amount  of 
multiple  scattered  radiation  whose  plane  of  polarization  is  orthogonal 
to  that  of  the  transmitted  wave,  and  (c)  a  negligible  amount  of  radiation 
whose  plane  of  polarization  is  parallel  to  that  of  the  transmitted  wave 
(some  of  the  radiation  with  parallel  polarization  is  not  completely 
attenuated  by  the  non  perfect  analyzer). 

A  more  vivid  and  explicit  representation  of  the  discrimination 
technique  is  depicted  in  the  photographs  shown  in  Fig.  17.  The  experi¬ 
mental  arrangement  illustrated  in  Fig.  16  was  used  to  obtain  the  photo¬ 
graphs  shown  in  Fig.  17.  For  this  situation  the  laser  beam  was  trans¬ 
mitted  along  the  long  dimension  of  the  cloud  chamber  and  impinged  upon 
the  diffuse  target  which  was  placed  opposite  a  side  viewing  port.  The 
target  was  mounted  on  a  tubular  holder  and  inclined  at  an  angJe  of  45° 
with  respect  to  the  transmitted  beam.  The  camera  was  positioned  opposite 
the  side  viewing  port  as  shown  in  Fig.  16.  Again,  as  in  the  previous 
set  of  photographs,  the  analyzer  was  placed  in  front  of  the  camera  with 
its  plane  of  polarization  perpendicular  to  that  of  the  transmitted  wave. 
The  transmitted  beam  enters  at  the  left  side  of  each  photograph  (not 
visible  in  Fig.  17b).  The  target  holder  is  visible  in  the  lower  part  of 
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each  photograph.  In  Fig.  17a  we  see  that  when  no  polarizer  was  used, 
the  intensity  of  the  scattered  radiation  completely  concealed  the  target 
from  view.  However,  in  fig.  17b  we  observe  that  the  use  of  the  analyzer 
attenuated  the  scattered  radiation  seen  by  the  camera  to  such  an  extent 
that  the  target  is  partially  visible. 

The  series  of  photographs  shown  in  Figs.  15  and  17  clearly  verify 
the  validity  of  the  polarization  discrimination  concept.  According  to 
the  Mie  theory  of  scattering,  if  a  linearly  polarized  electromagnetic 
plane  wave  is  incident  on  a  spherical  particle  of  arbitrary  composition 
and  size,  then  the  scattered  radiation  is  also  linearly  polarized  for 
all  scattering  angles  (i.e.  0 8  i  0  <  180°)  in  the  XZ  and  YZ  planes 
defined  in  Fig.  1.  On  the  other  hand,  the  radiation  reflected  (scattered) 
from  a  diffuse  (rough)  surface  will,  in  general,  be  depolarized.  For  the 
photographs  shown  in  Fig.  15  the  camera  was  in  the  XZ  plane  and  0  *  180°. 

We  see,  then,  that  the  observed  scattering  phenomena  shown  in  Fig.  15a 
and  b  are  in  good  agreement  with  the  Mie  and  rough  surface  theory  of 
scattering.  For  the  photographs  shown  in  Fig.  17  the  camera  was  in  the 
YZ  plane;  and  while  the  point  of  observation  was  on  a  line  normal  to 
the  direction  of  propagation  of  the  transmitted  wave,  we  can  assume 
that  the  scattering  angle  varied  over  a  wide  range  of  values.  This 
can  be  seen  from  a  comparison  of  Fig.  17a  and  Fig.  17b.  Figure  17a  is 
a  photograph  of  the  sum  of  the  individual  electromagnetic  waves  scattered 
out  of  the  beam  (within  the  field  of  view  of  the  camera)  along  the  direction 
of  propagation  of  the  transmitted  signal.  Each  of  these  component  waves 
had  a  different  direction  of  propagation  (scattering  angle)  in  the  YZ 
plane  which  was  determined  by  the  location  and  size  of  the  scattering 
particles.  As  shown  in  Fig.  17b,  the  analyzer  (its  plane  of  polarization 
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crossed  with  respect  to  that  of  the  transmitted  signal)  completely  atten¬ 
uated  the  scattered  radiation  surrounding  the  beam  while  transmitting  the 
perpendicular  component  of  the  radiation  reflected  from  the  target.  We  can 
assume,  then,  that  the  scattered  radiation  was  of  the  same  polarization  as 
that  of  the  transmitted  signal  and  wus  independent  of  the  scattering  angle 
(direction  of  observation).  This  series  of  photographs,  then  presents  a 
qualitative  justification  for  our  experimental  program. 

Originally,  the  experimental  verification  phase  of  this  effort  was 
to  be  performed  primarily  in  laboratory  and  if  the  results  warranted, 
would  be  expanded  to  include  a  flight  test  program  in  the  field  evaluation 
phase  of  the  effort.  However,  since  the  initiation  of  the  laboratory 
effort,  an  aerial  tramway  facility  was  constructed  in  the  Sandia  Mountains 
at  Albuquerque,  New  Mexico.  This  facility  has  provided  an  exceptional 
opportunity  for  performing  laboratory  experiments  in  an  actual  atmospheric 
environment.  The  program  was  justifiably  expanded  to  include  tramway 
experiments  and  the  results,  to  date,  are  presented  in  section  6.3. 

6.2  CLOUD  CHAMBER  EXPERIMENTS 

The  experimental  arrangement  shown  in  Fig.  18  was  used  for  the  He-Ne 
laser  polarization  discrimination  experiments.  The  well  collimated  light 
(A  »  6328  Angstroms)  output  from  a  He-Ne  gas  laser  was  beamed  into  an 
electro-optic  (KDP  crystal)  modulator,  the  output  of  which  was  directed 
through  a  Glan-Thompson  prism  and  special  mirror  into  and  along  the  long 
dimension  of  the  cloud  chamber.  Even  though  the  output  of  the  laser  is 
linearly  polarized,  we  used  the  Glan-Thompson  prism  to  ensure  a  high 
degree  of  polarization.  The  ratio  of  the  Intensities  of  the  orthogonally 
polarized  and  parallel  polarized  components  (I^/I^)  at  the  output  of 
the  Glan-Thompson  prism  was  measured  to  be  5.2  x  10  .  The  special  mirror 
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Figure  18.  Experimental  Arrangement  for  He-Ne  Laser  Cloud  Scattering  Experiments 


is  just  an  optical  flat  with  a  hole  cut  through  it  to  allow  the  passage 
of  the  transmitted  beam.  The  intensity  of  the  transmitted  beam  is  moni¬ 
tored  at  the  target  end  of  the  chamber  by  a  collimated  transmission 
detector  (RCA  PMT  No.  7102)  and  is  continuously  displayed  on  a  strip-chart 
recorder  (Hewlett-Packard  Model  680) .  Two  beam  splitters  are  used  to 
attenuate  (by  reflection)  the  signal  so  as  to  prevent  any  saturation  of 
the  detector.  They  each  have  a  rated  transmission  of  5  percent.  The  black 
felt  absorber  is  inclined  at  45°  with  respect  to  the  direction  of  propa¬ 
gation  of  the  reflected  beams  and  is  used  to  prevent  the  energy  reflected 
by  the  beam  splitters  from  entering  the  detection  system.  The  reflected 
radiation  from  the  target  and  the  backscattered  radiation  from  the  cloud 
were  collected  on  the  special  mirror,  collimated  by  the  thin  lenses,  and 
directed  through  another  Glan-Thompson  prism  onto  the  photocathode  of  an 
S-l  photomultiplier  tube  (RCA  No.  7102)  the  output  of  which  was  monitored 
and  photographed  on  a  Tektronix  Model  454  oscilloscope. 

Ue  chose  a  MgO  (magnesium  oxide)  coated  surface  as  a  standard  target 
because  of  the  high  and  almost  perfectly  diffuse  reflectance  of  MgO.  For 
our  particular  MgO  target,  the  coefficient  of  reflection  p  »  0.98  and  the 
average  value  of  the  ratio  of  the  measured  intensities  of  the  orthogonally 
polarized  and  parallel  components  Cl was  equal  to  0.324.  In  all 
experiments  where  the  target  was  used,  the  target  intercepted  all  of  the 
incident  bean  (which  was  well  collimated  at  the  input  of  the  chamber). 

For  every  experluent  the  total  intensity  (1^)  the  individual 
intensities  of  the  orthogonally  polarised  (1^)  *f>d  parallel  polarized 
(I  )  components  of  the  radiation  scattered  from  the  target  were  measured 
before  generating  a  cloud,  and  again  after  the  cloud  had  completely  dis¬ 
sipated.  The  average  values  of  these  measured  intensities  were  used  as  a 
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standard  reference  for  all  cloud  return  measurements.  In  measuring  the 
total  intensity,  I  . ,  Glan-Thompson  prism  No.  2  was  removed  and  the 
measured  value  had  to  be  corrected  for  a  7  percent  loss  In  the  prism. 

In  preliminary  experiments  we  attempted  to  determine  the  relation¬ 
ship  between  the  target  return  and  the  transmissivity  of  the  cloud.  Theo¬ 
retically,  the  transmissivity,  t,  is  defined  from  Equation  f68)  as 
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Thus,  if  the  target  is  a  perfect  diffuser,  we  see  from  Section  IV  and 
specifically  Equation  (285)  that  the  received  signal  (assuming  the  receiver 
is  located  adjacent  to  the  transmitter)  is  proportional  to  t2.  The  exper¬ 
imental  verification  of  the  theory  is  shown  in  Fig.  19.  The  plane  of 
polarization  of  Glan-Thompson  prism  No.  2  was  oriented  orthogonal  to  that 
of  prism  No.  1  to  reduce  tne  cloud  return  to  a  minimum.  The  transmissivity 
of  the  cloud  was  continuously  monitored  and  when  steady-state  conditions 
were  reached  for  a  specific  value  of  transmissivity,  the  MgO  target  was 
positioned  as  shown  in  Fig.  18  (which  causes  a  momentary  interruption 
in  the  monitoring  of  the  transmission)  and  the  amplitude  of  the  target 
return  was  recorded  from  the  oscilloscope.  The  target  was  then  removed 
and  the  return  from  the  cloud  was  recorded  along  with  the  transmissivity 
as  read  from  the  strip-chart  recorder.  Assuming  a  negligible  contri¬ 
bution  due  to  multiple  scattering,  the  log  of  the  difference  in  the 
oscilloscope  readings  is  plotted  against  the  log  of  the  corresponding  t. 

It  is  evident  from  the  resultant  slope  that  the  experimental  data  are 
in  close  agreement  with  theory  and  therefore  we  used  the  value  t 2  for 
subsequent  calculations  and  comparisons. 
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TRANSMISSIVITY  (t) 


Figure  19.  Target  Return  vs.  Cloud  Transmissivity 


For  the  polarization  discrimination  experiments  reference  measurements 

(It+,  ItA,  I  )  of  the  target  return  were  made  as  previously  discussed 

and  then  the  target  was  removed.  A  cloud  was  generated  and  the  individual 

pulse  returns  corresponding  to  the  total  intensity  I  and  the  component 

Intensities  I  .  and  I  .  of  the  radiation  scattered  from  the  clcud  were 
cl  cii 

photographed  on  the  oscilloscope.  The  corresponding  value  of  x  was 
appropriately  marked  on  the  strip  chart.  The  amplitudes  of  the  photo¬ 
graphed  returns  (both  target  and  cloud)  were  measured,  tabulated,  and 
corrected  for  Gian  prism  losses  when  appropriate.  All  values  were  then 
normalized  to  a  clear  air  (x  =  1)  target  return  where  It+  ■  1,  =  0./55, 

and  I  ■  0.245.  These  ratios  were  calculated  from  the  relationships 
I  ,  =  I  +  I  and  I  /I  =  G.324.  Plots  of  log  I  vs.  log  x  for  the 

»  '  LX  t  II  l»X  L|| 

cloud  returns  and  their  corresponding  target  returns  are  shown  in  Figs. 

20  through  22.  A  composite  plot  is  shown  in  Fig. 

It  is  easily  seen  from  Figs.  20,  21,  and  23  that,  when  no  analyzer 
was  used  on  the  receiver  or  when  the  plane  of  polarization  of  the  analyzer 
was  parallel  to  that  of  the  transmitted  signal,  the  cloud  return  was 
greater  than  the  target  return  for  x  <  0.5  (approximately).  On  the 
other  hand,  even  though  the  intensity  of  the  target  return  was  much 
lower  when  the  plane  of  polarization  of  the  analyzer  was  orthogonal  to 
that  of  the  transmitted  signal,  the  cloud  return  was  reduced  over  an 
order  of  magnitude  and  was  greater  than  the  target  return  only  when 
x  <  0.28.  We  see  from  Figs.  22  and  23  that  for  decreasing  values  of  r 
there  was  a  corresponding  increase  in  I  .  This  continuous  increase  in 
!ci  was  due  to  the  partial  depolarization  of  the  scattered  radiation 
caused  by  an  increase  in  multiple  scattering  at  the  higher  cloud  densities. 
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SIGNAL  RETURNS  (RELATIVE  SCALE) 


Figure  20.  Target  and  Cloud  Returns  vs.  Cloud  Transmissivity 
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SIGNAL  RETURNS  (RELATIVE  SCALE) 


Figure  21.  Target  and  Cloud  Returns  vs.  Cloud  Transmissivity 
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Figure  22.  Target  and  Cloud  Returns  va.  Cloud  Transmissivity 
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Figure  23.  Target  end  Cloud  Returns  vs.  Cloud  Transnissivlty 


It  is  apparent,  though,  that  I  approaches  an  asymptotic  limit  as  the 
density  increases  indicating  that  the  cloud  is  beginning  to  resemble  a 
solid  rather  than  an  extended  target. 

The  experimental  arrangement  shown  in  Fig.  24  was  used  for  the 
Ruby  laser  (X  *  6943  Angstroms)  polarization  discrimination  experiments. 
We  were  not  able  to  obtain  any  meaningful  data  from  this  phase  of  the 
program,  however,  because  the  amplitude  of  the  transmitted  signal 
varied  erratically  from  pulse-to-pulse.  The  amplitude  difference 
varied  as  much  as  20  percent.  This  phase  of  the  program  was  postponed 
until  the  source  of  the  malfunction  can  be  located  and  corrected. 

6.3  SANDIA  PEAK  AERIAL  TRAMWAY  EXPERIMENTS 

The  Sandia  Peak  Aerial  Tramway  provides  an  exceptional  capability 
for  performing  laboratory  experiments  in  an  actual  atmospheric  environ¬ 
ment.  It  is  the  longest  aerial  tramway  in  North  America,  the  distance 
between  the  lower  and  upper  terminals  being  2.7  miles.  The  difference 
in  elevation  between  the  lower  (elevation:  6550  feet)  and  upper  (ele¬ 
vation:  10378  feet)  terminals  is  3828  feet.  Two  60-passenger  cabins 
move  simultaneously  on  separate  cables  but  in  opposite  directions  (one 
moves  toward  the  upper  terminal  as  the  other  approaches  the  bottom) . 

The  distance  between  a  cabin  and  the  terrain  directly  below  varies  con¬ 
siderably  (from  50  to  1000  feet)  as  the  cabin  travels  next  to  and  ovei 
rocky  granite  spires,  cliffs,  and  pinnacles,  across  gullies  and  ravines 
formed  by  the  spring  runoffs,  and  over  steep  and  heavily-forested 
canyons  and  valleys.  In  traveling  from  the  base  to  the  top  of  the 
mountain,  the  tramway  passes  over  four  of  the  earth’s  seven  climato¬ 
logical  life  zones.  The  corresponding  species  of  trees  vary  from 
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Figure  24.  Experimental  Arrangement  for  Ruby  Laser  Cloud  Scattering  Experiments 


clusters  of  live  oak,  pinon,  and  juniper  at  the  lower  terminal  to  stands 
of  ponderosa  pine,  spruce,  fir,  and  aspen  between  the  lower  and  upper 
terminals.  In  addition,  many  varieties  of  wild  flowers,  shrubs,  grasses, 
and  succulent  plants  grow  on  the  west  side  of  the  mountain. 

Even  though  the  weather  in  the  Albuquerque  vicinity  is  usually 
clear  and  sunny,  there  are  many  days  throughout  the  year  when  Sandia 
Mountain  (including  the  tramway  facility)  is  enveloped  in  clouds.  While 
the  degree  of  cloudiness  does  not  appear  to  have  a  diurnal  variation, 
there  is  a  definite  seasonal  variation.  There  appears  to  be  a  predomin¬ 
ance  of  high  cumulus  type  rain  clouds  in  the  early  spring  and  late  summer, 
while  in  the  winter,  stratus  type  (including  fog)  clouds  appear  to  be  the 
most  numerous.  The  bases  of  the  spring  and  summer  clouds  are  usually 
higher  than  10,000  feet  while  in  the  winter  the  bases  are  generally 
much  lower,  varying  from  near  the  ground  (fog)  to  8000  feet. 

A  typical  view  of  one  of  the  tramway  cabins  and  the  surrounding 
terrain  on  a  partially  cloudy  morning  is  shown  in  Fig.  25. 

For  the  tramway  experiments  we  used  a  portable  battery-operated 
ladar  system  specifically  designed  for  use  in  the  field.  The  primary 
components  of  the  system  are  the  following: 

a.  Source;  GaAs  (gallium  arsenide)  semiconductor  laser  and  pulse 

driver  (IBM  Corp.) 

b.  Detector:  S-l  response  photomultiplier  tube  (RCA  Development 

type  C70102B) 

c.  Receiver:  solid  state  amplifier 

d .  Polarizers :  Clan-Thompson 

e.  Optics:  appropriate  mirrors,  lenses,  spectral  filters,  neutral 
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Figure  25.  Typical  View  of  Sandia  Tranway  Cabin  on  Partially  Cloudy  Day. 


density  filters,  etc. 

f.  Range  Indicator:  nixie  tube  digital  display  (XM23  Time  Interval 

Counter  Range  Finder,  IBM  Corp.) 

Additionally,  the  output  of  the  amplifier  was  displayed  on  a  Tektronix 
Model  454  oscilloscope.  A  Polaroid-Land  camera  was  used  for  photographing 
the  signal  returns. 

A  view  (inside  of  a  tram  cabin)  of  the  experimental  arrangement 
used  for  the  tramway  experiments  is  shown  in  Fig.  26.  A  complete  ladar 
system  including  the  XM23  nixie  tube  display  unit  is  mounted  on  the 
tripod.  A  similar  system  without  the  XM23  unit  is  mounted  on  a  tripod 
base  and  is  directed  through  an  opening  in  the  floor  of  the  tram  cabin. 

A  representative  sample  of  the  data  obtained  from  the  tramway 
experiments  is  shown  in  Figs.  27  through  29.  For  all  photographs  the 
tram  cabin  was  enveloped  in  clouds  with  the  type  of  cloud  varying 
from  very  dense  fog  to  very  dense  stratus.  The  visibility  varied  from 
10  to  20  feet.  A  representative  view  of  the  weather  conditions  is 
shown  by  the  photograph  in  Fig.  30.  The  ladar  unit  was  mounted  on  the 
tripod  base  located  at  the  opening  in  the  floor  of  the  tram  cabin,  and 
positioned  so  that  the  direction  of  propagation  of  the  transmitted  signal 
was  toward  the  snow-covered  terrain  below.  The  distance  between  trans¬ 
mitting  and  receiving  apertures  was  4  inches.  The  laser  was  operated  at  an 
output  power  (peak)  of  320  watts  at  a  pulse  repetition  frequency  of  20 
Hertz.  The  beamwidth  of  the  transmitter  and  the  field  of  view  of  the 
receiver  were  3  degrees.  A  300-Angstrom  spectral  filter  was  used  on  the 
receiver  to  reduce  the  back  ground  noise.  For  the  photographs  in  Figs. 

27a,  28a,  and  29a  the  plane  of  polarization  of  the  analyzer  (located  at 
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Figure  27.  Oscillogram  of  Cloud  and  Target  Returns,  (a)  Laser  Transmitter 
Output  Linearly  Folarize  Plane  of  Polarization  of  Analyzer 
(Located  at  Receiver)  Parallel  to  That  of  Transmitted  Signal; 

(b)  Laser  Transmitter  Output  Linearly  Polarized,  Plane  of  Pol¬ 
arization  of  Analyzer  (Located  at  Receiver)  Orthogonal  to  That 
of  Transmitted  Signal.  Vertical  Sensitivity:  .5v/div;  Sweep 
Rate:  100  nanosec/div.  Range:  Cloud  =  75  ft.,  Target  “  150  ft. 
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Figure  28.  Oscillogram  of  Cloud  and  Target  Returns,  (a)  Laser  Transmitter 
Output  Linearly  Polarized,  Plane  of  Polarization  of  Analyzer 
(Located  at  Receiver)  Parallel  to  That  of  Transmitted  Signal; 

(b)  Laser  Transmitter  Output  Linearly  Polarized,  Plane  of  Pol¬ 
arization  of  Analyzer  (Located  at  Receiver)  Orthogonal  to  That 
of  Transmitted  Signal.  Vertical  Sensitivity:  .2v/div;  Sweep 
Rate:  100  nanosec/div.  Range:  Cloud  -  100  ft.,  Target  s  175  ft. 
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Figure  29.  Oscillogram  of  Cloud  and  Target  Returns,  (a)  Laser  Transmitter 
Output  Linearly  Polarized,  Plane  of  Polarization  of  Analyzer 
(Located  at  Receiver)  Parallel  to  That  of  Transmitted  Signal; 

(b)  Laser  Transmitter  Output  Linearly  Polarized,  Plane  of  Pol¬ 
arization  of  Analyzer  (Located  at  Receiver)  Orthogonal  to  That 
of  Transmitted  Signal.  Vertical  Sensitivity:  .2v/div;  Sweep 
Rate:  200  nanosec/div.  Range:  Cloud  «  80  ft.,  Target  *  250  ft. 
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Figure  30.  Typical  View  of  Sandia  Tramway  Cabin  Enveloped  in  a  Dense 
Stratus  Cloud.  Camera  Was  Approximately  50  ft.  from  Cabin. 


the  receiver)  was  parallel  to  the  plane  of  polarization  of  the  trans¬ 
mitted  signal.  In  all  photographs  the  first  return  represents  the 
radiation  backscattered  from  the  cloud  while  the  second  return  represents 
the  radiation  reflected  by  the  terrain  (target).  The  difficulty  in 
discriminating  between  the  cloud  return  and  the  target  return  by  con¬ 
ventional  pulse  discrimination  techniques  is  obvious,  and  the  photo¬ 
graphs  clearly  illustrate  the  validity  of  the  polarization  discrimin¬ 
ation  concept.  In  addition,  a  comprehensive  evaluation  of  the  data 
revealed  that:  (a)  while  the  target  (terrain)  return  varied  in  range, 
the  cloud  return  remained  constant  (approximately)  in  range  and  appeared 
to  originate  at  a  distance  of  75(±  15)  feet  irrespective  of  the  density 
of  the  cloud,  and  (b)  the  terrain  return  appeared  to  be  depolarized 
for  all  varieties  of  vegetation  (whether  covered  with  snow  or  not) . 


SECTION  VII 


DISCUSSION 


7.1  SUMMARY  AND  CONCLUSIONS 

The  use  of  laser  detection  and  ranging  systems  in  the  atmosphere 
is  seriously  limited  by  the  severe  attenuation  of  the  signal  due  to 
absorption  by  atmospheric  gases  and/or  scattering  by  cloud  and  fog 
particles.  Furthermore,  the  cloud  and  fog  particles  may  scatter  the 
transmitted  radiation  in  the  direction  of  the  receiver,  thus  producing 
a  false  return.  This  false  return  must  he  distinguished  from  the  true 
target  return.  The  use  of  range-gating  techniques  to  eliminate  the  '  • 
cloud  returns  is  unacceptable  in  the  majority  of  applications  and  any 
electronic  discrimination  technique  based  on  determining  the  difference 
in  characteristics  between  the  cloud  and  target  returns  would  probably 
be  very  complicated,  if  not  impossible,  to  implement.  What  is  needed, 
then,  is  a  cloud/ target  discrimination  technique  that  while  being 
physically  realizable  is  also  reliable  and  inexpensive  to  implement. 

The  purpose  of  the  investigation  described  in  this  report  was  to 
determine  the  feasibility  of  using  an  optical  technique  for  discrim¬ 
inating  between  a  cloud  return  and  a  target  return.  Mie  [Ref.  42]  has 
shown  analytically  that  if  a  linearly  polarized,  monochromatic  electro¬ 
magnetic  plane  wave  propagating  through  a  homogeneous  medium  is  incident 
upon  a  single  spherical  particle  the  diameter  of  which  is  on  the  same 
order  of  magnitude  as  the  wavelength  of  the  incident  radiation,  then 
the  scattered  radiation  is,  in  general,  either  linearly  or  elliptically 
polarized  depending  on  the  plane  of  observation.  Other  investigators 
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have  extended  (but  not  analytically)  the  Mie  theory  to  a  homogeneous 
medium  (such  as  atmospheric  cloud)  containing  a  size  distribution  of 
independently  scattering  particles.  Beckmann  and  Spizzichino  [Ref.  4], 
on  the  other  hand,  have  shown  that  if  a  linearly  polarized,  monochromatic 
electromagnetic  plane  wave  is  incident  on  an  arbitrarily  inclined  rough 
surface,  then  the  scattered  radiation  will,  in  general,  be  depolarized 
if  the  radius  of  curvature  of  the  surface  irregularities  is  large  com¬ 
pared  to  the  wavelength  of  the  incident  radiation.  Since  the  major? ty 
of  useful  lasers  emit  radiation  in  the  visible  or  near  infrared  portion 
of  the  electromagnetic  spectrum,  atmospheric  clouds  or  fogs  actually 
represent  Mie  scattering  media  and  most  targets — earth  terrain  including 
man-made  irregularities,  airborne  objects  such  as  aircraft,  etc., — can 
be  considered  to  have  rough  surfaces  at  these  wavelengths.  One  possible 
method,  then,  of  distinguishing  a  cloud  return  from  a  target  return 
could  be  based  on  the  difference  in  the  polarization  characteristics 
of  each  return.  Specifically,  if  the  transmitter  output  of  a  laser 
detection  and  ranging  system  is  linearly  polarized  and  a  polarizer 
(analyzer)  is  placed  at  the  input  to  the  receiver  (assumed  to  be  located 
adjacent  to  the  transmitter)  with  its  plane  of  polarization  orthogonal 
to  that  of  the  transmitted  wave,  the  receiver  will  "see"  only  the  com¬ 
ponent  of  the  radiation  reflected  from  the  target  whose  electric  vector 
vibrates  at  right  angles  to  the  electric  vector  of  the  transmitted  wave. 
All  other  returns — due  to  scattering  from  cloud  particles  and  the  paral¬ 
lel  polarized  component  of  the  radiation  reflected  from  the  target — will 
be  eliminated  by  the  analyzer. 
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The  results  of  the  theoretical  and  experimental  investigation 
described  in  the  previous  chapters  verify  conclusively  the  validity  of 
the  cloud/target  polarization  discrimination  concept.  In  the  laboratory 
phase  of  the  program  a  linearly  polarized  signal  from  a  He-Ne  gas  laser 
was  propagated  through  an  artifically  produced  cloud.  The  total  inten¬ 
sity  and  the  corresponding  intensities  of  the  parallel  and  perpendicular 
electric  field  components  of  the  radiation  scattered  out  of  the  laser 
beam  by  the  cloud  particles  were  measured  as  a  function  of  the  transmis¬ 
sivity,  t,  of  the  cloud.  The  resultant  cloud  return  data  was  normalized 
to  a  clear  air  reference  target  return.  Plots  of  log  I  vs.  log  t  for  the 
cloud  returns  and  their  corresponding  target  returns  are  shown  in  Figs. 

20  through  23.  These  curves  show  that  for  t  >  0.75  the  entire  cloud- 
scattered  radiation  is  polarized  the  same  as  the  transmitted  signal 
(i.e.  I  *  I  ) ,  and  that  for  t  <  0.75  almost  all  of  the  cloud  return 
is  pdarized  the  same  as  the  transmitter  except  for  a  small  portion  that 
is  cross  polarized  relative  to  the  transmitted  signal  (i.e.  I  -  I£+) . 
These  results  indicate  that  for  t  >  0.75  single  scattering  prevailed 
while  for  t  <  0.75  the  backscattered  radiation  was  partially  depolarized 
due  to  multiple  scattering.  Thus,  these  curves  verify  that  for  less 
dense  clouds  the  cloud-scattered  radiation  is  almost  in  exact  agreement 
with  the  Mie  Theory  (all  of  the  backscatter  is  polarized  the  same  as  the 
transmitted  wave);  whereas,  for  more  dense  clouds  there  is  some  back- 
scattered  radiation  cross  polarized  relative  to  the  transmitter  polari¬ 
zation.  This  partial  depolarization  of  the  cloud  return  at  the  higher 
cloud  densities  was  caused  by  an  increase  in  the  multiple  scattering. 

In  Figs.  22  and  23  we  see  that  while  the  intensity  of  the  cross-polarized 
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component  (I  )  of  the  backscattered  radiation  should  be  zero,  it  in 
fact  increases  as  the  density  of  the  cloud  increases  indicating  the 
presence  of  depolarization  effects  due  to  multiple  scattering.  Addi¬ 
tional  evidence  supporting  our  conclusion  that  the  partial  depolarization 
of  the  cloud  returns  was  caused  by  multiple  scattering  is  given  by  van 
de  Hulst  [Ref.  54],  He  states  that  for  an  atmospheric  cloud  single 
scattering  prevails  for  less  dense  clouds  (t  >  0.9);  whereas,  the  effects 
of  multiple  scattering  may  be  observed  in  the  range  0.74  <  t  <  0.90, 
multiple  scattering  definitely  occurs  in  the  more  dense  clouds  when 
t  <  0.74.  We  conclude,  then,  that  the  laboratory  produced  cloud  actually 
represents  a  Mie  scattering  media  whose  scattering  characteristics 
simulate  a  naturally  occurring  atmospheric  cloud.  In  addition,  the 
size  distribution  of  cloud  particles,  as  determined  from  photomicro¬ 
graphs,  indicates  that  the  laboratory  produced  cloud  simulates  a  low 
altitude  atmospheric  stratus  cloud. 

We  see  also  from  Figs.  20,  21,  and  23  that,  when  no  analyzer  was 

used  on  the  receiver  or  when  the  plane  of  polarization  of  the  analyzer 

was  parallel  to  that  of  the  transmitted  signal,  the  cloud  return  was 

greater  than  the  target  return  for  t  <  0.5  (i.e.  I.  >  It+  and  Ic))  > 

I  for  t  <  0.5).  On  the  other  hand,  in  Figs.  22  and  23  we  see  that 
t  it 

when  the  plane  of  polarization  of  the  analyzer  was  orthogonal  to  that  of 

the  transmitted  signal,  the  cloud  return  was  reduced  approximately  94 

percent  and  was  greater  than  the  target  return  only  when  t  <  0.28  (i.e. 

I  >  I  for  t  <  0.28).  The  data  of  Arnulf,  et  al,  [Ref.  3],  Carrier, 
cl  tl 

et  al,  [Ref.  11],  and  other  researchers  indicate  that  the  transmissivity 
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per  meter  of  most  clouds  and  fogs  is  usually  between  70  to  90  percent. 
While  the  data,  then,  in  Figs.  20  and  23  indicate  that  for  transmis¬ 
sivities  above  48  percent  the  target  return  will  be  greater  than  the 
cloud  return,  even  when  the  analyzer  is  not  used  on  the  receiver,  the 
target  in  this  case  (a  block  of  MgO)  was  almost  both  a  perfect  reflector 
and  a  perfect  diffuser.  Most  natural  earth  objects  (water  surfaces, 
coniferous  forest,  grass  meadows,  plowed  fields,  etc.)  are  very  nearly 
perfectly  diffuse  targets  at  laser  wavelengths;  however,  when  compared 
to  MgO  they  usually  are  poor  reflectors  (except  possibly  when  covered 
with  snow).  The  reflection  coefficients  of  various  earth  surfaces  is 
tabulated  in  the  Smithsonian  Meteorological  Tables  [Ref.  37]  and  the 
values  vary  from  a  low  of  0.03  -  0.10  for  water  surfaces  to  a  high  of 
0.03  -  0.25  for  ground  surfaces.  The  reflectance  of  fresh  fallen  snow 
varies  from  0.70  to  0.86.  The  coefficient  of  reflection  for  man-made 
objects  (roads,  buildings,  etc.)  varies  from  0.08  to  0.35  [Ref.  51]. 

The  value  of  the  discrimination  technique  when  used  in  a  realistic 
situation  is  illustrated  in  Fig.  31.  Here  we  replotted  the  cloud 
returns,  Ic+  and  I  ,  from  Fig.  23  substituting  for  the  MgO  returns, 
two  representative  earth  targets  (p  *»  0.10  and  p  «  0.25).  The  earth 
targets  are  assumed  perfect  diffusers  and  their  returns  are  normalized 
to  the  reflectance  of  the  MgO  target.  These  curves  show  that  when  no 
analyzer  is  used  on  the  receiver,  the  cloud  return  will  be  greater  than 
the  return  from  a  typically  good  reflecting  target  (p  -  0.25)  for 
t  <  0.65  (i.e.  I£+  >  I  for  x  <  0.65).  Similarly,  for  an  average 
reflecting  target  (p  -  0.10)  the  cloud  return  will  be  greater  than  the 
target  return  (i.e.  I  >  I')  for  x  <  0.75.  On  the  other  hand,  if  a 

C*r  lt 
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Figure  31.  Target  and  Cloud  Returns  vs.  Cloud  Transmissivity  for  Raallstic 
Earth  Targets 
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polarizer  is  placed  at  the  input  to  the  receiver  with  its  plane  of  polar¬ 
ization  orthogonal  to  that  of  the  transmitted  signal,  we  see  that  the 
two  target  returns  (I^.and  1'^)  are  greater  than  the  cloud  return  (1  ) 

for  t  >  0.36  and  t  >  0.48,  respectively.  These  transmissivity  levels 
correspond  to  very  dense  clouds  which  do  not  occur  very  frequently  in 
the  atmosphere.  We  conclude,  then,  that  the  use  of  the  polarization 
technique  will  greatly  improve  the  ability  of  a  laser  detection  and 
ranging  system  to  discriminate  between  the  returns  from  an  atmospheric 
cloud  and  a  realistic  target. 

The  data  collected  in  the  field  evaluation  phase  of  the  experi¬ 
mental  program  clearly  confirms  the  results  of  the  laboratory  investi¬ 
gation.  The  field  tests  were  performed  in  one  of  the  two  cabins  of  the 
Sandia  Peak  Aerial  Tramway.  The  cabin  was  used  as  a  transmitter/receiver 
site  and  the  terrain  below  the  cabin  was  used  as  a  target.  The 
tramway  experiments  were  conducted  with  a  portable  battery-operated  GaAs 
(gallium  arsenide)  rangefinder  system  specifically  designed  for  use  in 
the  field.  A  representative  sample  of  the  data  obtained  in  the  field 
evaluation  phase  is  presented  in  Figs.  27  through  29.  In  Fig.  27,  the 
relative  intensity  of  the  cloud  return  was  reduced  80  percent  when  the 
plane  of  polarization  of  the  analyzer  was  rotated  90  degrees;  that  is, 
the  oscilloscope  display  of  the  cloud  is  reduced  from  I  «  lv.  (peak) 

C  I! 

in  Fig.  27a,  to  I  ■  0.2v.  (peak)  in  Fig.  27b.  Similarly,  in  Figs.  28 
and  29,  the  relative  intensities  of  the  cloud  returns  were  reduced  75 
and  90  percent  respectively.  The  average  reduction  in  the  relative 
intensity  of  the  cloud  returns  for  all  of  the  tramway  data  ranged  between 
75  to  90  percent.  In  addition,  a  comprehensive  evaluation  of  the  data 
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revealed  that;  (a)  while  the  target  (terrain)  return  varied  in  range, 
the  cloud  return  remained  constant  in  range  (approximately)  and  appeared 
to  originate  at  a  distance  of  75  (±15)  feet  irrespective  of  the  density 
of  the  cloud,  and  (b)  the  terrain  return  appeared  to  be  depolarized  for 
all  varieties  of  vegetation  (whether  covered  with  snow  or  not) >  The 
expected  signal  returns  from  the  cloud  and  the  terrain  below  the  tram 
car  were  calculated  from  Equation  (285) 

pP,_A_e”2cfrcosY 

P  „  __L_5 -  (285) 


and  Equation  (301) 


where 


6(tt)ArIo2tt(1  -  cos  ~)E2(20) 


and 


r  -2crx 

E2(2o)  ”  /  - - dx 

r  x2 
o 

The  receiver  output  voltage,  V  ,  was  calculated  from 

K 

VR  "  NR  PR  S  ^  \ 


(301) 

(281) 


(300) 


(309) 


Where  N  is  the  efficiency  of  the  receiver  optics,  S  is  the  sensitivity 
R 

of  the  detector,  is  the  load  resistance,  and  is  the  gain  of  the 
receiver.  The  calculated  values  for  the  cloud  and  target  returns  were 
in  good  agreement  with  the  measured  values  for  an  assumed  terrain  (target) 
reflectivity  p  •  0.15. 
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The  following  calculations  are  presented  to  illustrate  the  close 
correlation  between  the  measured  and  calculated  values  of  the  cloud  and 
target  returns.  Specifically,  the  measured  values  of  the  receiver  output 
voltage  for  the  cloud  and  target  returns  shown  in  Fig.  29  are  compared 
to  the  calculated  values  where 


P  " 

PT  “ 

ar“ 

o  - 


0.15 

320  watts 

2 


144 


ft.' 


0.015m  1  *  .0045  ft. 


y  -  45° 
r  -  250  ft. 


r  -  2  ft. 
o 


0T  *  60  rad’ 


eR  “  60  rad‘ 


6(h)  *  6  x  10  4m  X-ster  1 


3.84  x  10-7  ft."1 


E2(2o)  =0.49 


N_  -  0.0315 

K 


3100 


amps 

watt 


-  316  n 
Ay  “  368 

Substituting  these  values  into  the  appropriate  equations,  we  obtain 
the  following  values  for  the  target  return. 
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PR  -  3.77  x  10" '  watts 

V  «  0.4253  volts 
iR 

Similarly,  for  the  cloud, 

P„  -  1.3  x  10  b  watts 

K 

*  1.47  volts 

The  values  of  the  receiver  output  voltages  measured  from  Fig,  29 
are 

VR  -  0.4  volts  (peak) 

VD  =  1.0  volt  (peak) 

K 

The  calculated  value  of  the  receiver  output  voltage  for  the  target 
return  is  in  good  agreement  with  the  measured  value.  The  measured  value 
of  the  receiver  output  voltage  for  the  cloud  return  is,  however,  almost 
50  percent  greater  than  the  calculated  value.  The  reason  for  this  differ¬ 
ence  i s  that  the  detector  was  probably  saturated  by  the  cloud  return. 

The  data  collected  from  both  the  laboratory  verification  phase 
and  the  field  evaluation  phase  of  the  program  has  definitely  shown  that 
the  denominator  of  the  polarization  discrimination  coefficient 


PDC  « 


800 

Bj.00 


(307) 


can  be  physically  made  very  small.  This  substantiates  the  assumption 
that  the  Mie  theory  of  scattering  (originally  derived  for  a  single 
sphere)  is  applicable  to  a  homogeneous  medium  consisting  of  many  scattering 
particles,  not  necessarily  of  the  same  size,  such  as  the  clouds  that  were 
generated  in  the  laboratory  and  that  were  encountered  in  the  field  tests 
on  board  the  tram. 
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7.2  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

(a)  Mie  presented  his  analytical  solution  for  the  field  components 
of  the  radiation  scattered  from  a  linearly  polarized  monochromatic  plane 
wave  incident  on  a  single  homogeneous  spherical  particle  in  1908.  Since 
then  many  reserachers  have  extended  the  Mie  theory  to  a  homogeneous 
medium  containing  a  size  distribution  of  independently  scattering  par¬ 
ticles.  These  expanded  versions  of  the  Mie  theory  are,  in  reality,  only 
approximations  based  on  simplifying  assumptions  that  may  not  be  valid. 

The  ultimate  use  of  the  laser  in  present  and  proposed  applications  where 
the  laser  energy  may  be  severely  attenuated  when  propagating  through  a 
medium  containing  many  scatterers  will  depend  on  the  results  of  an  accur¬ 
ate  analysis  of  the  interaction  between  the  propagating  wave  and  the 
scattering  medium.  An  analytical  solution  based  on  electromagnetic  theory 
is  required,  then,  for  the  diffraction  of  a  collimated  beam  of  electro¬ 
magnetic  energy  propagating  through  a  medium  consisting  of  many  spherical 
(or  similar  shaped)  particles,  not  necessarily  of  the  same  diameters. 
Furthermore,  the  effects  of  multiple  scattering  definitely  should  be 
included  in  the  analysis. 

(b)  While  the  validity  of  the  discrimination  technique  has  been 
experimentally  verified  for  two  optical  wavelengths  (one  in  the  visible 
and  one  in  the  infrared  portion  of  the  electromagnetic  spectrum),  a 
thorough  analytical  and  experimental  analysis  of  the  wavelength  depen- 
dance  (if  any)  of  the  concept  has  not  been  performed.  It  has  been  shown 
by  several  researchers  [Refs.  3,  11,  13]  that  the  transmissivity  of 
clouds  and  fogs  in  virtually  independent  of  wavelength.  While,  on  the 
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other  hand,  the  transmissivity  of  haze  Increases  rapidly  with  increasing 
wavelength  until  at  10  microns  the  haze  becomes  almost  transparent 
[Refs.  3,  13].  Deirmendjian  [Ref.  13]  and  Carrier,  et  al  [Ref.  11]  have 
shown  that  the  backscattered  radiation  is  definitely  wavelength  sensitive 
with  the  amount  of  backscatter  decreasing  as  a  function  of  increasing 
wavelength. 

(c)  Arnulf,  et  al  [Ref.  3]  has  reported  that  the  optical  properties 
of  artificial  smoke  are  very  similar  to  those  of  haze,  that  is,  the 
transmissivity  of  artificial  smoke  Increases  with  increasing  wavelength 
and  approaches  100  percent  at  10  microns.  Additionally,  physical 
measurements  of  the  composition  of  artificial  smoke  demonstrated  that  it 
consists  of  spherical  droplets  and  almost  never  solid  particles.  The 
results  of  recent  artificial  smoke  experiments  conducted  in  our  labora¬ 
tory  are  in  good  agreement  with  Arnulf's  findings  and  indicate  that  the 
polarization  discrimination  technique  works  equally  well  for  smoke  clouds. 
The  theoretical  and  experimental  program  should  be  expanded  to  include 
an  evaluation  of  the  effectiveness  of  the  discrimination  technique  *n 
other  aerosols  such  as  smog  and  smoke. 
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